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THE  SEVEN-HOLE  PRESSURE  PROBE 

Lawrence  Reed*,  Jack  D.  Mattingly**,  and  Frederick  Jonas** 

Abstract 

This  paper  documents  recent  and  past  developments  with  respect  to 
the  seven-hole  pressure  probe.  Included  are  discussions  on  probe 
design,  construction,  calibration,  and  measurement  capabilities.  The 
effects  on  probe  measurements  in  shear  flows,  as  well  as  methods  of 
correction,  are  also  included.  Finally,  research  applications  in  the 
measurement  of  unknown  flowfields  are  presented.  As  shown,  the 
seven-hole  pressure  probe  is  a  valuable  and  highly  accurate  device  for 
quantitatively  documenting  unknown  flowfields. 

I.  Introduction 


The  purpose  of  this  paper  is  to  document  recent  and  past 
developments  with  respect  to  the  seven-hole  pressure  probe  developed  at 
the  United  States  Air  Force  Academy  in  conjunction  with  the  NASA  Ames 
Research  Center.  As  with  all  pressure  probes,  the  seven-hole  probe  can 
be  used  to  map  unknown  internal  or  external  flowfields,  giving  valuable 
information  to  the  aerodynami c i s t .  The  advantage  of  the  seven-hole 
probe  over  other  flowfield  measuring  devices  is  its  widely  expanded 
measurement  range  and  flexibility.  To  document  this  effort,  Section  II 
presents  brief  background  material  about  the  probe's  design  and  use. 
Section  III  includes  an  overview  of  probe  calibration  in  uniform  flows, 
both  compressible  and  incompressible.  In  Section  IV,  recent  results  of 
investigations  into  probe  measurement  capablities  in  shear  flow 
environments  are  presented.  Measurement  capabilities  as  well  as 
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measurement  sources  of  error  are  discussed  in  Section  V.  Finally, 
recent  research  efforts  on  the  application  of  seven-hole  probe 
measurements  in  unknown  flows  are  presented  in  Section  VI. 

I  I .  Background 

Numerous  techniques  are  available  to  measure  unknown  flow  fields. 
While  methods  such  as  tufts,  streamers,  and  vanes  are  primarily  used 
for  flow  visualization,  they  are  insufficient  for  quantitative 
information.  Detailed  data  on  flow  size,  direction,  and  pressure 
usually  require  direct  obtrusive  flow  measurement. 

One  of  the  oldest  known  quantitative  techniques  may  be  found  in 
the  pressure  probe.  Although  other  techniques,  such  as  hot  wire 
anemometers  and  laser  doppler  anemometers,  have  been  developed, 
pressure  probes  are  desirable  from  a  flight  vehicle  standpoint  because 
of  their  simplicity  and  ruggedness  (Ref.  I).  Coupled  with  relatively 
low  cost,  pressure  probes  are  an  excellent  measurement  device  for 
research,  development,  and  industrial  applications. 

One  of  the  classic  problems  in  obtrusive  flow  measurement  is  the 
disturbance  caused  by  a  probe  to  the  flowfield  it  is  measuring.  For 
certain  applications,  such  as  aircraft  pitot-static  tubes,  this 
disturbance  or  perturbation  is  relatively  unimportant.  These  probes 
may  be  large,  for  the  relative  perturbation  of  the  probe  to  that  of  the 
aircraft  is  negligible;  further,  installation  of  the  probes  is 
accounted  for  in  the  airspeed  measurement.  To  determine  the  effect  of 
flow  disturbances  for  other  applications  requires  varying  levels  of 
knowledge  about  the  flowfield  being  measured  as  well  as  probe  geometry 
and  installation.  The  accuracy  and  measurement  ranges  of  such  probes 
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also  depend  on  the  size,  geometry,  number  of  holes,  and  installation. 
The  seven-hole  probe  was  developed  at  the  United  States  Air  Force 
Academy  to  increase  the  measurement  range  of  an  non-nulling  pressure 
probe  while  minimizing  probe  size. 

A.  Nulling  Versus  Stationary  Probes 

The  increased  measurement  range  of  a  multi-hole  probe  revolves 
on  its  ability  to  sense  more  pressures  at  the  face  of  the  probe.  Two 
procedures,  nulling  and  stationary,  may  be  used  to  acquire  the  pressure 
data.  A  nulling  probe  is  rotated  in  one  or  two  planes  until  opposing 
peripheral  ports  measure  equal  pressures.  The  corresponding  angle  of 
rotation  determines  the  flow  angle.  For  the  non-nulling  approach,  the 
probe  is  held  stationary  as  pressures  in  opposing  peripheral  ports  are 
recorded.  The  pressure  differences  are  then  transformed  to  flow  angles 
through  previously  known  calibration  relationships.  Although  the 
nulling  procedure  allows  analysis  of  high  angle  flow,  it  is 
mechanically  complex,  time  consuming,  and  hence,  may  not  be  cost 
effective.  On  the  other  hand,  stationary  probes  with  up  to  five  holes 
are  incapable  of  accurate  measurements  at  high  local  flow  angles 
(greater  than  30  degrees)  because  of  flow  separation  around  the  probe 
tip.  Explained  in  detail  in  a  later  section,  the  seven  hole  probe  is 
the  only  non-nulling  probe  that  can  determine  flow  angles  up  to  70 
degrees  relative  to  its  axis.  When  combined  with  a  computerized  data 
acquisition  system,  the  probe  may  record  nearly  two  data  points  per 
second.  This  rate  is  significantly  faster  than  current  nulling 
devices,  which  require  considerable  time  for  the  balancing  of  probe  tip 
pressures  before  each  measurement  can  be  taken  (Ref.  2). 
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B.  Analytical  Model  of  Pressure  Probe 

Pressure  probes  used  for  measurement  of  flow  direction  and 
magnitude  normally  consist  of  an  aerodynamic  body  with  a  symmetrical 
arrangement  of  sensing  holes  (Figure  1).  Some  typical  geometries  for 
pressure  probes  are  reviewed  in  References  3  through  17  and  pressure 
probes  are  treated  in  more  general  terms  in  References  18  and  19. 


TOTAL  PRESSURE 
PORT 


Figurs  1.  TYPICAL  FLOW-DIRECTION  PROBE 

(EACH  UNIT  ON  SCALE  REPRESENTS  0.6mm) 

Since  the  pressure  probes  of  interest  are  to  be  used  without 
rotation,  the  sensitivity  to  flow  angularity  is  extremely  important.  A 
general  analytical  model  was  developed  by  Huffman  using  slender  body 
theory  (Ref.  20),  which  permits  synthesizing  a  sensor  geometry  for  a 
special  application  and  analysis  of  probe  behavior. 

While  numerically  solving  the  three-dimensional  potential  flow 
equations  for  the  flow  around  an  aerodynamic  probe  may  be  more  accurate 
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And  when  evaluated  on  the  body  the  following  expression  results: 

( 

cp  -  [  f  +  f  2/A+(R')2]  cos2a  c  o  s 2  B  +sin2B  [l-4sin20  -4(R')2  cos20] 

+  sin2a  c o 8 2  B  [ l-4cos 2 0- 4 ( R ' ) 2 s in  2 0 ] +2 s ina  sin2B  sin2Q  [ 1 - ( R * ) 2 ] 
-2R'(l+f/2)  [cosa  sin2B  cos0+sin2a  cos2  8  sinO] 

where  f  is  to  be  evaluated  on  the  body  surface,  i.e.  r=R.  Equation  4 
represents  the  major  result  of  Huffman's  analysis  and  can  be  used  to 
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determine  probe  angular  and  static  pressure  sensitivity  for  arbitrary 
geometr ies . 

The  pressure  difference  between  ports  on  opposite  sides  of  the 
probe  is  used  to  determine  the  flow  angularity.  Since  the  side  ports 
of  the  probe  are  at  the  same  z  and  r  locations,  the  pressure  difference 
is  proportional  to  the  difference  of  Cp  values  at  different  e  locations 
and 


CP  “Cd  “  4s1,,28  [8in20.-8in20.+(R')2(cos20-cos20.)l 

i  pj  J  1  j  i  < 

+4sin2o  c os 2 B  [cos20,-cos20  +(R* )2 (sln20  -sin20  ) ] 

J  J  1 

+2sina  sin2$  [l-(R')2]  [ sin201-sin20 ^ ] 

-2R'  (l  +  f/  2)  [  coso  sin26  (cosOj^-cos©  )+8in2a  co82B(sin0  .-sin©  )] 

J  1  J 


Figure  3.  SEVEN-HOLE  PROBE 
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The  average  pressure  coefficient  can  be  used  to  determine  the 
flowfield  static  pressure,  P,, .  Cp  is  evaluated  at  a  number  of  theta 
values  and  the  results  summed  and  divided  by  the  number  of  theta 
values.  For  the  seven-hole  probe  of  Figure  3,  this  process  yields 


[  1-f-f 2/4+(R')2J  cos 2 a 


cob  2  (5-  { 1+2  (R  ’ )  2  ] 


(10) 


Huffman  developed  an  expression  for  a  quasi  total  pressure  from 
his  analysis  by  integrating  the  pressure  from  the  body  surface, 
resolving  this  force  into  an  axial  component,  and  dividing  by  the  cross 
sectional  area.  His  relationship  is 

l 

(cp)z-  C— -  /(R2)'[f+f2/4+(R*)  ]dz-slnWn2a  cos2P  <U> 

*o 

where  l  denotes  the  integration  length. 

For  a  conical  shaped  body,  R*  is  constant  and  Equation  11  reduces 
to 


[  1-f-f  2/4-(R')2]  cos2a  cos20-l 


(12) 


A  probe  pressure  coefficient  that  is  proportional  to  the  dynamic 
pressure  can  be  obtained  by  subtracting  Equation  10  from  Equation  12 
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yielding  a  dynamic  pressure  coefficient. 


(Cp)z"(Cp)"  2^R' ) 2 [ l-cos2a  coa2ej 


(13) 


Some  insight  into  probe  calibration  can  be  obtained  by  studying 


the  plot  of  a  and  0  versus  ACu  and  ACe  in  Figure  5,  the  plot  of  <Cp> 


versus  a  and  B  in  Figure  6,  the  plot  of  <Cp>z  versus  a  and  B  in  Figure  7, 


and  the  plot  of  £<C  >z-<C  >] versus  a  and  (5  in  Figure  8.  Figures  1  through 


P  P 

8  were  calculated  for  R'=.268  and  f=0.  Note  the  near  linear  relation 


between  each  angle  and  its  respective  pressure  coefficient  as  shown  in 


Figure  5.  Also,  the  nearly  circular  contours  of  constant  <C  >,  <C  >z,  and 

P  P 


<Cp>z-<Cp>^Jas  shown  in  Figures  6,7,  and  8,  respectively.  Probe 


P  ~P 

manufacturing  imperfection  and  angular  misalignment  of  the  pressure 
ports  will  cause  these  curves  to  become  somewhat  distorted. 
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Figure  6.  CONTOURS  OF  CONSTANT<Cp> 
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CONTOURS  OP  CONSTANT  <C  p  >, 
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C.  Seven-Hole  Probe  Design 

The  seven-hole  probe  is  characterized  by  six  periphery  ports 
surrounding  one  central  port  (Fig.  9).  The  probe  is  constructed  by 
packing  seven  properly  sized  stainless  steel  tubes  into  a  stainless 
steel  casing.  For  the  current  probe  used  at  the  Air  Force  Academy,  the 
inner  seven  tubes  have  an  outside  diameter  of  .028  inches  with  a 
.0005-inch  wall  thickness.  Once  assembled  in  the  order  shown  below, 
the  tubes  are  soldered  together  and  machined  to  the  desired  half  angle 
(usually  25  or  30  degrees)  at  the  tip  (Ref.  2).  It  should  be  noted 
that  the  manufacture  of  seven-hole  probes  (versus  four/five  hole 
probes)  is  much  simpler  due  to  this  packing  arrangement. 


STAINLESS  TUBING 
K  SOLDER 


& 


0.109  IN. 
DIAM. 


Figure  9.  PROBE  GEOMETRY 
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III.  Seven-Hole  Probe  Calibration  Theory 

Because  of  their  small  size  and  individual  construction,  all 
probes  have  inherent  manufacturing  defects  that  require  the  unique 
calibration  of  each  probe.  Gallington  (Ref.  24)  developed  the 
calibration  theory  required  for  incompressible,  uniform  flows.  His 
power  series  method  produces  explicit  polynominal  expressions  for  the 
desired  aerodynamic  properties  and  is  easily  programmed.  The  following 
section  presents  a  synopsis  of  Gallington's  scheme  for  incompressible, 
uniform  flow  calibration  as  presented  by  Gallington  (Ref.  20)  and  by 
Gerner  and  Mauer  (Ref.  2). 

A.  Low  Versus  High  Flow  Angles,  Incompressible  Flow 

In  this  calibration  technique  for  incompressible,  uniform 
flow,  the  probe  face  is  divided  into  two  sectors.  The  inner  flow 
sector  deals  with  low  flow  angles  in  which  the  angle  between  the 
probe's  axis  and  the  freestream  velocity  vector  is  less  than  30 
degrees.  The  other  or  outer  sector  deals  with  high  flow  angles  of  over 
30  degrees.  Thirty  degrees  is  the  dividing  point  because  flow 
typically  begins  to  detach  over  the  top  surface  of  the  probe  tip  around 
this  local  flow  angle.  The  hole  numbering  system  used  for  both  sectors 
and  the  remainder  of  the  report  is  as  follows: 


FRONT  VIEW 


SIDE  VIEW 


Figure  10.  PORT  NUMBERING  CONVENTION  AND  PRINCIPAL  AXES 
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As  the  numerator  measures  changes  in  flow  angularity,  the  denominator 
nond imens i ona 1 i zes  the  expression  with  the  apparent  dynamic  pressure. 
The  center  port  pressure,  P^,  approximates  the  local  total  pressure 
while  the  average  of  the  circumferential  port  pressures  Pj  ^ 
approximates  the  local  static  pressure.  To  transform  these  pressure 
coefficients  to  the  tangential  reference  system,  Gallington  (Ref.  24) 
formulates  the  following  relationships: 


C  =  ^  ( 2  C  +  C  -  C  ) 
a  3  a l  a2  as 


1 


C.  =  (C  +  C  ) 

6  yy  a  2  U  3 


(15) 


It  is  important  to  realize  that  and  C&  are  not  independent  of  each 
other;  that  is,  Cu  is  a  function  of  all  six  peripheral  ports  while  C6 
is  a  function  of  all  but  ports  1  and  4. 

Besides  the  two  angular  pressure  coefficients,  two  other  low  angle 


pressure  coefficients,  CQ  and  Cq,  are  defined  in  Gallington's  work 


(Ref.  24) 


C 

o 


P?  -  P  1  ~  6 


c 

q 


P7  -  Pl “6 


oL 


P 


“-L 


(16) 


CQ  is  the  apparent  total  pressure  coefficient  with  respect  to  each  hole 

and  is  a  means  to  convert  actual  pressures  measured  by  the  probe  to 

accurate  values  of  local  total  pressure.  The  numerator  measures  the 

difference  between  the  approximate  total  pressure  measured  by  the 

center  port  And  the  actual  local  total  pressure  PQ^.  As  with  the 

previous  coefficients,  the  denominator  nond imens i ona 1 i zes  the 

expression  with  the  apparent  local  dynamic  pressure. 

The  velocity  pressure  coefficient,  Cq ,  serves  a  similar  conversion 

function  as  C  except  that  it  relates  the  probe  pressures  to  the  actual 
o 

dynamic  pressure.  The  numerator  in  this  coefficient  represents  the 
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probe's  approximation  of  the  local  dynamic  pressure  while  the 
denominator  represents  the  actual  dynamic  pressure  of  the  freestream 
test  cond i t i ons . 

0.  Probe  Axis  System  for  High  Flow  Angles 

The  real  advantage  of  using  the  non-nulling  seven-hole  probe 
over  other  multi-hole  probes  appears  in  the  ability  to  measure 
high-angle  flows.  At  high  angles  of  attack  (greater  than  30  degrees), 
the  flow  detaches  over  the  upper  surface  of  the  probe,  and  pressure 
ports  in  the  separated  wake  are  insensitive  to  small  changes  in  flow 
angularity  (Fig.  12). 


Figure  12.  FLOW  PATTERN  OVER  SEVEN-HOLE  PROBE 
AT  HIGH  ANGLE  OF  ATTACK 
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Consequently,  high  angle  measurements  must  only  be  made  from  ports  in 
the  attached  flow.  Five-hole  and  other  pressure  probes  have  been 
ineffective  in  this  regime  due  to  this  flow  separation  and  lack  of 
sensing  ports  in  the  attached  flow  region.  The  seven-hole  probe  avoids 
this  problem  because  at  least  three  ports  remain  in  the  attached  flow 
region,  allowing  sufficient  data  to  be  recorded  to  document  the  flow 
angle. 

According  to  Gallington,  the  tangential  reference  system  is 
inappropriate  at  high  flow  angles  because  of  indeterminate  angles  and 
singularities  (Ref.  24).  Instead,  the  polar  reference  system  is  used 
(Fig.  13)  where  e  represents  the  pitch  angle  and  $  represents  the  roll 
angle. 


POLAR 

TANGENT 

u  =  V  cos  6 

v  =  V  sln6  sin  $ 

»*V  sine  cos  $ 

a_  =  arctan  1! 

1  u 

6-  =arctan  1 

1  u 

Figure  13.  HIGH  ANGLE  REFERENCE  SYSTEM 

More  specifically,  e  denotes  the  angle  that  the  velocity  vector  forms 
with  respect  to  the  probe's  x-axis,  and  *  signifies  the  azimuthal 
orientation  of  the  velocity  vector  in  the  y-z  plane.  We  should  note 
that  a  positive  $  is  measured  counterclockwise  from  the  negative  z  axis 
when  the  probe  is  viewed  from  the  front. 
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Based  on  Gallington's  original  work  with  the  seven-hole  probe  at 
high  angles  of  attack,  Gerner  and  Mauer  accomplished  the  following 
analysis  (Ref.  2).  Kuethe  and  Chow's  Foundations  of  Aerodynamics: 
Bases  of  Aerodynamic  Design,  explains  that  the  separation  points  of  a 
cylinder  in  turbulent  flow  are  over  100  degrees  from  the  frontal 
stagnation  point.  Flow  around  a  conical  body  such  as  a  probe  tip  was 
likely  to  remain  attached  longer,  and  the  u-velocity  component  was 
prone  to  extend  the  separation  points  downstream  even  further. 
Therefore,  for  the  high  angle  flow  shown  in  Figure  14,  ports  3,  4,  and 
7  lie  in  attached  flow,  and  port  1  lies  in  separated  flow.  The  flow 
over  ports  2  and  6  is  unpredictable,  and  their  readings  are  discarded 
(Ref.  2). 


Flour*  14.  FLOW  OVER  PROBE  AT  HIGH  ANGLE  OF  ATTACK 

E.  Pressure  Coefficients  for  High  Flow  Angles 

As  with  low  flow  angles,  high  flow  angularity  is  measured  by 
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dimensionless  pressure  coefficients  that  are  functions  of  port  pressure 
differences.  Unlike  low  flow  angles,  however,  the  center  port  pressure 
is  now  the  most  dependent  port  on  local  flow  angle.  Consequently,  a 
pitch  angular  pressure  coefficient  should  measure  the  difference 
between  the  center  port  and  a  new  stagnation  port.  In  the  example 
below,  the  pressure  at  P^  approximates  the  local  total  pressure. 


-  P 


+  P' 


(17) 


The  expression  is  again  nond imens i ona 1 i zed  by  the  apparent  dynamic 
pressure.  The  average  of  and  P^  approximates  the  static  pressure 
and  is  relatively  independent  to  changes  in  roll. 

Although  the  average  is  independent  of  4,  the  difference  is 
sensitive  to  roll  angle.  As  the  probe's  azimuthal  orientation  changes, 
the  windward  pressure  rises  and  the  leeward  pressure  falls.  The  result 
is  a  roll  angular  pressure  coefficient  (example  for  port  4  is  Equation 
8)  that  is  also  nond imens i ona 1 i zed  by  the  apparent  dynamic  pressure. 


Pi  +  Ps  (18) 

2 

The  high  flow-angle  CQ  and  coefficients  are  translated  from 

their  low  flow  angle  counterparts  using  the  same  rationale  found  in  the 

development  of  CG  and  C 4, .  The  low-angle  CQ  and  coefficients  are 

changed  to  account  for  the  different  ports  that  represent  total  and 

static  pressures  in  the  high-angle  regime.  The  equations  for  CQ,  , 

C  ,  and  C„  are  as  follows: 
o’  q 


USAFA-TR-86-2 


F.  Division  of  Angular  Space 

The  most  obvious  difference  between  low  and  high  pressure 
coefficients  is  that  six  high  coefficients  are  needed  for  the  high 
angles.  This  fact  leads  to  the  question  of  what  factors  determine  when 
a  certain  set  of  coefficients  should  be  employed.  Gallington  proposes 
the  "division  of  angular  space"  shown  in  Figure  15  (Ref.  24). 


SECTOR  2 
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This  method  separates  probe  measurements  into  seven  sectors  --  a 
central  low  angle  sector  and  six  circumferential  high  angle  sectors. 
Data  points  are  placed  in  a  given  sector  based  on  the  highest  port 
pressure  measured  on  the  probe. 

G.  Polynomial  Power  Series  Expansion 

Once  the  pressure  coefficients  are  calculated,  they  must  be 
converted  to  ur ,  ,  CQ,  or  for  low  angle  flows  and  B  ,  ♦  ,  CQ,  or  Cq 

for  high  angle  flows.  This  conversion  is  accomplished  by  solving  the 
following  fourth  order  power  series: 

Order  of  Tens 


+  Oth 

+  “’Si +  l,t 

“kJi  *  “■SiSi +  +  !nd  (21) 

^ca,  ♦  tfSiSi ♦  “Vei +  “"Si +  M 

^‘Cn  *  rf<Sl  +  Kf’CaiC8l  +  +  *''Ch  “h 


"A"  signifies  the  desired  output  quantity  with  the  subscript  denoting 
the  ith  such  quantity.  "K"  denotes  the  presently  unknown  calibration 
coefficients.  The  calibration  process  entails  finding  these 
calibration  coefficients  with  the  A' s  as  known  conditions  and  the  C' s 
as  measured  pressure  coefficients  at  these  known  conditions. 
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H.  Determination  of  Calibration  Coefficients 

Rewriting  in  matrix  notation  for  n  terms,  the  power  series 

becomes 
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or  simply 


[A]  =  fC]  fK  ] 

Solving  for  Qk]  by  the  least  squares  curve  fit  outlined  by  Netter  and 
Vasserman  in  Gallington  (Ref.  24),  the  following  is  obtained: 


»* 


l. 


[K]  -  [CTC]_1[C]T[A]  (24 

With  known  calibration  coefficients,  the  probe  can  be  used  to 
determine  output  quantities  (  ur,  Sr,  CQ,  C^,  or  0  , 4  ,  CQ,  C^) 
explicitly  for  an  unknown  flowfield.  It  is  important  to  remember  that 
to  this  point  the  flowfield  has  been  uniform  and  incompressible. 
Temperature  or  velocity  gradients  have  not  been  accounted  for. 

Finally,  while  the  local  total  and  dynamic  pressures  are  not  found 
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explicitly  from  the  power  series,  they  are  easily  determined  by 
substituting  CQ  and  into  Equation  16  (low  angles)  or  Equation  20 
(high  angles). 

I.  Extension  to  Compressible  Flows 

Gerner  and  Mauer  (Ref.  2)  expanded  the  technique  to  subsonic 
compressible  flows  with  the  introduction  of  a  nondimens ional  pressure 
coefficient  representative  of  compressibility  effects.  This 
coefficient,  C^,  had  to  become  insignificant  at  very  low  Mach  numbers; 
thus,  any  terms  bearing  in  the  power  series  expansion  would  have  to 
approach  zero  as  Mach  approached  zero.  This  would  leave  essentially  an 
incompressible  power  series  with  the  regular  two  angular  pressure 
coefficients  as  previously  described.  Because  pressure  probes  are 
unable  to  determine  Mach  number  in  the  hypersonic  range,  had  to 
approach  a  finite  limit  at  very  high  Mach  numbers.  Consequently,  large 
changes  of  Mach  numbers  in  this  region  would  have  negligible  effect  on 
the  compressibility  coefficient. 

Gerner  and  Mauer  (Ref.  2)  found  that  these  requirements  were 
satisfied  by  the  apparent  dynami c- to- total -pressu re  ratio  (low  angles). 


r  _  p?  ~  pi~6 
M  “  P7 


(25) 


Similarly,  for  high  angle  flow,  the  compressibility  coefficient  in  each 
sector  is: 
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One  of  the  primary  problems  of  extending  the  seven-hole  probe's 
capabilities  into  compressible  flow  lies  in  the  mathematics.  First,  by 
adding  a  third  coefficient  to  the  fourth  order  power  series,  the  number 
of  terms  and  hence  calibration  coefficients  jumps  from  15  to  35. 

Gerner  and  Mauer  say  that  approximately  80  data  points  in  two  variables 
(C  and  CD  )  for  each  of  the  seven  sectors  are  needed  for  a  complete 

-  p 

incompressible  calibration.  That  fact  results  in  560  data  points,  and 
the  addition  of  a  Mach  number  compressibility  coefficient  makes  the 
data  set  unwieldy  (Ref.  2).  These  two  factors  create  complex  and  time 
consuming  matrix  operations.  Additionally,  the  amount  of  run  time 
needed  to  operate  the  wind  tunnel  for  all  of  the  calibration  data 
points  is  prohibitive  as  well  as  costly. 

Gerner  and  Mauer  suggest  a  two-part  solution  to  make  the  addition 
of  feasible.  By  reducing  the  fourth  order  power  series  to  a  third 
order,  the  number  of  calibration  coefficients  is  reduced  to  20.  Gerner 
and  Mauer  also  decrease  the  number  of  data  points  required  for 
calibration  by  employing  a  6x6  Latin  Square  technique  (Ref.  2)  for 
purposes  of  obtaining  calibration  data.  The  Latin  Square  is  a 
numerical  method  that  ensures  a  homogeneous,  random  sample  of  a 
three-dimensional  parameter  space. 

Seven  hole  probe  work  incorporating  these  two  changes  results  in 
the  following  conclusions  for  compressible  flow.  The  third  order 
polynominal  expansion  accurately  represents  the  parameter  space  up  to 
but  not  beyond  70  degrees  of  pitch.  Next,  based  on  reasonably  close 
correlations  between  the  standard  deviations  of  compressible 
calibrations  and  incompressible  calibration,  Latin  Squares  produce  a 
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sample  which  accurately  represents  a  very  large  three-dimensional 
parameter  space. 

With  these  changes,  the  seven-hole  probe  may  now  be  used  in  an 
unknown,  uniform  compressible  flow.  Since  the  seven-hole  probe  has 
never  been  subjected  to  supersonic  applications,  no  shocks  will  lie 
ahead  of  the  probe,  and  the  isentropic  flow  relationship  is  employed  to 
find  Mach  number  (Ref.  2). 

2  2  P  -P  (1-y)/y 

Mi  =  — (  [1  -1)  (27) 

Y  P0 

Knowing  CQ,  ,  and  the  seven  port  pressures,  the  local  Mach  number  may 
be  determined  explicitly  from  the  dynamic  to  total  pressure  ratio. 

This  ratio  for  the  inner  sector  is  represented  in  Equation  28  below 
(Ref.  2). 


[Cq 


'P7 


(28) 


Similar  equations  for  the  outer  sectors  may  be  derived,  using  the 
appropriate  port  pressures  for  approximate  total  and  dynamic  pressures. 
IV .  Seven-Hole  Probe  in  Shear  Flow 

Under  the  present  calibration  scheme,  the  existence  of  a  velocity 
gradient  in  the  flow  will  cause  the  probe  to  measure  an  erroneous 
uniform  flow  at  fictitious  flow  angle.  For  example,  a  flow  having  the 
properties  M= . 3 ,  u=0°,  0=0°,  with  a  velocity  gradient  (Fig.  16)  might 
cause  the  probe  to  "think  that  it  sees  a  flow  with  a  certain 
angularity,  saye=10°  (see  Fig.  17). 


Fiflur#  17.  APPARENT  FLOW  CONDITION 


Although  the  probe  has  seven  pressure  ports,  each  with  a  different 
pressure  coefficient,  all  seven  values  compute  a  single  point  property 
in  the  flow.  The  use  of  these  point  properities  in  a  uniform  flow 
calibration  produces  an  apparent  flow  condition  that  does  not  truly 
represent  the  actual  flow.  However,  by  comparing  adjacent  points  in  a 
flowfield,  a  velocity  gradient  can  be  determined.  Then  using  the 
"backs tepping"  technique  developed  in  the  following  sections,  the 
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apparent  measured  values  can  be  corrected  to  the  actual  flow  values. 

A.  Slender  Body  Theory 

The  slender  body  theory  developed  by  Huffman  is  utilized  in 
this  section  to  model  flow  around  the  probe  (Ref.  20).  With  this 
model,  the  probe's  port  pressure  coefficients  can  be  estimated  as  a 
function  of  flow  angularity.  Using  Huffman's  method  of  viewing  the 
probe  as  a  small  perturbation  to  an  otherwise  uniform  stream,  the 
pressure  coefficients  are  determined  by  analytic  means.  Further,  the 
model  determines  the  variation  of  pressure  coefficients  with  flow 
direction.  Not  included  are  viscous  or  flow  separation  effects  as  the 
theory  is  based  on  potential  flow. 

The  body-geometry  function,  f,  is  a  parameter  used  to  define  the 
geometric  characteristics  of  a  slender  body  and  is  given  in  Equation  2. 
Huffman  incorporates  the  body  shape  function  into  calculations  for  the 
pressure  coefficient  as  given  in  Equation  4.  A  simplified  form  of 
Equation  4  was  used  in  this  analysis  as  given  below. 

C  =  cos^  a  cos^  8[-f  -  (H')^J  +  sin^  8(1-4  sin^  0] 

+  sin^  a  cos^  8(1-4  cos^  0]  +  cos  a  sin  2 8 f  —  2 R '  cos  G]  (29) 

2 

+  sin  2a  cos  8[-2R'  sin  0]  +  sin  2a  cos  8(4  sin  0  cos  0] 

where  R'  is  dR/dX  with  R  as  the  radius  of  the  probe,  and  (u)  is  the 
angle  measured  from  the  y-axis  to  the  pressure  port  in  question.  The 
subscript  "c"  signifies  that  the  approximation  is  valid  for  uniform 
incompressible  flow. 
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The  function  that  Huffman  develops  for  the  pressure  coefficient  at 
port  7  simplifies  to  Equation  30. 

C  =  cosset  cos^B  -  sin^B  -  sin^a  cos^B  (30) 

p7 

c 

One  must  realize  that  slender  body  theory  is  only  an  approximate 
model  of  flow  around  the  seven-hole  probe.  The  reason  for  its  use 
stems  from  the  conditions  that  the  body  radius  is  much  less  than  the 
body  length  (slender  body),  and  that  it  is  easily  applied  to 
compressible  flows.  The  only  problem  with  this  approach  is  that  it 
assumes  that  the  rate  of  change  of  body  radius  with  respect  ot  body 
length  must  be  small.  The  seven-hole  probe  does  not  completely  conform 
to  this  assumption  for  two  reasons.  First,  the  probe  has  a  blunt  tip 
with  a  slope  discontinuity,  while  slender  body  theory  assumes  an 
aerodynamical ly  smooth  body  coming  to  a  point.  Because  the  tip  is 
blunt,  flow  disturbances  originate  at  that  tip.  Second,  due  to  the 
orientation  of  the  peripheral  hole  surfaces,  separation  occurs  on  the 
back  holes  at  moderate  angles  of  attack.  Slender  body  theory  assumes 
non - separ ated  flow  over  the  body  for  all  flow  angles.  Even  with  these 
discrepancies,  slender  body  theory  is  effective  in  predicting  basic 
trends  for  coefficient  of  pressure  data  as  shown  in  Figure  18.  As  a 
quick,  computational  technique,  it  generates  reasonable  pressure 
coefficient  approximations  for  the  following  shear  flow  analysis.  It 
can  also  be  used  to  show  that  in  a  uniform  flow  with  no  flow  angularity 
the  probe  tip  disturbances  (to  within  1%  of  the  freestream)  are 
confined  to  within  five  body  radii  regardless  of  probe  fineness  ratio 
or  Mach  number  (Mi  0.6). 


PREDICTED  AND  MEASURED 
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B.  The  Shear  Gradient 

In  order  to  work  with  a  two-dimensional  shear  gradient,  an 

alpha-beta  plane  is  defined  normal  to  the  direction  of  the  actual 

velocity  and  through  the  center  point  of  port  number  7  (Fig.  19a).  The 

velocity  at  port  7  is  used  as  a  reference  velocity;  hence,  =V^  for 

c 

all  gradient  analysis. 


Figure  19a.  SEVEN-HOLE  PROBE  IN  A  SHEAR  GRADIENT 


The  shear  gradient  is  defined  on  this  normal  plane.  The  gradient 
distance  is  measured  from  port  7  to  the  intersection  of  another 
velocity  vector  parallel  to  the  one  through  port  7  and  the  normal  plane 
(point  x,y,z). 
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By  determining  this  gradient  distance,  the  following  general  equation 
can  be  used  to  define  the  gradient’s  effect  on  the  velocity. 


.  dV 

+  y 

dy  J 


dV 

dz 


(31) 


where  y  and  z  are  the  components  of  the  gradient  distance  and  dV/dz  and 
dV/dy  are  indicators  of  the  magnitude  of  the  alpha  and  beta  shear 
gradients  respectively.  To  put  this  equation  in  a  more  usable  form, 
the  gradients  are  non-d imens i ona 1 i zed : 


V 


V7  d+ 


s  dV  2. 
V7  dy  s 


+ 


£  dV  z 
V;  dz  s 


(32) 


where  s  is  the  distance  between  the  centers  of  opposite  ports.  With 
this  new  equation,  the  terms  ~  11  and  11  can  be  assigned  a  magnitude 


of  relative  gradient. 
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C.  Gradient  Effects  on  Port  Pressure  Coefficients 

For  a  uniform  calibration  the  pressure  at  the  peripheral  holes 
may  be  found  by  the  following  expression. 


pi 


(1+1  IV 
V7  dy 


(1)  +  = 

s'  V7  dz  V' 


(33) 


Since  C  is  a  uniform  coefficient  based  on  the  geometry  of  the 

Pic 

slender  body  theory,  V.  is  the  only  term  that  accounts  for  gradient 

effects.  Consequently,  in  shear  flow  differs  for  each  peripheral 
hole.  If  the  velocity  term  is  kept  constant  at  and  the  gradient 

effects  are  accounted  for  by  the  coefficients  of  pressure,  the 
following  expression  results: 

Pi  -  P  +  C  (*ap V?)  (34) 

1  “  PA  1 

c 


Solving  Equation  33  and  34  for  C  ,  it 

^i 

proportional  to  the  square  of  velocity 


is  evident  that  pressure  is 
in  incompressible  flow. 


Pi  =  p«o  +  cPi  (*spv5) 


(35) 


Since  the  velocity  varies  proportionally  to  the  gradient  correction 
factor  (li-gradient  effects)  in  Equation  32,  the  pressure  coefficient 
varies  in  a  shear  gradient  by  the  following  relation: 

C  -c  (^)2  (36) 

Pi  Pl  V7 


D.  Determination  of  Point  (x,  y,  z) 

To  solve  the  previous  equation,  y  and  z  must  be  found  through 
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the  determination  of  point  (x,y,z).  The  vector  analysis  behind  the 
location  of  point  (x,y,z)  is  explained  in  detail  in  Johnson  and  Reed 
(Ref.  25).  The  results  are  as  follows: 

x  *  +  (xj  -  x cosa  cosB  -  cosa  cosB  [ y i  sinB  +  z*  sina  cosB] 

y  *  yi  +  (x  ~  xp  tanB/cosa  ^ 

z  «  +  (x  -  x^)  tana 


Port 


(i) 

xi 

yt 

2i 

1 

s/2tan  c 

(s/2)cos0j 

(s/2)sin0j 

2 

s/2tan  c 

(s/2)cos02 

(s/2)sin02 

3 

s/2tan  c 

(s/2)cos03 

(s/2)sln03 

4 

s/2tan  e 

(s/2)cos0^ 

(s/2)sin0^ 

5 

s/2tan  c 

(s/2)cos03 

(s/2)sin03 

6 

s/2tan  e 

(s/2)cos0g 

(s/2)sin0^ 

7 

d/2tan  c 

0 

0 

Since  x.,  y.,  z.,  and  Xj  are  known  values,  the  knowledge  of  alpha  and 
beta  will  locate  the  point  (x,y,z). 

E.  Shear  Flow  Measurement  Corrections 

At  this  point,  the  mathematical  relationships  necessary  to 
compare  the  effects  of  velocity  gradients  are  available.  Equation  15 
provides  the  ability  to  calculate  Cu,  C  .  and  Cb,  and  Cgc  realizing 
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By  comparing  the  actual  to  the  apparent  angular  pressure  coefficients, 
the  actual  angles  can  be  incremented  until  the  apparent  and  actual 
coefficients  are  equal.  This  procedure  indicates  the  flow  angle  errors 
inherent  in  the  measurement  of  velocity  gradients. 

In  order  to  determine  the  effects  of  different  shear  gradients  at 
varying  alphas  and  betas,  Johnson  and  Reed  (Ref.  25)  calculated  the 
apparent  flow  angles  resulting  from  a  span  of  actual  flow  angles  and 
shear  gradients.  The  calcuations  show  the  errors  between  the  apparent 
angles  and  actual  input  values.  Consequently,  the  actual  flow  angles 
can  be  determined  by  backs tepp i ng  from  the  apparent  angles.  Although 
these  calculations  only  allow  manual  corrections  to  apparent 
measurements,  the  theory  of  this  section  can  be  used  to  generate  a 
family  of  data  points.  These  points  may  then  be  input  into  a  computer 
surface  fitting  scheme  for  highly  accurate,  near  real-time  correction 
to  the  apparent  data. 

V.  Seven-Hole  Probe  Measurement  Capabilities 

Three  key  elements  of  the  seven-hole  probe's  measurement 
effectiveness  in  unknown  flows  are  discussed  further  in  order  to 
provide  an  understanding  of  the  basics  behind  its  actual  and  potential 
capabilities.  These  three  elements  to  be  discussed  are  the  probe's 
calibration,  construction,  and  measurement  error  sources. 

A.  Calibration 

Every  pressure  probe  must  be  uniquely  calibrated  due  to 
manufacturing  imperfections.  The  specifics  of  the  seven-hole  probe's 
calibration  scheme  has  been  outlined  in  great  detail  in  the  previous 
sections.  As  in  similar  schemes  for  other  multi-hole  probes, 
seven-hole  probe  calibration  allows  for  the  relative  flow  angles,  total 
pressure,  and  static  pressure  to  be  determined  explicitly.  With 
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computer  processing,  these  calculations  are  performed  in  near  real 
time. 

As  discussed  earlier,  problems  with  the  present  calibration  arise 
when  attempting  to  make  measurements  in  shear  flow.  This  is  because 
the  present  calibration  and  associated  calibration  coefficients  are 
determined  in  a  uniform  flow.  Jonas  recorded  measurement  discrepancies 
in  his  attempts  to  use  the  seven-hole  probe  to  map  unknown  flowfields 
with  suspected  shear  (Ref.  26).  Specifically,  he  compared  seven-hole 
and  total  pressure  probe  measurements  in  a  vortex  wake  created  by  the 
wing  leading-edge  extensions  of  a  Northrop  VATOL  model. 
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In  a  vortical  flow,  local  total  pressure  should  decrease  through 

the  wake  because  of  viscous  effects.  This  requires  the  coefficient 
P  P 

^TOTAL  (  °L  ~ _ -)  t0  a  negative  value.  The  results  of  Jonas '  pressure 

O  ® 

measurement  comparison  are  shown  in  Table  1  below  (Ref.  26). 

Total  Pressure  Probe 


Seven-Hole 
Pressure  Probe 

Ai  i gned  with 
Freestream 

Al  igned  with 
Local  Flow 

Number  of 

Readings  (50 

Data  Samples 

Per  Reading) 

1 

10 

69 

Average  Value 

+0.197 

-0.046 

-0.013 

Dev iat i on  ,o 

+0.003 

+0.0237 

+0.023 

Maximum  Value 

-- 

-0.003 

+0.040 

JRtotal  . 

Minimum  Value 

-- 

-0.067 

-0.055 

Pe9IeKtage  of 

Positive  CTQTAL 
Read i ngs 

-- 

0% 

26% 

Table  1 


Jonas  completed  the  following  discussion  based  on  the  results  of 
Table  1  (Ref. 26).  As  can  be  seen  from  the  data,  the  total  pressure 
probe  measurements  do  not  agree  with  seven-hole  probe  pressure 
measurements.  For  the  total  pressure  probe  aligned  with  the  freestream 
flow,  no  measurements  existed  where  was  positive.  However, 

these  values  may  not  be  an  accurate  measure  of  local  total  pressure 
because  flow  angularity  may  have  caused  separation  at  the  probe  tip. 

For  the  total  pressure  probe  aligned  with  the  local  flow,  26  percent  of 
the  measurements  resulted  in  a  local  total  pressure  greater  than  that 
of  the  freestream  value.  These  measurements,  however,  are  lower  than 
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those  obtained  using  the  seven-hole  probe.  One  must  realize  that  the 
total  pressure  probe  was  bent  at  an  angle  based  on  measurements  made 
using  the  seven-hole  pressure  probe.  Since  the  seven-hole  measurements 
were  made  in  a  region  of  high  shear  flow,  the  local  flow  angularity 
reading  may  also  be  erroneous.  Consequently,  the  flow  around  the  tip 
of  the  total  pressure  probe  supposedly  aligned  with  the  local  flow 
might  also  be  separated,  again  leading  to  false  pressure  readings.  The 
average  value  for  Cygy^  of  the  bent  total  pressure  probe,  however, 
remains  negative.  This  evidence  seems  to  confirm  our  theoretical 
expectations  of  viscous  effects. 

Although  the  contr ad ict i ve  probe  readings  seem  to  indicate 
incorrect  probe  measurements  in  the  high  shear  regime,  the  possibility 
of  positive  CyOTAL  regions  may  not  be  totally  excluded.  A  vortex  flow 
is  very  complex  and  not  completely  understood.  Comparisons  in  a  low 
shear  flow  environment  (airfoil  wake)  of  seven-hole  probe  measurements 
and  hot-wire  anemometer  measurements  are,  however,  very  favorable  and 
no  discrepancies  or  positive  Cy^y^  measurements  exist  (Ref.  26).  This 
seems  to  indicate  that  the  problem  of  positive  Cy^y^’s  are  only 
associated  with  high  shear  flows,  and  are  either  due  to  the  probe 
(manufacture,  calibration,  etc.)  or  actually  represent  a  local  flow 
phenomena  not  clearly  understood.  Whatever  the  case,  care  must  be 
taken  when  making  measurements  in  high  shear  flow  environments. 

B.  Construction 


A  pressure  probe  should  be  as  small  as  possible  to  keep  the 
flow  it  is  measuring  undisturbed;  however,  this  is  accompanied  by 
several  construction  disadvantages.  For  example,  smaller  probes 
require  delicate,  precise  machining  techniques.  As  size  decreases, 
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burrs  and  other  manufacturing  imperfections  cause  greater  probe  tip 
construction  errors,  and  the  port  tubes  are  more  likely  to  leak. 

Although  the  five-hole  probe  may  seem  easier  to  construct  than  the 
seven-hole,  the  exact  opposite  is  true.  The  construction  of  the 
seven-hole  probe  is  greatly  simplified  since  the  port  tubes  are 
arranged  in  the  only  geometric  fitting  possible.  Gallington  and 
Hollenbaugh  (Ref.  27)  explain  three  difficulties  that  this  geometric 
arrangement  overcomes,  difficulties  which  still  exist  in  the  five-hole 
probe. 


Figure  21.  SEVEN-HOLE  PROBE  CONSTRUCTION 
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SOLDER 


Figure  22.  CONVENIENT  PROBE  DESIGN 


First,  it  is  difficult  to  hold  five-hole  tubing  parallel  with  the 
center  tube  and  in  perfect  azimuthal  position  for  soldering.  Second,  a 
high-powered  magnifier  is  required  to  assure  equal  chamfer  on  the  four 
side  tubes.  Finally,  there  is  no  guarantee  that  the  four  solder 
fillets  shown  below  will  be  equal  or  even  nearly  equal. 


Figure  23.  FIVE-HOLE  PROBE  ASSEMBLED 
FROM  TUBES 
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As  a  quick  note,  a  simpler  geometric  arrangement  for  five  tubes  is 
shown  below.  This  order,  however,  lacks  the  required  center  tube. 


Figure  24.  FIVE  TUBES  IN 
THE  SMALLEST 
CIRCUMSCRIBED 
CIRCLE 


Over  the  years,  multi-hole  probe  construction  has  uncovered 
certain  features  that  increase  the  accuracy  of  pressure  measurement. 
Total  pressure  measurement  can  be  largely  desensitized  to  alpha  and 
beta  by  flaring  the  stagnation  port  at  an  angle  of  30  to  60  degrees. 
Yet,  the  probe's  overall  sensitivity  to  flow  angularity  can  be 
increased  by  decreasing  the  tip  half-angle  (e).  The  drawback  to  a 
small  half-angle  is  that  the  flow  will  separate  from  the  leeward 
surface  when  ur  exceeds  e  (Ref.  1). 

C.  Measurement  Error  Sources 

Once  smaller  probes  are  built,  they  are  more  sensitive  to 
measurement  errors  caused  by  flow  debris  and  damage.  These  two  effects 
can  be  minimized  prior  to  use  by  reverse  airflow  through  clogged  ports, 
by  dust  covers  over  port  entrances,  and  by  simple,  careful  handling  of 
sma 1 1 er  pr obes . 


48 


USAFA-TR-85-2 


Huffman  identifies  three  other  error  sources  that  adversely  affect 
a  probe's  measurement  of  pressure,  orientation,  and  velocity  (Ref.  1). 
These  sources  include  the  time  lag  between  sensed  and  actual  pressure, 
the  pressure  transducer  resolution  and  frequency  response,  and  the 
resolution  of  the  ana  1 og- to-d i g i ta 1  conversion.  The  pressure 
transducer  resolution  and  ana  1 og - to - d i g i t a  1  resolution  are  both  beyond 
the  scope  of  this  report;  however,  the  other  error  sources  are 
discussed  in  the  following  paragraphs  (Ref.  1). 

As  flowfield  conditions  change,  the  probe  surface  pressures  change 
nearly  instantaneously.  These  pressure  changes  are  transmitted  to  the 
pressure  transducers  by  a  finite  amount  of  air  through  a  connection 
tube.  Huffman  discusses  the  fact  that  the  pressure  lag  caused  by  this 
finite  travel  time  is  related  to  the  speed  of  pressure  propagation  and 
the  pressure  drop  due  to  the  viscous  effects  of  the  tubing  (Ref.  1). 
Three  physical  characteristics  directly  affect  this  time  lag: 

1 .  Tubing  d i ame t er 

2.  Tubing  length 

3.  Transducer  cavity 

Decreasing  the  tubing  diameter  and  increasing  the  tubing  length  both 
increase  the  viscous  effects  on  the  air  and  increase  the  pressure  time 
lag.  For  a  larger  transducer  cavity,  more  time  is  required  for  the 
pressure  pulse  to  propagate  from  the  tubing  exit  to  the  transducer 
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Schlapkohl  and  Buzzell  (Ref.  28)  define  frequency  response  as 

the  inverse  of  the  maximum  dwell  time  necessary  for  the 
probe  system  to  react  to  the  maximum  pressure  difference 
expected  where  the  pressure  sensed  by  the  pressure 
transducer  reaches  99  percent  of  the  actual  (surface) 
pressure. 

The  frequency  response  of  the  transducer,  the  computational  time 
required  to  convert  the  measured  pressure  to  the  desired  output,  the 
pressure  differential  between  the  probe  tip  and  the  transducer  face, 
and  the  fluid  density  all  affect  the  probe's  overall  frequency  response 
(Ref.  28). 

Probe  frequency  response  rates  are  generally  of  a  few  Hertz. 
Because  modern  transducers  -and  microprocessors  require  only 
milliseconds  to  operate,  their  operation  is  essentially  instantaneous, 
and  little  can  be  accomplished  in  these  areas  to  decrease  the  overall 
frequency  response.  The  final  two  factors  are  associated  with  the 
fluid  dynamic  properties  of  the  flow.  The  difference  between  the  fluid 
pressure  at  a  port  entrance  and  the  pressure  at  the  transducer  face 
acts  as  a  driving  potential,  and  a  greater  differential  increases  the 
frequency  response  (Ref.  28).  Since  a  higher  density  increases  the 
collision  rate  between  molecules,  and  therefore,  increases  the 
propagation  rate  of  a  pressure  pulse,  greater  density  is  synonymous 
with  greater  frequency  response.  In  considering  the  effect  of  density, 
one  must  remember  the  assumption  of  a  constant  static  pressure  across 
the  face  of  the  probe.  Remembering  P=  PRT,  one  should  realize  that  the 
density  of  the  flow  is  dependent  upon  pressure. 
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lications:  Measurements  of  Unknown  Flows 


Research  efforts  with  the  seven-hole  probe  center  on  two 
objectives:  1)  to  establish  the  capabilities  and  limitations  of  the 
seven-hole  probe  and  2)  to  map  unknown  flowfields.  The  following 
section  explains  how  these  two  objectives  are  accomplished  in  the 
analysis  of  a  simple  vortex  flow. 

A.  Vortex  Flows 

To  create  the  vortex  flow,  a  finite  wing  section  of  a  NACA 
0008  airfoil  was  placed  in  the  two  by  three  foot  subsonic  wind  tunnel 
at  the  United  States  Air  Force  Academy  (Figure  25).  The  tests  were  run 
with  the  wing  at  an  angle  of  attack  of  8  degrees  and  a  uniform  flow 
velocity  of  150  ft/sec. 


iso  n/i»c 

uniform 

flow 


Figure  26.  WING  VORTEX  TEST 
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Computer  graphing  techniques  were  used  to  map  out  the  flow  velocities 
and  direction  for  all  the  points  in  a  series  of  parallel 
two-dimensional  planes,  all  of  which  were  perpendicular  to  the  uniform 
f  1  ow. 


Figure  26.  WINQ  VORTEX  OATA  PLANES 
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The  two-dimensional  planes  or  test  grids,  started  twelve  inches  aft  of 
the  wing's  trailing  edge  and  moved  forward  (toward  the  wing)  by  two 
inch  increments  for  every  data  set.  From  these  test  runs,  capabilities 
and  limitations  of  the  seven-hole  probe  under  mild  gradient  conditions 
are  examined  and  a  wingtip  vortex  is  mapped. 

B.  Crossflow  Velocities 

The  first  step  in  vortex  analysis  is  to  examine  the  crossflow 
velocities  in  each  plane.  Cross  flow  velocity  plots,  which  are  scaled 
by  a  factor  of  two  are  shown  in  Figures  26a  to  27h  (Ref.  29).  With 
past  errors  in  gradient  regime,  one  would  expect  errors  in  seven-hole 
readings.  Suspect  data  of  this  nature  (Ref.  29)  are  identified  by 
inconsistencies  in  the  crossflow  plots,  such  as  an  inordinately  large 
or  small  arrow  (Fig.  28)  or  an  arrow  in  the  wrong  direction.  However, 
the  actual  crossflow  plots  do  not  confirm  these  expectations.  The 
velocities  are  well-behaved,  and  no  incons i s tancies  are  seen  near  the 
vortex  core  where  high  shear  exists. 
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Figure  28.  CROSSFLOW  DATA  ERRORS 


The  complete  set  of  velocity  data  may  be  shown  on  a  velocity 
magnitude  centerline  plot  for  ease  of  analysis  (Figure  29). 
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C.  Definition  of  Flow  Field  Pressure  Coefficients 

Because  pressure  measurements  are  affected  by  slight  shifts  in 
wind  tunnel  velocity  and  temperature,  pressure  coefficients  are  used 
when  the  local  and  total  tunnel  pressures  are  compared  and 
non-d imens i ona 1 i zed  by  the  tunnel  dynamic  pressure.  These  three 
coefficients  are  defined  below: 


^STATIC 

CTOTAL 


CDYN  =  CTOTAL  ~  CSTATIC 


(39) 


o  00 

Probe  measurement  errors  may  be  found  by  determining  the  correct 
values  for  these  three  pressure  coefficients,  either  theoretically  or 
experimentally  with  another  device,  and  comparing  these  values  with 
actual  seven-hole  probe  data.  Vortex  pressure  coefficient  trends  may 
be  found  through  analysis  of  a  two-dimensional  vortex. 


D’  CTOTAL 

For  an  ideal  two-dimensional  vortex,  angular  velocity 
increases  inversely  as  the  distance  to  the  vortex  center  decreases 
(Figure  31).  Real  vortices,  however,  are  subject  to  viscous  effects. 
Outside  of  point  A,  the  flow  is  essentially  inviscid.  Inside  that 
point,  viscous  forces  reduce  angular  velocity  until  it  reaches  zero  at 


the  vor  tex  center . 
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INVISCID  VISCID 

FLUID  FLUID 


Flflure  31.  ANGULAR  VELOCITY  OF  VORTEX 


These  two  flow  regions  directly  affect  values.  Because 

total  pressure  is  constant  in  incompressible,  inviscid  flow, 
zero  in  the  inviscid  region.  On  the  other  hand,  the  viscous  forces 
found  within  the  vortex  core  decrease  the  flow's  fluid-mechanical 
energy.  PQ  decreases  causing  to  become  negative.  Figure  32 

diagrams  CTnTA|  's  expected  behavior. 
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Figure  32.  VORTEX  CT0TAL 


The  expected  shape  was  derived  by  Jonas  (Ref.  26)  from  theoretical 
considerations  of  a  2-D  vortex  decaying  with  time  (Equation  40). 
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As  discussed,  the  bucket-shaped  phenomenon  is  present  in  the 
experimental  vortex  data  shown  in  Figure  33. 
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As  the  flow  progresses  downstream,  lower  velocity  gradients  produce 
smaller  viscous  forces,  smaller  deficits  of  f 1 u i d -mechan i ca 1  energy, 


and  small  bucke t - s haped  area  plots.  The  two  areas  of  positive 
in  each  of  the  centerline  plots  is  unexpected.  The  fact  that  the 
magnitude  of  the  positive  Cy^y^y  values  is  always  greater  on  the 
outboard  side  of  the  wing  than  on  the  inboard  side  is  inexplicable. 

The  sets  of  contour  plots  and  axonometric  plots  (Appendix  A)  show  that 

Ctqtal  *s  P05*1'^  only  on  the  two  sides  or  the  vortex  parallel  to  the 
wing.  The  vortex  data  behaved  normally  above  and  below  the  wing.  In 
duplicating  Jonas's  previous  positive  Cygy^y  results  the  vortex  data 


seem  to  indicate  a  limitation  in  the  seven-hole  probe's  measurement 


capabilities.  It  is  still  possible,  however,  that  the  positive  C 
values  may  be  the  result  of  a  transfer  mechanism  that  is  not  yet 
unde r  s  t ood . 


TOTAL 


E.  C 


STATIC 


TOTAL 


values  of  C 


Like  CJOTAL,  the  presence  or  absence  of  viscosity  affects  the 
values  of  C^yAy((-..  At  freestream  conditions,  CsyA7]C  equals  zero. 
Progressing  inward  from  the  freestream  conditions  to  point  A,  the 
angular  velocity  increases,  static  pressure  drops,  and  C^y^yj^-. 
decreases  from  zero  to  a  negative  value.  To  determine  what  happens  in 
the  viscous  region,  the  following  mathematical  analysis  is  required, 
dp  cu,! 

dF= 


For  a  free  vortex  u  =  —  so,  iL  _  il£. 

0  r  d  r  r 


For  solid  body  rotation  u  *  r  so ,  ^  =  mW  r 
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Integrating  for  p  results  in  the  following: 


P  = 


(41 


P 


c2  +' 


2  2 

pw  r* 


(41 


Equations  4la  and  41b,  when  superimposed,  reveal  a  bucket-shape  like 
that  of  The  entire  curve  is  shown  in  Figure  34. 


Figure  34.  VORTEX  C  STATIC 


F*  CDYNAMIC 

DYNAMIC  is  simP,y  the  subtraction  of  CSTATIC  from  CTQTAL  as 
shown  in  Equation  39.  Again,  at  the  freestream,  yn /VM I C  's  *c*ea**y 
zero.  ^DYNAMIC  equals  the  absolute  value  of  csyA7jC  from  the 
freestream  conditions  inward  to  point  A  since  C70TAL  is  zero.  At  the 
center  of  the  vortex,  local  velocity  and  therefore  local  dynamic 
pressure  are  zero  causing  ^DYNAMIC  t0  ec^ua*  From  point  A  to  the 

center,  CDyNAMIC  changes  as  in  Figure  35. 
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C  DYNAMIC 

Figure  35.  VORTEX  C  DYNAMIC 


Experimental  data  support  the  theoretical  analysis  of  CDYN .  This 
seems  somewhat  surprising  in  lieu  of  the  questionable  positive 
values.  The  important  factor  missing  from  the  coefficient  plots  are 
the  relative  magnitudes.  Because  values  are  very  close  to  zero 

the  contribution  of  CyoTAL  is  not  very  significant  in  the  calculation 

of  CDYNAMIC’ 
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G.  Research  Applications 

Today's  flow  mapping  efforts  are  much  more  complex  than  the 
vortex  example.  Besides  Jonas'  VATOL  flowfield  measurements, 
seven-hole  probe  work  at  the  Air  Force  Academy  has  examined  canard 
wakes,  lifting  surface  wakes  of  canard/swept  wing  aircraft,  and 
flowfield  characteristics  of  square  cross-sectional  missle  bodies. 

Probe  calibration  and  measurement  of  unknown  flowfields  have  also  been 
conducted  at  the  NASA  Ames  Research  Center  in  the  2x2  foot  and  14  foot 
transonic  wind  tunnels  (Ref.  30).  Measurements  of  wing  and  canard 
jet-flap  effects  as  well  as  the  effects  of  prop  fan  installations  have 
been  made  with  multiple  and  single  probe  installations.  Vet,  these 
more  complex  research  efforts  are  founded  on  the  basic  analysis  and 
data  presentation  explained  in  the  vortex  example. 

Griffin  (Ref.  31)  used  cross  velocity  and  pressure  contours  in  his 
study  of  canard/ forward  swept  wing  aircraft.  Figures  37  and  38  show 
the  two  types  of  data  plots  in  the  same  spatial  relationship  that  their 
data  points  have  to  the  model. 


Figure  37.  TYPICAL  VELOCITY  DATA  Figure  38.  TYPICAL  PRESSURE  DATA 

PLANE  LOCATION  PLANE  LOCATION 

RELATIVE  TO  MODEL  RELATIVE  TO  MODEL 


Jonas  (Ref.  26)  superimposed  these  two  types  of  plots  in  his 
investigation  of  (see  Figure  39). 
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Figure  39.  REGIONS  OF  POSITIVE  C  T0TAL  .  DATA  PLANE  ZMP75 


Jonas  also  used  a  series  of  axonometric  projections  to  analyze  the 
growth  (diffusion)  and  decay  (dissipation)  of  vortices.  Evidence  of 
mixing  is  found  in  the  relative  flattening  of  the  C jq-j- ^  values 
(negative  CT~T..  being  up  or  out  of  the  plane  of  projection). 
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Gerner  and  Durston  at  NASA  Ames  took  Jonas'  data  and  produced  a 
series  of  12  color  contour  photographs  that  span  the  development 
through  decay  of  the  VATOL  vortices.  The  study  was  limited  to  the 
examination  of  local  total  pressure  or  CfoTAL ’  The  color  contours 
allow  more  detail  in  data  representations  by  making  pressure 
differences  easier  to  see.  Regions  of  positive  CjoTAL  are  clearly 
distinguished  from  other  points  in  the  flowfield.  As  seen  in  the  other 
data  schemes,  the  color  contours  reveal  that  the  VATOL  is  in  a  slight 
sideslip.  This  causes  the  right  vortex  to  burst  prior  to  the  left 
vortex. 

VII.  Cone  1  us i ons 

As  shown,  the  seven-hole  pressure  probe  remains  a  valuable 
measurement  tool  for  the  documentation  of  unknown  flowfields.  The 
device  has  a  greater  measurement  range  and  flexibility  than  other 
similar  obtrusive  f 1 ow  measu r i ng  devices.  The  seven-hole  probes 
themselves  are  easily  constructed  and  calibrated  for  use  in  subsonic 
compressible  flows.  Measurements  in  adverse  shear  flows  can  be 
corrected  to  give  actual  flow  conditions  based  on  the  methodology 
developed  in  this  paper.  Finally,  as  shown  in  the  application  section, 
the  probe  is  not  only  a  valuable  research  device  but  educational  as 
well  in  demonstrating  fundamental  flow  properties. 
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A. 

i 


Symbol s 

the  ith  value  of  a  particular  data  point  either  known 
determined  from  calibration  equations 


CDYN 

CApM 

C 

o 

CP 

S 

Static 

CTOTAL 

Cu 

ce 

C, 

c0 

f 

Ki 

1 

M 

P 

q 

( r ,0 ,z) 

R 

s 

u 

V 

(x,y  ,z) 


H 


local  dynamic  pressure  coefficient 

coefficient  representative  of  compressibility  effects 
apparent  total  pressure  coefficient 
coefficient  of  pressure 
apparent  dynamic  pressure  coefficient 
local  static  pressure  coefficient 
local  total  pressure  coefficient 
angle  of  attack  pressure  coefficient 
angle  of  sideslip  pressure  coefficient 
roll  angle  pressure  coefficient 
pitch  angle  pressure  coefficient 
body  geometry  function 
calibration  coefficients 
length 
Mach  number 
Pressure 

dynamic  pressure 
cylindrical  coordinates 
radius 

distance  between  centers  of  opposing  holes 
local  velocity,  perturbation  velocity  component 
velocity 

Cartesian  coordinates 
angle  of  attack 
angle  of  sideslip 
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6 

e  probe  tip  half  angle 

K  ratio  of  specific  heats 

f  ro 1 1  ang  1  e 

p  dens i ty 

9  pitch  angle 

subscr i pts 

c  uniform  incompressible  flow 

(i,j)  local  property  or  condition 

o  total  or  stagnation  condition 

00  freestream  conditions 
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WIND  TUNNEL  STUDY  ON  THE  DRAG  REDUCING  EFFECTS  OF  AN  AFT-MOUNTED 

VENTRAL  FIN  ON  A  T-38  MODEL 
Gregory  W.  Bice* 

Editor's  Note 

This  article  is  the  result  of  the  author's  Aero  499  Independent  Study 
Project.  It  will  also  be  presented  at  the  1985  American  Institute  of 
Aeronautics  and  Astronautics  Region  V  Student  Paper  Competition. 

Abstract 

The  effects  of  an  aft-mounted  ventral  fin  on  a  1/18.4  scale  T-38 
model's  subsonic  drag  characteristics  were  investigated  in  the  2'x3' 
Subsonic  Wind  Tunnel  at  the  USAF  Academy.  Wind  tunnel  tests  were 
performed  on  the  baseline  model  using  boundary  layer  trips  and  with  a 
ventral  fin  and  clipped  vertical  tail  as  configuration  variables.  Data 
were  taken  to  determine  the  effects  of  the  configuration  changes  on  the 
aerodynamic  coefficients  and  stability  derivatives.  Data  were  obtained 
by  varying  angle  of  attack  and  sideslip  angle  of  the  model  in  separate 
test  runs  at  100  fps.  A  six  component  internal  balance  was  used  to 
measure  force  and  moment  data.  Test  results  showed  that  a  consistent 
reduction  in  drag  coefficient  of  up  to  5%  resulted  from  the  addition  of 
the  ventral  fin  to  the  baseline  model  at  low  angles  of  attack.  The 
drag  coefficient  increased  up  to  17%  over  the  baseline  model  with  the 
ventral  fin  when  the  vertical  tail  was  clipped.  Results  of  oil  flow 
tests  showed  no  significant  vortex  or  turbulent  wake  patterns  for  any 
con f i gu r at i on  tested.  The  exact  nature  of  the  cause-effect 
relationship  by  which  ventral  fin  reduced  the  subsonic  drag  coefficient 
remains  unknown. 


I .  Objective 

The  purpose  of  this  study  was  to  measure  the  reduction  in  drag  due 
to  the  addition  of  an  af t -moun t ed  ventral  fin  on  a  1/18.4  scale  T-38 
model.  In  addition,  1  a t er al / d i r ec t i ona I  static  stability  derivatives 
were  compared  for  model  configurations  with  a  normal  vertical  tail  and 
with  a  clipped  vertical  tail.  Finally,  the  cause-effect  relationship 
of  the  ventral  fin's  drag  reducing  properties  was  analyzed  and 
documen  ted . 
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I  I .  Theor  y 

Aerodynamic  drag  can  be  divided  into  three  major  categories: 
parasite,  induced,  and  wave  drag.  Since  our  wind  tunnel  tests  were 
conducted  at  subsonic  speeds,  wave  drag  is  not  present.  Induced  drag, 
or  drag-due-to  lift,  is  primarily  dependent  on  the  wing  aspect  ratio 
and  wing-body  combination.  Since  all  test  runs  were  conducted  at 
identical  alpha  sweeps  with  no  changes  to  the  wing,  the  induced  drag 
can  be  considered  a  constant  for  all  alpha  testing.  This  leaves  only 
parasite  drag  to  consider. 

At  angles  of  attack  below  stall,  drag  data  for  many 
configurations,  such  as  the  T-38,  can  be  expected  to  fit  the  classic 
drag  curve  equation, 


CD  =  CD 


KC, 


(1) 


where  KC^  is  the  induced  drag  coefficient  and  CQ  is  the  parasite  drag 

o 

coefficient  (Ref.  1).  Parasite  drag  is  further  divided  into 

interference  drag  and  profile  drag.  Both  are  functions  of  the 

airplane's  shape,  surface  area,  and  configuration.  Profile  drag  is 

composed  of  drag  due  to  skin  friction  and  flow  separation.  Since  the 

skin  friction  coefficient,  Cn  ,  is  a  function  of  Reynolds  number,  it 

uf 

can  be  assumed  that  Cn  remains  essentially  constant  because  all  tests 

u{ 

were  conducted  at  almost  the  same  velocity  and  temperature  (Ref.  I). 

One  purpose  of  modifying  this  aircraft  by  adding  the  ventral  fin  is  to 
try  and  reduce  the  drag  due  to  flow  separation  in  the  tail  area  of  the 
T-38.  Looking  at  a  figure  of  the  T-38,  a  slight  upward  slope  or 
boattail  can  be  seen  at  the  tail.  The  ventral  fin  is  to  be  added  in 
this  location.  By  adding  this  ventral  fin,  it  is  hoped  that  the 
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airflow  will  somehow  be  changed,  causing  a  decrease  in  separation  or 
form  drag . 

As  the  flow  moves  aft  on  the  surface  of  the  plane,  it  will  lose 
kinetic  energy  in  the  boundary  layer  due  to  skin  friction.  This 
decrease  in  kinetic  energy  in  the  form  of  reduced  velocity,  combined 
with  the  adverse  pressure  gradient,  causes  the  flow  to  work  "uphill"  in 
the  tail  area.  Consequently,  the  boundary  layer  may  separate  somewhere 
along  the  plane,  causing  a  separated  flow  region  and,  hence,  separation 
drag. 

Besides  using  the  ventral  fin  as  a  variable,  boundary  layer  trips 
were  also  used.  Many  models  have  such  smooth  surfaces  that  the 
boundary  layer  flow  remains  laminar  and  never  transitions  to  turbulent. 
The  main  reason  it  remains  laminar,  however,  is  due  to  the  Reynold's 
numbers  obtained  in  the  wind  tunnel  which  are  typically  one  to  two 
orders  of  magnitude  lower  than  those  of  the  actual  plane  in  flight. 

The  purpose  of  these  trips,  which  are  seen  in  Figure  3  as  the  thin 
strips  of  white  tape,  is  to  transition  the  boundary  layer  flow  from 
laminar  to  turbulent  at  the  same  position  that  it  would  occur  on  the 
actual  plane,  simulating  more  of  a  "real"  situation.  Tripping  the 
boundary  layer  to  turbulent  flow,  more  typical  of  actual  flight,  should 
result  in  higher  skin  friction  drag  but  lower  separation  drag. 

In  adding  the  ventral  fin  to  the  T-38,  elements  other  than  drag 
had  to  be  addressed,  such  as  the  static  stability  of  the  aircraft. 
Stability  considerations  are  very  important  because  they  affect  the 
flight  control  system  and  the  pilot's  ability  to  fly  the  plane.  In  the 
area  of  longitudinal  stability,  the  fin  should  not  affect  the  plane,  at 
least  theoretically,  because  it  does  not  contribute  to  any  longitudinal 
pitching  moments  about  the  center  of  gravity.  However,  in  the  area  of 
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1  a t er a  1 -d i r ec t i ona 1  stability,  the  fin  will  have  an  effect.  Since  it 
is  mounted  aft  of  the  plane's  center  of  gravity,  the  fin  should  have  a 
stabilizing  effect  (Cn^>0)  on  weathercock  stability,  meaning  a  positive 
sideslip  will  cause  a  positive  yawing  moment.  The  fin  should  also 
cause  a  decrease  (C,  <0)  in  the  lateral  or  roll  stability  of  the  plane. 

1  e 

This  implies  a  positive  sideslip  will  cause  a  positive  rolling  moment 
which  is  destabilizing  according  to  the  conventions  of  aerodynamics 
(Ref.  2). 

III.  Apparatus  and  Procedu  re 

This  study  was  conducted  in  the  2'x3'  Subsonic  Wind  Tunnel  at  the 
U.S.  Air  Force  Academy.  The  tunnel,  shown  in  Figure  1,  consists  of  a 
200  HP  synchronous  motor  with  a  Gyrol  fluid  drive  unit,  a  four-bladed 


compressor,  flow  conditioning  regions  such  as  the  stilling  chamber,  a 
converging  nozzle,  a  test  section,  and  finally  a  diffuser.  Velocities 
between  50  fps  and  400  fps  can  be  reached  in  the  clean  test  section 


along  with  dynamic  pressure  ranging  from  1.8  psf  to  130  psf.  During 
this  study,  a  velocity  of  100  fps  was  used  which  resulted  in  a  test 
Reynolds  number  of  about  0.2  x  lO*’.  This  airspeed  does  not  appear  to 
be  representative  of  T-38  flight  speeds.  However,  at  subsonic  Mach 
numbers,  which  is  the  flight  regime  most  used  by  the  T-38,  aerodynamic 
force  coefficients  stay  relatively  constant,  making  airspeed 
irrelevant.  The  velocity  was  limited  due  to  model  support  system 
dynamics  because  serious  vibrations  resulted  at  higher  airspeeds. 

The  T-38  model  was  supported  by  a  sting  and  model  support  system 
shown  in  Figure  2.  Variations  in  the  angle  of  attack  were  achieved 
using  a  pitch  sector,  strut  and  boom  assembly,  and  an  electric  motor. 
Sideslip  angle  tests  were  accomplished  by  rolling  the  model  and  balance 
90°  and  by  varying  the  sideslip  angle  while  holding  a  constant  zero 
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degree  angle  of  attack. 


TURNING  VANES  (TYPICAL) 


CHAMBER 


Figure  1.  General  Schematic  of  fh  the  2'  x  3'  Subsonic  Wind  Tunnel 
at  USAFA. 
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The  model  used  in  the  test  section  was  a  1/18.4  scale  T-38  with  a 
NACA  65A004.8  (mod)  airfoil  for  the  wing  sections.  Some  of  the  model 
dimensions  are  as  follows: 

1  eng  th  =  28.76  in 
wing  span  =  16.45  in 
MAC  =  5.04  in 

2 

Wing  planform  area  =  72.17  in 
CG  position  =  19.4%  MAC 

The  configurations  tested  included  combinations  of  boundary  layer 
trips,  the  ventral  fin,  and  the  clipped  vertical  tail.  Two 
configurations  are  shown  in  "igure  3.  In  order  to  get  both 
longitudinal  and  lateral  directional  data,  angle  of  attack  (alpha) 
sweeps  and  sideslip  angle  (beta)  sweeps  were  performed  using  the  model 
suppor  t  sys  tern. 


Figure  3.  T-38  Configurations  (a)  normal, 

(b)  with  Ventral  Fin  and  Clipped  Tail. 


The  purpose  of  testing  the  airplane  with  a  clipped  vertical  tail 
was  to  try  to  return  Cn  to  its  original  value.  In  other  words,  the 
tail  was  clipped  to  compensate  for  the  increased  weathercock  stability 
of  the  airplane  resulting  from  the  addition  of  the  ventral  fin. 
Clipping  the  vertical  tail  should  also  cause  a  decrease  in  roll 
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stability,  Cj  ,  which  will  probably  be  compounded  by  the  ventral  fin 
because  these  two  surfaces  are  180  degrees  apart  on  the  airplane. 
Because  of  the  reduced  surface  area  of  the  tail,  the  overall  drag  on 
the  T-38  should  decrease,  unless  the  tail  serves  some  purpose  in 
reducing  separation  drag. 

Force  and  moment  data  were  acquired  using  a  0.75"  diameter,  six 
component  steady-state  internal  strain-gage  balance,  seen  in  Figure  4. 
The  force  balance  is  mounted  on  the  end  of  the  sting  support  and  is 
inserted  into  the  model.  It  contains  six  internal  strain-gages,  which 
measure  forces  and  moments  on  the  model,  and  it  transfers  this  data  to 
a  computer  in  the  form  of  an  electrical  output.  The  specifications  of 
the  0.75"  force  balance  are  shown  in  Figure  5.  The  test  procedure 
included  taking  tare  data,  such  as  model  weight  and  amplifier  zero 
settings  prior  to  each  run.  In  addition,  tunnel  total  and  static 
pressure  transducers  were  initially  calibrated  to  determine  the  psi/mV 
ratio.  Using  this  ratio  along  with  the  atmospheric  conditions,  tunnel 
velocity  could  be  calculated.  All  data  was  acquired  using  a  Digital 
PDP  11/45  computer  data  acquisition  system. 
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Figure  4.  Six  Component  Force  Balance  with  Positive  Balance  Loads. 


IV .  Discuss  ion 

In  this  section,  I  will  discuss  the  results  of  the  drag 
measurements  and  the  effects  that  the  ventral  fin  had  on  drag 
coefficient.  In  addition,  I  will  discuss  the  static  stability  results 
and  try  to  explain  why  the  ventral  fin  apparently  reduced  the  parasite 
drag  coefficient.  Some  of  the  results  obtained  did  not  agree  well  with 
what  theory  predicted  while  others  did. 

A.  Drag  Measurements 

The  resulting  reduction  in  drag  coefficient  due  to  the  ventral 


fin  can  be  seen  from  the  summary  of  alpha-sweep  data  in  Table  1  and 
Figures  6  through  15. 
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Figure  5 


Figure  5. 
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Figure  10.  CQ  Vs  C.  with  Boundary  Trips,  Clipped  Tail  and  Fin. 
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The  values  of  and  K  were  calculated  using  linear  regression 

o 

techniques  and  then  compared  to  the  values  of  determined  by  an 

o 

eyeball  curve  fit  of  the  data  points.  For  the  -5°  to  +10°  alpha  sweeps 
(runs  1,2, 3, 7,  and  8)  between  a  14  and  18  drag  count  reduction  in 

o 

exists  due  to  the  addition  of  the  ventral  fin.  One  drag  count  equals 
.0001  Cp.  Note  the  differences  in  the  value  of  K,  differences  which 
indicate  that  the  lift  distribution,  or  span  efficiency  factor,  is 
being  altered  in  some  way.  1  attributed  this  to  the  effects  of  the 
walls  of  the  wind  tunnel  interacting  with  the  wings  at  the  higher 
values  of  alpha  (6°  to  10°). 

For  this  reason,  I  recalculated  the  values  of  CD  and  K  and  got  a 

o 

smaller,  more  consistent  range  of  values  for  K  in  the  -5°  to  +5°  alpha 
region.  Since  the  purpose  of  this  project  was  to  determine  the  drag 
coefficient  at  low  to  intermediate  CL  values,  I  decided  that  this  would 
also  give  a  more  realistic  idea  of  the  drag  count  reduction.  Analyzing 
the  -5°  to  +5°  region,  the  drag  reduction  varied  between  2  and  14  drag 
counts.  A  look  at  the  drag  polars  for  all  test  runs  shows  a  lot  of 
variance  in  at  low  values  of  C^.  This  is  due  to  the  fact  that  drag 
forces  on  the  order  of  only  0.2  lbs  were  being  measured  by  the  force 
balance  which  can  handle  up  to  35  lbs  in  axial  loading.  In  other 
words,  the  balance  wasn't  being  exercised.  The  values  for  drag 
reduction  stated  in  the  two  paragraphs  above  were  for  the  T-38 
configured  with  a  normal  vertical  stabilizer.  Further  experimentation 
was  done  with  the  tail  clipped  off  at  the  top. 

Contrary  to  my  expectations  of  a  reduction  in  drag,  the  drag  count 
increased  significantly  when  the  vertical  tail  was  clipped  (runs  7  and 
8),  leading  me  to  believe  that  the  vertical  tail  somehow  reduces  the 
amount  of  separation  drag  in  the  upper  tail  area.  From  the  standpoint 
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of  mere  speculation,  reducing  separation  may  also  be  the  function  of 
the  ventral  fin  in  reducing  drag  on  the  bottom  side  of  the  T-38. 

B.  Stat ic  Stabi 1 i ty 

The  reason  for  clipping  the  tail  was  to  compensate  for  the 


increased  directional  stability  of  the  plane  due  to  the  ventral  fin. 
Increased  stability  is  not  desirable  after  a  certain  point  is  reached 
because  the  maneuverability  of  the  plane  is  degraded.  Looking  at  Table 
2,  which  summarizes  the  stability  derivatives  calculated  during 
beta-sweep  tests  (runs  5  and  6),  I  found  that  the  ventral  fin  did  cause 


an  increase  in  Cn^  as  predicted,  from  0.0038/deg  to 
Clipping  the  tail  reduced  this  value  to  0.0027/deg. 


0 . 0043/deg . 
indicating 


that 


too  much  of  the  vertical  stabilizer  was  cut  off  in  trying  to  return  the 
C  to  its  original  value. 
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The  initial  roll  stability  results  were  unexpected  because  the 
data  indicated  that  the  addition  of  a  ventral  fin  was  making  C,  more 
negative,  or  stabilizing  the  aircraft  in  the  roll  axis  --  something 
that  is  improbable  according  to  theory.  However,  data  for  runs  5,  6, 


9,  and  10  consistently  show  that  the  addition  of  the  ventral  fin  causes 

the  value  of  Cj  to  decrease  and  become  more  negative.  In  other  words, 

theory  is  being  defied.  How  or  why  this  is  happening  is  still  unclear. 

One  trend  which  causes  concern  is  the  fact  that  the  aircraft  goes  from 

being  stable  in  the  roll  axis  to  being  unstable  when  the  vertical  tail 

is  clipped.  Again,  this  indicates  that  the  tail  was  clipped  too  much, 

just  as  the  C  data  indicated. 
nti 

Upon  analysis  of  the  data  at  zero  alpha  and  beta,  I  found  that  the 

stability  derivatives,  C,  and  C  ,  indicated  that  the  plane  was 

‘(i 

asymmetrically  loaded  and  mounted  in  the  wind  tunnel  with  a  slight 
amount  of  negative  sideslip.  Looking  at  the  data  collected  from  all 
tests,  I  realized  that  this  was  a  consistent  trend  at  zero  beta  and 
a  1  pha . 

The  central  questions  of  this  study  on  drag  reduction  remain 
unanswered.  Why  does  the  ventral  fin  reduce  the  drag  on  the  T-38? 

What  is  the  cause-effect  relationship  between  the  fin  and  the  flow? 

Upon  analysis,  graphing,  and  reduction  of  the  data  acquired  during 
testing,  I  came  to  several  conclusions  about  the  effects  of  the  ventral 
fin  on  the  T-38's  drag  and  static  stability.  However,  the  reasons  that 
some  of  these  effects  occur,  especially  the  reduction  in  drag,  are 
uncertain.  Even  after  having  talked  to  several  aer od ynami c i s t s ,  I  can 
at  best  simply  speculate  as  to  why  this  ventral  fin  causes  a  drag 
reduction.  However,  I  am  certain  that  this  is  a  local  phenomenon 
present  only  in  the  tail  region.  The  ventral  fin  is  either  reducing 
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the  amount  of  separation  in  the  tail  area  or  moving  the  point  of 
separat i on  back. 

C.  Theor i es /Specu 1  a t i ons 

To  further  investigate  the  phenomenon  and  try  to  discover  what 
was  happening,  oil  flow  tests  were  performed  on  the  model  in  the  same 
conditions  as  the  wind  tunnel  data  runs.  Figure  16  shows  the  results 
of  two  oil  flow  tests. 

From  the  pictures,  it  is  evident  that  no  strong  vortices  or  turbulent 
wakes  are  being  shed  by  the  ventral  fin  nor  were  any  present  before  the 
fin  was  added.  I  have  speculated,  however,  that  a  slight  portion  of 
the  wake  from  the  tip  of  the  fin  may  be  interacting  with  the  wake  of 

the  body  causing  the  airflow  to  separate  further  aft  and  the  drag  to 

decrease.  My  reasoning  here  is  that  although  the  fin  wake  may  nc  be 
strong  enough  to  move  the  oil  and  show  a  vortex  pattern,  it  may  be 
influencing  the  flow  in  the  tail  region.  1  would  suggest  using  smoke 
streams  to  get  a  better  visualization  of  the  phenomenon  taking  place  in 
the  tail  area  because  smoke  will  be  more  sensitive  to  small  changes  in 

the  flow  pattern.  Althougii  this  method  would  be  sufficient  as  a 

qualitative  test,  a  seven  hole  pressure  probe  or  some  type  of  external 
flow  measuring  device  should  be  used  to  gather  quantitative  data  about 
the  flow  field  in  this  region. 

The  second  speculation  1  have  formed  is  that  the  ventral  fin  is 
acting  as  a  f l ow- s t r a i gh ten i ng  device  by  preventing  spanwise  flow  from 
occurring  across  the  bottom  of  the  tail.  This  spanwise  flow  could  be 
caused  by  the  vorticies  generated  by  the  engine  inlets  or  wing-root 
intersections.  Although  the  oil  flow  tests  did  not  show  any  evidence 
of  spanwise  flow,  it  is  possible  that  weak  vortices  were  shed  that  did 
not  show  up  in  the  oil  tests.  While  the  amount  of  drag  reduction  due 
to  the  fin  is  small,  it  is  consistent  and  should  be  investigated 
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further  to  find  the  exact  nature  of  the  fin's  effect. 

V .  Cone  1  us i ons 

A.  The  ventral  fin  on  the  tail  section  of  the  T-38  was 
responsible  for  reducing  the  parasite  drag  on  the  airplane. 

B.  Lat er a  1 -d i r ec t i ona 1  stability  derivatives  are  affected  by  the 
ventral  fin.  Directional  stability  is  increased,  as  is 
lateral  stability  (which  is  contrary  to  theoretical 
predict! ons ) . 

C.  No  theory  has  been  conceived  as  to  why  the  roll  stability 
derivative  is  more  stable  with  the  addition  of  the  ventral 
fin. 

D.  The  vertical  tail  was  clipped  too  much,  causing  a  decrease  in 
directional  stability,  the  appearance  of  an  unstable  rolling 
moment  coefficient,  and  an  increase  in  parasite  drag. 

VI .  Recommendat i ons 

A.  Take  more  force  and  moment  data  to  verify  the  ventral  fin's 
potential  in  reducing  drag. 

B.  Perform  additional  testing  on  the  T-38  using  a  seven  hole 
pressure  probe  to  plot  the  flow  field  around  the  tail  region, 
with  and  without  the  ventral  fin. 

C.  Do  further  flow  visualization  studies  using  smoke  streams. 

D.  Lengthen  the  clipped  tail  or  keep  the  regular  sized  tail  in 
order  to  keep  the  stability  derivatives  close  to  normal 
T-38  values. 

S ymb  o 1 s 

Eng  1  i s h  Symbo 1 s 

A  aspect  ratio 

b  wingspan(ft) 
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K 

L 
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s 
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Greek  Symbols 


wing  mean  aerodynamic  chord,  MAC,  (ft) 

drag  force  (lbf) 

span  efficiency  factor 

dr ag-due-to-lift-f actor 

lift  force  ( lbf ) 

dynamic  pressure  (lbf/ft  ) 

Reynolds  number  (  V  /u) 

2 

wing  planform  area  (ft  ) 
velocity  (  f  t / s  ) 
total  drag  coefficient  (D/qS) 
parasite  drag  coefficient 
friction  drag  coefficient 
lift  coefficient  (L/qS) 
rolling  moment  coefficient 
yawing  moment  coefficient 


a  angle  of  attack  (deg) 

e  angle  of  sideslip  (deg) 

p  dens ity  (slug/ft^) 

u  coefficient  of  viscosity  (slug/ft  s) 
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AN  EXPERIMENTAL  INVESTIGATION  OF  AN  AIRFOIL  UNDERGOING 
LARGE  AMPLITUDE  PITCHING  MOTIONS 
3.M.  Walker*,  H.E.  Helin**,  AND  J.H.  Strickland*** 

Abstract 

Flow  visualization  and  near-surface  hot-wire  experiments  were  performed 
in  the  USAF  Academy  Aeronautics  Laboratory  2ft  x  3ft  subsonic  wind 
tunnel  on  an  airfoil  undergoing  large  amplitude  pitching  motions.  The 
experiments  were  conducted  using  a  6-in.  NACA  0015  airfoil  at  an 
airfoil  Reynolds  number  of  45000  by  pitching  from  0°  to  60°at  various 
constant  angular  pitch  rates.  The  two  cases  presented  represent  two 
different  nond imen s i ona 1  pitching  rates,  a  +  ,  equal  to  the  pitch  rate, 

<i  ,  nondimens  i  ona  1  i  zed  by  the  chord,  c,  and  the  freestream  velocity,  U*,. 

Data  for  the  two  cases,  in  which  values  of  a  +  were  equal  to  0.2  and  0.6 
show  the  dramatic  effect  of  pitch  rate  on  flow  structure.  Large  scale 
vortical  structures  are  seen  in  both  cases  but  appear  in  different  form 
at  higher  angles  of  attack  for  the  larger  a+  value.  These  structures 
are  very  energetic,  producing  reverse  flow  velocities  near  the  airfoil 
surface  of  one  to  two  and  one-half  times  the  freestream  velocity. 

I .  Introduct ion 


In  recent  years  the  study  of  unsteady  aerodynamics  has  intensified 
due  to  the  need  to  more  fully  understand  the  aerodynamics  of 
helicopters,  axial  flow  turbines  and  compressors  with  inlet 
distortions,  vertical  axis  wind  turbines,  and  missiles  and  fixed  wing 
aircraft  undergoing  rapid  maneuvers.  More  recently  the  need  to 
understand  the  unsteady  flow  behavior  around  lifting  surfaces  which 
have  significant  stalled  regions  on  their  surfaces  has  become  apparent. 

A  number  of  studies  have  been  done  with  regard  to  analytical 
approaches  to  the  "dynamic  stall"  problem  as  evidenced  by  the  reviews 
of  McCroskey  (Ref.  1  and  2).  In  general,  these  approaches  include 
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semi -emp i r i ca 1  models  (Ref.  2-6),  boundary  layer  prediction  schemes 
which  can  be  used  to  predict  some  of  the  behavior  associated  with 
unsteady  stall  (Ref.  7-9),  models  which  represent  the  separated  shear 
layers  as  discrete  vortices  (Ref.  10-13),  and  full  Nav i er -Stokes 
solutions  (Ref.  14-15).  To  date,  however,  none  of  these  analytical 
methods  have  been  proven  to  be  totally  satisfactory  for  predicting  the 
lift  and  drag  on  arbitrary  airfoil  sections  undergoing  arbitrary 
airfoil  motions  with  the  potential  for  dynamic  stall.  Part  of  this 
lack  of  credibility  is  rooted  in  the  fact  that  dynamic  stall  data, 
which  can  be  used  to  test  these  techniques,  are  still  somewhat  limited. 
In  addition,  the  lack  of  data  may  be  hiding  a  certain  amount  of  the 
"physics"  of  the  phenomenon  which  should  be  included  in  the  less 
sophisticated  approaches  to  make  them  more  reliable. 

Most  of  the  experimental  data  obtained  for  unsteady  separated  flow 
over  airfoils  have  been  obtained  from  oscillating  airfoils  undergoing 
relatively  small  sinusoidal  pitch  oscillations  (1°  to  10°)  about  a 
relatively  low  mean  angle  of  attack  (0°  -15°)  as  typified  by  the 
experiments  reported  by  McCroskey  and  Philippe  (Ref.  16),  McAlister  and 
Carr  (Ref.  17),  Martin,  et  al  (Ref.  18),  and  Robinson  and  Luttges  (Ref. 
19).  The  types  of  data  obtained  in  these  studies  include  flow 
visualization  data,  hot-film  and  hot-wire  data,  and  surface  pressure 
measurements.  Such  airfoil  motions  are,  of  course,  applicable  to  many 
of  the  fluid  devices  studied  in  the  past.  On  the  other  hand,  few 
experimental  data  are  available  for  situations  in  which  an  airfoil 
undergoes  very  large  angle  of  attack  excursions  taking  it  into  deep 
dynamic  stall.  Applications  such  as  the  recently  conceived 
"supermaneuverability"  of  fighter  aircraft  (Ref.  20)  require  a  more 
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thorough  understanding  of  dynamically  stalled  airfoils  at  angles  of 
attack  which  may  exceed  45°.  In  addition,  the  motion  of  the  airfoil  for 
this  application  will  perhaps  be  more  closely  related  to  a  constant 
motion  as  opposed  to  a  sinusoidal  motion. 

A  very  limited  amount  of  experimental  data  has  been  obtained  for 
airfoils  undergoing  constant  pitching  rate  motions  up  to  moderate 
angles  of  attack  of  at  least  30°.  These  works  include  the  study  of 
Harper  and  Flanigan  (Ref.  21),  who  obtained  force  balance  data  on  a 
small  aircraft  model  pitching  up  to  30°,  the  work  of  Ham  and  Garelick 
(Ref.  22),  who  obtained  surface  pressure  measurements  on  an  airfoil 
pitching  up  to  30°,  and  the  work  of  Francis,  et  al  (Ref.  23),  who 
obtained  surface  pressure  measurements  on  an  airfoil  pitching  up  to  60°. 
None  of  these  works  have  included  any  flow  visualization  data  or  hot 
wire  data.  The  flow  visualization  data  and  surface  hot  wire  data 
presented  in  this  paper  are  intended  to  partially  fill  this  void. 

I  I .  Experimental  Arrangement 

All  data  were  obtained  in  the  USAF  Academy  Aeronautics 
Laboratory's  low-speed  2ft  x  3ft  subsonic  wind  tunnel  on  a  6in.  chord 
NACA  0015  airfoil.  For  this  particular  series  of  experiments  the 
freestream  velocity  was  maintained  between  19.0  and  20.0  ft/sec  which 
yields  a  Reynolds  number  based  on  the  airfoil  chord  of  approximately 
45000.  A  computer  controlled  pitch  oscillator  described  by  Francis,  et 
al  (Ref.  24)  was  used  to  impart  constant  pitch  rate  motions  to  the 
airfoil  between  0°  and  60°  angles  of  attack.  For  the  two  cases  studied 
herein  u  was  held  constant  at  440°/sec  and  13S0°/sec  which  yielded  values 
of  u+ 5  .ic/Uro  equal  to  0.2  and  0.6  respectively. 

The  basic  flow  visualization  scheme  makes  use  of  a  "smoke  wire" 
stretched  across  the  tunnel  normal  to  the  span  and  upstream  of  the 


i  ie 


USAFA-TR-85-2 


pitching  airfoil.  A  smoke  producing  oil  (theoretical  fog  fluid)  was 
coated  on  a  .005  in.  diameter  tungsten  wire  which  was  in  turn  heated 
electrically  to  produce  a  large  number  of  fine  smoke  streaks  in  the 
flow.  These  streaks  were  spaced  fairly  uniformly  due  to  the  regular 
spacing  of  oil  droplets  formed  when  the  wire  was  coated.  The  smoke  was 
illuminated  by  a  high  intensity  strobe  light  placed  downstream  of  the 
airfoil.  A  PDP  11  computer  was  programmed  to  control  the  smoke  wire, 
the  airfoil  pitching  motion,  and  the  strobe  flash  in  the  proper  timing 
sequence.  A  35mm  camera  aimed  parallel  to  the  pitch  axis  was  used  to 
record  the  visual  data. 

In  order  to  obtain  near  surface  velocity  measurements,  seven 
hot-wires  were  mounted  on  the  upper  surface  (suction  side)  of  the 
airfoil  as  shown  in  Figure  1.  The  hot-wire  sensing  elements  (TSI-10 
hot  films)  were  soldered  to  pairs  of  number  9  sewing  needles  which 
protruded  above  the  airfoil  surface  approximately  3  percent  (0.2  in.) 
of  the  chord  length.  The  needle  supports  in  turn  were  mounted  in 
electrically  insulated  plugs  machined  flush  with  the  airfoil  surface. 
The  probes  were  located  along  the  chord  line  at  7,  20,  33,  47,  60,  78, 
and  87  percent  chord.  Velocity  signals  were  obtained  from  the  probes 
using  a  TS I  model  1050  hot-wire  anemometer  system  at  a  sampling  rate  of 
1  KHz.  Twenty  five  repetitions  of  each  case  were  run  to  obtain 
ensemble  averages  of  the  velocity  signal  from  each  probe. 


USAFA-TR-85-2 


0.4  in  6  spaces  at  0.8  in  -  4.8  in  0.8  in 


Fig.  1  Surface  Hot-Wire  Configuration 

III.  Experimental  Results 

The  instantaneous  angle  of  attack,  a,  of  the  airfoil  as  a  function 
of  time  is  shown  in  Figure  2  for  the  two  pitching  rates  reported 
herein.  The  increase  of  a  with  time  is  very  linear  for  a*  =0.6  and 
fairly  linear  for  a+=0.2,  the  maximum  deviation  being  2%  and  1% 
respect i ve 1 y . 


0  10  20  30  40  50  60  70  60  90 

a t (degrees) 

Fig.2  Variation  of  Pitch  Angle  With  Time 
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A.  Flow  Visualization 

Flow  visualization  data  are  shown  in  Figures  3  and  4  at 
selected  angles  of  attack  for  the  two  pitching  rates.  The  differences 
in  flow  structure  are  rather  dramatic  at  similar  angles  of  attack  for 
the  two  cases.  Data  at  the  lower  a+  value  in  Figure  3  show  a 
separation  bubble  beginning  to  form  at  the  leading  edge  at  an  angle  of 
attack  of  about  20°.  This  bubble  quickly  evolves  into  a  leading  edge 
vortex.  This  vortex  increases  in  size  until  its  diameter  is  on  the 
order  of  half  the  chord  length.  While  continuing  to  grow,  the  vortex 
begins  to  move  away  from  the  upper  (suction)  surface  of  the  airfoil 
which  has  by  then  reached  a  high  angle  of  attack  (approximately  40°). 
Francis,  et  al  (Ref.  23),  whose  data  were  taken  under  slightly 
different  conditions  from  Figure  3,  suggest  that  the  maximum  lift  at 
the  lower  rate  is  produced  at  angles  of  attack  between  30°-35o.  The  flow 
visualization  data  in  turn  suggests  that  maximum  lift  occurs  when  the 
vortex  on  the  suction  side  is  well  developed  and  still  relatively  close 
to  the  surface. 

Data  at  the  higher  u+  value  in  Figure  4  shows  the  formation  of  two 
distinct  vortices  which  appear  to  be  of  comparable  importance  although 
not  necessarily  of  comparable  strength.  Initially  there  appears  to  be 
a  region  of  boundary  layer  separation  which  moves  from  the  trailing 
edge  toward  the  leading  edge.  At  an  angle  of  attack  near  30°,  a  leading 
edge  vortex  begins  to  form.  At  around  33°  the  leading  edge  vortex  has 
grown  large  enough  to  retard  the  forward  movement  of  the  trailing  edge 
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separation  zone.  At  an  angle  of  attack  near  40°,  the  forward  progress 
of  the  trailing  edge  separation  zone  is  halted  entirely  due  to 
"reattachment"  of  the  flow  caused  by  the  leading  edge  vortex. 
Furthermore,  the  clockwise  vorticity,  shed  at  the  trailing  edge  due  to 
increasing  bound  circulation,  reattaches  the  flow  at  the  trailing  edge 
and  produces  a  counter  clockwise  vortex  in  the  trailing  edge  separation 
region  similar  to  the  leading  edge  vortex.  As  the  angle  of  attack 
increases,  both  the  leading  edge  and  trailing  edge  vortices  grow  to 
nearly  a  half  chord  in  diameter  and  remain  relatively  stable  in  their 
positions  until  the  airfoil  reaches  its  maximum  angle  of  60°.  McAlister 
and  Carr  (Ref.  17)  have  noted  similar  behavior  on  airfoils  undergoing 
harmonic  oscillations  although  the  effects  were  much  less  pronounced 
due  to  the  low  angles  of  attack  in  their  study.  They  denote  the 
leading  edge  vortex  as  the  "dynamic  stall  vortex"  and  the  trailing  edge 
vortex  as  the  "shear  layer  vortex."  They  also  note  that  the  "shear 
layer  vortex"  seems  to  disappear  at  high  Reynolds  numbers,  leaving  only 
the  "dynamic  stall  vortex."  The  boundary  layer  separation  and  the 
disappearance  of  the  trailing  edge  vortex  at  high  Reynolds  numbers 
simply  indicates  that  the  flow  does  not  separate  in  the  trailing  edge 
region  prior  to  the  formation  of  the  leading  edge  vortex. 
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B.  Surface  Velocity  Data 

Selected  "near-surface"  velocity  data  are  shown  in  Figures  5 
and  6  for  the  corresponding  cases  given  in  Figures  3  and  4 
respectively.  The  primary  purpose  of  the  data  taken  from  the  seven 
hot-wire  probes  was  to  obtain  insight  into  the  strength  of  large  scale 
vortices  on  the  suction  side  of  the  airfoil  by  measuring  the  reverse 
flow  near  the  wall  (3%  chord  away  from  the  wall).  The  direction  of  the 
flow  is  not  apparent  from  examination  of  the  single  element  hot-wire 
data  without  the  aid  of  the  flow  visualization  data.  In  some  cases, 

even  comparison  of  both  types  of  data  does  not  allow  one  to  make  a 

definitive  decision  concerning  flow  direction.  In  addition,  the  data 
are  inaccurate  for  cases  in  which  the  velocity  is  close  to  zero  due  to 
the  high  turbulence  intensity  levels  under  those  conditions.  However, 
significant  forward  and  reverse  flow  velocities  are  easily  dectectable 
and  are  measurable  with  reasonable  accuracy. 

For  the  a+  value  of  0.2,  maximum  reverse  flow  velocities  on  the 
order  of  140%  of  freestream  were  measured.  The  maximum  reverse  flow 
velocity  occurs  directly  under  the  center  of  a  vortex.  The  arrows  in 
Figure  5  represent  the  approximate  chordwise  locations  of  the  center  of 
the  vortex  as  it  moves  along  and  away  from  the  airfoil.  At  an  angle  of 

attack  of  24®  a  suction  peak  and  regions  of  leading  edge  separation, 

reattachment,  and  trailing  edge  separation  are  indicated  in  both  Figure 
3a  and  5a.  At  an  u  of  27°,  Figures  3b  and  5b  show  leading  edge 
separation,  a  probable  reverse  flow  region,  and  a  reattachment  region. 
The  small  reverse  flow  velocity  obtained  at  sensor  4  is  apparently  due 
to  the  fact  that  the  hot-wire  is  sticking  up  somewhat  into  the  vortex 
core.  Figures  3c  and  5c  show  leading  edge  separation,  a  dramatic 
reverse  flow  region,  and  a  short  reattachment  zone  at  i  =  36°.  Figures 
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3d,e,f,  and  5d,e,f  show  the  vortex  convecting  downstream  with  the 
subsequent  reduction  in  reverse  flow  velocities  near  the  surface.  It 
is  interesting  to  note  that  since  the  airfoil  is  still  pitching  up, 
this  vortex  actually  moves  slightly  forward  in  a  reference  frame  fixed 
to  airfoil. 

For  the  a+  value  of  0.6,  maximum  reverse  flow  velocities  on  the 
order  of  210%  of  the  freestream  velocity  are  observed.  As  before,  the 
peak  reverse  flow  velocities  occur  under  the  centers  of  the  vortices, 
and  the  arrows  in  Figure  6  note  the  chordwise  locations  of  the  vortex 
centers.  Figure  6a  shows  what  would  appear  to  be  a  typical  surface 
velocity  distribution  for  flow  over  an  airfoil  in  which  separation  is 
just  beginning  to  occur  at  the  rear  of  the  airfoil.  It  should  be 
noted,  however,  that  under  steady  flow  conditions  this  situation  would 
occur  at  angles  of  attack  of  less  than  10“  whereas  in  the  present  case 
the  angle  of  attack  is  28®.  Figures  4b  and  6b  show  the  development  of 
the  two  vortices  as  well  as  the  stagnation  or  reattachment  zone  between 
them.  Note  that  the  leading  edge  vortex  builds  rather  rapidly  and 
initially  produces  a  much  higher  reverse  flow  velocity  than  does  the 
trailing  edge  vortex.  By  the  time  that  the  airfoil  pitches  to  its 
maximum  angle  of  60°,  however,  the  trailing  edge  vortex  produces  a 
reverse  flow  velocity  of  the  same  order  as  the  leading  edge  vortex. 
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a)  or  -  24°  d)  ot  =  40° 


b)  o  *  27°  e)  ol  =  47° 


c)  <*  -  36°  0  a  -  59° 


Fig. 5  Absolute  Velocity  Ratio  jul/u*  Obtained  From  Surface 
Hot-Wires  Re-45,000,  a+0.2 
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IV .  Cone  1  us i ons 

From  the  work  presented  herein,  it  can  be  seen  that  the  flow 
structure  is  a  strong  function  of  the  nond imens i ona 1  pitching  rate,  u+ . 

The  higher  the  value  of  u+,  the  higher  the  angle  of  attack  reached 
before  the  beginning  of  flow  separation,  the  more  energetic  the  suction 
peak,  the  more  energetic  the  leading  edge  vortex,  and  the  higher  the 
reverse  flow  velocities  near  the  surface  of  the  airfoil.  In  addition, 
at  the  higher  u+  value,  significant  secondary  vortical  structures 
appear.  For  the  two  cases  presented,  u+  =0.2  and  u  =0.6,  large  reverse 
flow  velocities  on  the  order  of  1  to  2  times  the  freestream  velocity, 
respectively,  near  the  suction  surface  of  the  airfoil,  indicate  the 
relative  magnitudes  of  vortical  energy  present.  Maximum  reverse  flow 
velocities  occur  for  any  given  angle  of  attack  directly  under  the 
center  of  the  vortex.  Other  parameters  which  were  not  varied  in  this 
study,  but  which  may  also  significantly  affect  the  flow  structure,  are 
the  pitch  axis  location  (maintained  at  25%  chord  in  the  present  study), 
Reynolds  number,  and  initial  angle  of  attack. 
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INTERRELATED  EFFECTS  OF  PITCH  RATE  AND  PIVOT  POINT  ON  AIRFOIL 

DYNAMIC  STALL 

Hank  E.  He  1  i n *  and  3ohn  M.  Walker** 

Abstract 

Experimental  investigations  were  conducted  to  study  energetic 
dynamic  stall  vortices  and  the  associated  unsteady  aerodynamics 
generated  by  a  pitching  NACA  0015  airfoil.  The  airfoil  model  was 
pitched  from  0°  to  60°  at  constant  rates  of  460°/s,  920°/s,  and  1380°/s 
about  its  quar t er -chor d ,  half  chord, and  three  quarter  chord  positions. 
Extensive  35mm  still  photographs  and  16mm  high-speed  movies,  both 
employing  smoke  wire  flow  visualization,  documented  the  initiation  and 
development  of  the  time  dependent  dynamic  stall  phenomena.  In 
addition,  hot-wire  anemometry  measurements  were  made  which  provided  for 
more  quantitative  analysis  of  the  unsteady  separated  flowfields.  Pitch 
rate  and  pivot  point  were  shown  to  have  interrelated  effects  on  the 
development  of  the  dynamic  stall  flowfield.  In  many  cases  similar 
"looking"  flowfields  were  generated  by  different  combinations  of  pitch 
rate  and  pivot  point.  However,  significant  differences  were  observed  in 
the  near-surface  velocity  profiles. 

I  .  I n  t  r  oduc  t i on 


The  energetic  nature  of  the  unsteady  flowfields  generated  by 
pitching  airfoils  has  been  a  topic  of  study  for  much  of  the  20th 
Century  (Ref.  1).  The  vast  number  and  interaction  of  variables 
involved  have  significantly  complicated  both  theoretical  and 
experimental  efforts  to  understand  the  fluid  mechanics  and  the 
development  of  reliable  prediction  methods.  Most  theoretical  studies 
revolve  around  attempts  to  relax  assumptions  postulated  in  steady  thin 
airfoil  theory.  In  general,  the  results  from  these  studies  have  been 
somewhat  restricted  and  significant  progress  has  been  elusive.  This  is 
due  in  part  to  the  lack  of  broad  scope  experimental  data  which  are 
needed  to  set  modeling  parameters  and  as  a  basis  for  checking 
theoretical  calculations.  The  general  thrust  of 
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experimental  unsteady  aerodynamic  research  has  been  to  reduce  the 
undesirable  effects  associated  with  dynamic  stall.  Many  of  the 
parameters  involved  have  been  restricted  to  model  specific  areas  and 
types  of  problems.  However,  with  the  more  recent  realization  that  the 
extremely  energetic  nature  of  unsteady  flows  might  be  exploited  and 
utilized  to  enhance  performance,  the  scope  of  research  efforts  has  been 
expanded  significantly  (Ref.  2). 

Unsteady  aerodynamic  effects  generated  by  pitching  or  oscillating 
airfoils  are  generally  classified  under  the  heading  of  "Dynamic  Stall." 

Inherent  in  this  phenomenon  is  the  development  of  a  dynamic  stall 
vortex  which  occurs  as  the  lifting  surface  dynamically  surpasses  its 
static  stall  angle  of  attack.  Large  unsteady  aerodynamic  forces  are 
generated,  from  which  the  lift,  drag,  and  moment  coefficients  greatly 
exceed  their  maximum  static  counterparts  (Ref.  3).  The  unsteady 
effects  of  dynamic  stall  are  dominated  by  turbulent  flow  and  the 
production  of  large  scale  vortices.  The  more  general  methods  employed 
in  analyzing  dynamic  stall  involved  tests  with  airfoils  driven  through 
moderately  large  amplitude,  slow  oscillations  in  angle  of  attack.  This 
was  consistent  with  attempts  to  understand  and  eliminate  the 
undesirable  effects  associated  with  dynamic  stall  on  helicopter  rotors. 

In  contrast,  recent  efforts  are  now  also  exploring  the  possibilities 
of  utilizing  the  unsteadiness  of  the  flowfield  to  enhance  performance. 
These  studies  have  examined  new  parameter  combinations  including  much 
higher  osc i 1 1  a t i on /p i t ch  rates  (Ref.  2).  Studies  have  been  performed 
to  analyze  the  repetitive  interaction  of  the  dynamic  stall  vortices  as 
a  means  of  maintaining  flow  attachment  (and  hence,  greatly  increasing 
lift)  at  high  angles  of  attack  (Ref.  <0  .  In  addition,  extensive 
studies  have  also  been  done  to  correlate  this  phenomenon  as  a  function 
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of  the  driving  parameters  involved,  that  is,  pivot  location,  airfoil 
shape,  Reynolds  number,  and  mean  angle  of  attack  (Ref.  5). 

Large,  energetic  vortical  structures  have  been  shown  to  be  an 
important  correlate  of  oscillating/pitching  airfoils,  and  current 
research  activities  have  documented  their  dramatic  impact  on  airfoil 
lift  and  drag  characteristics  (Ref.  6  and  7).  In  addition,  the  impetus 
to  exploit  the  energetic  nature  of  large  vortices  as  a  potential  to 
enhance  performance  has  already  been  demonstrated  (Ref.  1).  Clearly, 
before  such  a  realistic  utilization  is  possible,  extensive  studies  must 
be  performed  to  expand  our  knowledge  concerning  fundamental  aspects. 
Questions  regarding  the  direct  relationship  between  the  time  dependent 
fluid  dynamics  and  associated  a i r f o i l / f 1 owf i e 1 d  interaction,  have  yet 
to  be  answered.  Before  attempts  to  utilize  these  phenomena  are  made, 
the  flowfield  dependence  upon  the  driving  parameters  of  airfoil 
geometry  and  dynamics  must  be  fully  understood.  The  present  study 
focuses  on  only  a  small  part  of  this  problem:  the  relationships  between 
the  pitch  rate,  pitch  axis,  and  resultant  vortex  development. 

I  I .  Methods 

Experiments  were  conducted  using  a  NACA  0015  airfoil  model  in  the 
USAF  Academy's  subsonic  wind  tunnel.  Two  experimental  techniques  were 
employed:  I)  smoke-wire  flow  visualization  using  35mm  still  and  16mm 
high-speed  movie  cameras,  and  2)  near-surface  hot-wires  mounted  in  a 
staggered  array  on  the  airfoil's  surface. 

The  airfoil  model  was  a  NACA  0015  extruded  aluminum  section  with  a 
6"  chord  and  22"  span.  The  airfoil  was  pitched  from  0°  to  60°  by  a 
s i ng 1 e -deg r ee -o f - f r eedom  oscillator  employing  a  DC  stepping  motor  and  a 
computer  driven  index  controller.  Airfoil  angular  position  was 
measured  by  a  linear  potentiometer  mounted  on  the  oscillator.  System 
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performance  checks  documented  cycle  to  cycle  repeatability  with  a 
maximum  standard  deviation  of  0.5°  over  a  25  run  ensemble  average. 

The  USAF  Academy's  subsonic  wind  tunnel  has  a  2'x3'  test  section. 
All  testing  was  conducted  at  a  f  reestream  velocity  of  20  ft/s,  which, 
with  the  6"  NACA  0015,  resulted  in  an  airfoil  Reynolds  number  of 
45,000.  The  turbulence  intensity  level  of  the  tunnel  in  this  velocity 
range  has  been  documented  at  less  than  0.1%.  The  low  turbulence  level 
and  the  low  velocity  were  necessary  for  high  quality  detailed  flow 
visual i za t i on . 

Flow  visualization  was  obtained  using  a  .005"  tungsten  wire  coated 
with  theatrical  fog  fluid.  An  electrical  current  was  applied  to  the 
wire,  which  produced  a  set  of  fine  streaklines  across  the  test  section. 

These  streaklines  were  photographed  using  a  high  intensity  arc-lamp 
strobe  light,  with  a  flash  duration  of  7  ps,  coupled  with  a  35mm  camera 
and  a  16mm  high-speed  movie  camera.  The  35mm  photographs  are 
phase-locked  single  exposures  on  400  ASA  TRI-X  film  developed  at  1200 
ASA.  Only  one  photograph  could  be  taken  during  a  pitch  motion,  a 
procedure  which  required  multiple  motion  sequences  with  each  successive 
photograph  at  an  incremented  angle  of  attack.  The  16mm  high-speed 
camera  system,  however,  allowed  multiple  photographs  over  one  motion 
cycle,  which  provided  for  detailed  studies  of  the  flowfield  evolution 
as  a  function  of  time.  The  entire  osc i I  1  a t or /smoke-wi r e/ f J as h/camer a 
system  was  controlled  and  synchronized  using  a  PDP  11/45  and  PDP  11/03 
computer  system. 

Near  surface  velocity  measurements  were  made  using  seven  hot-wires 
mounted  on  the  upper  (suction)  surface  of  the  airfoil.  The  hot-wire 
sensing  elements  (TSI-10  Hotfilms)  were  fixed  at  .2"  above  the  surface 
and  equally  spaced  at  .8"  increments.  The  sensor  positions  were:. 4", 
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1.2",  2.0",  2.8",  3.6",  4.4",  and  5.2"  on  the  6"  chord  airfoil. 

Velocity  signals  were  obtained  using  a  TS I  1050  anemometer  system 
coupled  with  the  Laboratory  Peripheral  System  package  of  the  PDP  11/45 
computer.  The  velocities  at  the  various  sensor  locations,  with  the 
airfoil  at  0°  angle  of  attack,  were  compared  with  a  freestream  hot-wire 
to  obtain  normalizing  coefficients.  The  measured  velocities  over  the 
airfoil  surface  were  in  turn  normalized  to  a  percent  of  freestream,  for 
analysis.  The  sampling  rate  varied  with  the  pitch  rate  employed  but 
was  always  above  l  KHz. 


0.8  in 


Fig.  1  Surface  Hot-Wire  Configuration 


III.  Resu 1 ts 

A  well  known  effect  of  oscillating  or  pitching  an  airfoil  is  that 
the  static  stall  angle  can  be  exceeded  by  some  amount  before  dynamic 
stall  occurs.  Dynamic  stall  is  preceded  by  the  separation  of  flow  near 
the  trailing  edge.  This  forms  a  separated  region  over  the  airfoil 
which  grows  rapidly  as  its  initiation  point  moves  up  from  the  trailing 
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edge  toward  the  leading  edge.  The  initiation,  development,  and  even 
the  strength  of  the  dynamic  stall  vortex,  have  been  shown  to  be  a 
complex  function  of  the  driving  parameters.  This  function  is  far  from 
being  fully  understood.  Rather  than  complicate  the  scheme  further  by 
using  periodic  motions,  which  set  up  a  hysteresis  loop,  constant  pitch 
motions  from  0°  up  to  an  angle  of  attack  of  60°  were  employed.  Once  the 
maximum  angle  was  reached,  the  airfoil  was  stopped,  allowing  the 
unsteady  flowfields  to  develop  freely. 

A  steady  state  study  of  static  stall  conditions  for  the  NACA  0015 
airfoil  was  conducted.  As  the  airfoil's  angle  of  attack  was  slowly 
incremented  toward  the  static  stall  angle  under  steady  flow  conditions, 
a  trailing  edge  separation  zone  formed  and  moved  up  the  suction  surface 
of  the  airfoil  until,  at  an  angle  of  attack  of  12°,  flow  was  separated 
over  fifty  percent  of  the  airfoil  chord.  This  was  readily  determined 
by  smoke  flow  visualization.  This  steady  flow  stall  criterion  will  be 
used  to  compare  the  delay  in  onset  of  dynamic  stall  with  regard  to 
angle  of  attack  and  the  initiation  of  the  dynamic  stall  vortex. 
Qualitative  analyses  of  flow  visualization  are  based  upon  the 
initiation  of  the  dynamic  stall  vortex,  its  temporal  position  with 
respect  to  the  airfoil  angle  of  attack  --  which  is  dependent  on  the 
pitch  rate  employed  --  and  its  spatial  position  with  respect  to  the 
airfoil  chord.  These  analyses  are  compared  and  correlated  with  the 
near  surface  velocities  measured  with  the  hot-wires. 

Table  1  shows  the  combination  of  pitch  rates  and  pitch  axes  that 
were  employed.  These  combinations  of  driving  variables  allowed  for 
analyses  of  the  initiation  and  development  of  unsteady  dynamic  stall 
flowfields  based  upon  chord  angular  displacement  with  respect  to  time, 
airfoil  leading  edge  angular  velocities,  and  the  interrelated  effects 


136 


USAFA-TR-85-2 


of  chord  and  leading  edge  angular  velocities. 

PITCH  RATE 

460  ° / s  920°/ s  1  380  ° / s 

0.25c  XXX 

0.5c  X 

0.75c  X 

Table  1.  Pitch  Rates  and  Pitch  Axes 

A.  Increased  Pitch  Rate,  Fixed  Pitch  Axis 

With  the  pitch  axis  fixed  at  0.25  c,  the  dramatic  effects  of 
increased  pitch  rate  are  visually  documented  in  Figure  2.  As  the  pitch 
rate  was  increased  from  460°/s  to  920°/s,  and  finally  to  !380°/s,  the 
initiation  of  the  dynamic  stall  vortex  was  delayed  to  higher  angles  of 
attack.  This  delay  is  not  a  linear  function  of  the  pitch  rate.  The 
effects  of  an  increment  in  pitch  rate  are  more  prominent  at  lower  rates 
than  at  higher  rates.  At  a  pitch  rate  of  460°/s,  the  dynamic  stall 
vortex  began  to  form  around  a  =22°.  At  double  the  pitch  rate,  920°/s,  it 
began  forming  at  a  =31°.  However,  an  increase  to  1380°/s  only  pushed  the 
dynamic  stall  vortex  initiation  to  a  =34°.  Another  dominant  feature  is 
the  spatial  position  of  the  vortex  with  respect  to  the  airfoil.  At  the 
lower  pitch  rates  the  vortex  formed  much  earlier  in  the  pitch  cycle  and 
moved  rapidly  downstream.  As  the  pitch  rate  increased,  the  dwell  time, 
that  is,  the  time  that  the  vortex  remained  close  to  and  over  the 
airfoil,  also  increased.  In  all  cases,  the  dynamic  stall  vortex 
remained  energetic  and  continued  to  rotate  until  well  downstream  of  the 
airfoil.  This  is  not  readily  discernible  from  still  photographs  taken 
over  multiple  cycles  but  was  clearly  evident  in  the  high-speed  16mm 
movies  taken  over  a  single  cycle. 


h)  45  .  1380  /s 


g)  30°,  1380°/s 


Fig.2  continued 
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The  16mm  movies  taken  at  500  frames  per  second  visualized  the 
dramatic  differences  between  streaklines  and  streamlines  for  a  time 
dependent  flowfield.  Many  of  the  contorted  streaklines  developed  early 
as  the  vortex  was  forming,  and  then  were  carried  along  as  historical 
patterns.  However,  the  energetic  rotation  of  the  vortex  was  easily 
seen  from  frame  to  frame  and  the  streaklines  are  indicative  of  the 
vortex  position  and  temporal  development.  In  addition,  the  movies 
document  the  direction  of  the  flow  at  different  points  in  the 
flowfield.  This  was  very  useful  in  determining  rotational  directions 
and  reverse  flow  areas  over  the  surface. 

When  the  movies  are  compared  to  the  still  photographs  and  hot-wire 
data,  with  respect  to  angle  of  attack,  a  quantitative  analysis  of  the 
flowfield  and  vortex  development  is  easily  made.  Since  the  hot-wire 
sensors  are  close  to  the  surface,  they  indicate  the  magnitude  of  the 
flow  which  is  essentially  parallel  to  the  surface.  And  when  correlated 
with  the  flow  visualization,  the  direction  of  the  flow  along  the 
surface  is  readily  discernible  much  of  the  time.  Peaks  in  the  profiles 
match  with  the  spatial  positions,  as  observed  in  the  flow  visualization 
of  the  vortex  centers  over  the  airfoil.  Figure  3  shows  the  velocities 
in  percent  freestream  at  the  first  sensor  location  for  the  three 
different  pitch  rates.  Peaks  in  the  velocity  profiles  at  the  first 
sensor  agree  with  the  vortex  initiation  angles  documented  in  the  flow 
visualization.  Interestingly,  up  to  the  point  of  vortex  initiation, 
the  slope  of  the  velocity  profiles,  as  a  function  of  the  angle  of 
attack,  is  relatively  independent  of  the  pitch  rate.  Of  further 
interest,  the  maximum  velocity  of  the  first  sensor  was  also  independent 
of  the  pitch  rate  as  evidenced  by  the  near  equal  magnitudes  of  the 
peaks  in  Figure  3.  However,  the  relative  magnitudes  of  the  velocities 
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at  stations  further  aft  are  different  for  changes  in  pitch  rate. 

Correlation  of  peaks  in  the  velocity  profiles  with  observed  vortex 
positions  from  the  flow  visualizations  showed  that  the  maximum 
velocities  increase  as  the  pitch  rate  was  increased.  An  example  of 
this  is  in  Figure  4,  which  shows  the  magnitude  of  the  velocity  at 
sensor  number  2  on  the  suction  surface  of  the  airfoil  for  two  different 


pitch  rates.  These  higher  velocities  and  the  fact  that  the  flow 
visualization  also  showed  more  cohesive  appearing  vortices  as  the  pitch 
rate  was  increased,  are  a  direct  inference  that  vortex  strength 
increases  as  a  function  of  increasing  pitch  rate.  This  readily 
suggests  that  the  increased  dwell  time  at  the  higher  pitch  rates  is 
directly  coupled  to  the  increased  energetic  and  cohesive  nature  of  the 
dynamic  stall  vortices  at  the  higher  pitch  rates. 


Fig.4b  Magnitude  of  the  Velocities  over  the  Airfoil  Surface,  in 
Percent  Freestream,  at  Sensor  Number  Two.  Pitch  Axis  is  .25c. 
Pitch  Rate  is  a)  920°/s,  b)  1380°/s. 
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B.  Increased  Pitch  Axis,  Pitch  Rate  Constant 

As  shown  by  the  flow  visualization  photographs  in  Figure  5, 
the  effects  of  moving  the  pitch  axis  further  along  the  chord  are 
similar  to  those  discovered  for  increasing  pitch  rate.  As  the  distance 
from  the  leading  edge  was  increased  from  0.25  c  to  0.5  c  and  finally 
0.75  c,  the  onset  of  dynamic  stall  was  again  delayed.  Once  the  leading 
edge  vortex  was  generated,  its  evolution  and  growth  characteristics 
during  the  remaining  portion  of  the  pitch  cycle,  and  with  time, 
appeared  constant. 

Although  the  initiation  of  the  dynamic  stall  vortex  was  delayed  to 
very  high  angles  of  attack,  the  subsequent  rapid  development  and 
movement  of  the  vortex  over  the  airfoil  was  very  similar,  with  no 
appreciable  increase  in  dwell  time  regardless  of  the  pitch  axis 
location.  This  is  in  contrast  to  the  coupled  effect  of  delayed 
initiation  and  increased  dwell  time  as  pitch  rate  was  increased.  In 
addition,  there  is  not  a  perceivable  increase  in  coherent  appearance 
with  increased  pitch  axis  distance  as  there  was  with  increased  pitch 
rate  (Figures  2  and  5).  Again,  the  high-speed  I6mm  movie  documented 
the  energetic  nature  of  the  dynamic  stall  flowfield,  showing  continued 
growth  and  significant  rotation  throughout  the  time  the  vortex  was  in 
the  vicinity  of  the  airfoil. 

Dynamic  stall  at  the  first  hot  wire  sensor,  evidenced  by  a  sharp 
decrease  in  the  magnitudes  of  the  velocities  shown  in  Figure  6,  occurs 
at  30°,  36°,  and  42°  for  pitch  axies  of  0.25  c,  0.5  c  and  0.75  c, 
respectively.  This  represents  a  relatively  linear  increase  in  stall 
delay  as  a  function  of  increased  pitch  axis  distance.  Figure  6  is  also 
illustrative  of  significant  differences  in  the  velocity  magnitude 
profiles  at  the  first  sensor  location.  There  is  a  definite  decrease  in 
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the  slope  with  respect  to  the  angle  of  attack  as  the  pitch  axis 
distance  from  the  leading  edge  was  increased.  Also,  the  peak  velocity 
at  the  first  sensor  decreased  as  pitch  axis  was  increased  --  more  so 
f r om  0.25  c  to  0.5  c  than  from  0.5  c  to  0.75  c.  These  two  phenomena 
are  completely  different  from  observations  at  increasing  pitch  rate 
(Figure  3 ) . 

As  mentioned  in  the  flow  visualization  discussion,  the  development 
of  the  stall  flowfield  appears  quite  similar,  except  for  an  increase  in 
the  dynamic  stall  angle,  regardless  of  the  pitch  axis  employed.  The 
remarkable  similarities  in  the  velocity  profiles  as  a  function  of  time, 
at  sensor  number  three,  are  shown  in  Figure  7  for  the  three  pitch  axes. 

As  pitch  axis  distance  was  increased,  the  initial  slopes  decreased 
slightly  (noted  before  at  sensor  number  one),  the  maximum  magnitude  of 
the  peaks  decreased  slightly,  and  dynamic  stall  was  delayed.  However, 
the  general  shapes  of  the  profiles  are  overwhelmingly  similar  and  are 
merely  shifted  an  increment  of  time  as  the  pitch  axis  is  moved  aft. 
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c)  60°,  .25c 


e)  45  ,  .5c 


d)  30°,  .5c 


f)  60°,  .5c 


Fig. 5  Pitch  Axes  of  0.25  c,  a  0.5  c  and  0.75  c  with  Pictures  at 
30°,  45°and  60°  for  each  axes.  Pitch  Rate  Constant  at  920°/s. 
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Fig.6  Magnitude  of  the  Velocities  at  Sensor  Number  One,  in 
Percent  Freestream,  for  Pivot  Axes  of  0.25  c,  0.5  c  and  0.75  c 
Pitch  Rate  is  920°/s. 
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IV .  Discuss  ion 

Studies  of  an  airfoil  pitched  from  0°  to  60°  angle  of  attack  at 
various  pitch  rates  about  different  pitch  axes  have  shown  the 
development  of  a  complex  interactive  flowfield.  The  most  prevalent 
feature  was  the  formation  of  an  energetic  dynamic  stall  vortex  which 
was  both  spatially  and  temporally  dependent  on  the  driving  parameters. 
The  dynamic  stall  phenomenon  was  extremely  reproducible.  In  light  of 
this  reproducible  nature,  single  shot  photographs,  rather  than  motion 
cycles,  and  ensemble  averaged  hotwire  profiles  were  used  to  quantify 
flowfield  development.  In  addition,  high  speed  movies  documented  the 
energetic  rotation  and  travel  of  the  dynamic  stall  vortex  as  it 
developed.  The  presence  of  the  dynamic  stall  vortex  was  observed  to 
have  significant  effect  on  the  velocity  distribution  over  the  airfoil 
surface  at  angles  of  attack  well  beyond  the  normal  static  stall. 

As  the  dynamics  of  the  airfoil  motion  were  changed,  the  following 
impacts  on  the  dynamic  stall  process  were  observed: 


1  )  Pi tch  Rate 


ELOCITY 
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Fig.7c  Magnitude  of  the  Velocities  over  the  airfoil  surface,  in 
Percent  Freestream,  at  Sensor  Number  Three.  Pitch  Rate 
is  920°/s.  Pitch  Axis  is  a)  0.25c,  b)  0.5c,  c)  .7 5c. 


USAFA-TR-85-2 


It  is  clearly  evident  that  as  the  pitch  rate  is  increased, 
the  angle  of  attack  at  which  dynamic  stall  occurs  is  also  increased. 
This  effect  is  nonlinear  and  more  prevalent  at  increments  in  low  pitch 
rates  than  at  the  higher  rates.  In  fact,  although  using  a  different 
range  of  parameters,  Gormont  suggests  that  the  angle  of  attack  at  which 
dynamic  stall  occurs  varies  with  the  square  root  of  the  pitch  rate  (Ref 
8).  In  addition,  the  slope  of  the  velocity  profile  with  respect  to 
angle  of  attack,  close  to  the  surface  and  near  the  leading  edge,  does 
not  change  as  the  pitch  rate  is  increased.  Also,  the  corresponding 
peak  velocities  near  the  leading  edge  change  very  little.  There  are 
direct  indications  that  increased  pitch  rates  form  more  coherent  and 
energetic  dynamic  stall  vortices.  As  the  vortices  become  more 
energetic,  (as  was  shown  by  increased  velocity  magnitudes  in  the 
profiles  aft  of  the  leading  edge),  they  remain  over  the  airfoil  for  a 
longer  period  of  time.  This  would  dramatically  influence  the  time 
averaged  aerodynamic  forces  on  the  airfoil  at  high  angles  of  attack. 

2 )  Pitch  Axis 

Moving  the  pitch  axis  toward  the  trailing  edge  emulates  some 
of  the  effects  of  an  increased  pitch  rate.  Dynamic  stall  occurs  at 
correspondingly  higher  angles  of  attack,  and  the  large  dynamic  stall 
vortex  again  has  significant  impact  on  the  velocities  over  the  airfoil 
surface.  However,  differences  do  exist  in  the  velocity  profiles  and 
vortex  dynamics  from  those  noted  for  an  increased  pitch  rate.  The 
slope  of  the  velocity  magnitude  profiles  at  the  first  sensor  changed 
substantially;  the  velocity  increased  at  a  slower  rate  and  did  not 
reach  as  high  a  peak  when  the  pitch  axis  was  moved  toward  the  trailing 
edge.  Also,  the  dynamic  stall  vortex,  once  initiated,  evolved  rapidly 
and  moved  quickly  downstream  regardless  of  the  pitch  axis.  Substantial 
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similarities  were  seen  in  the  velocity  magnitude  profiles  of  the 
mid-chord  region  for  the  various  pitch  axis  locations.  The  profiles 
had  the  same  relative  peaks  and  magnitudes  with  the  only  difference 
being  a  slight  temporal  increment  as  the  pitch  axis  was  moved  aft. 

These  observed  phenomena  show  the  dramatic  effects  of  pitch  rate 
and  pitch  axis  on  the  dynamic  stall  process.  Clearly,  more 
quantitative  aerodynamic  measurements  are  needed  to  precisely  select  an 
optimal  parameter  combination  to  enhance  performace.  Experimental 
studies  to  quantify  these  and  other  effects  are  underway  using  more 
extensive  hot-wire  techniques  and  minature  surface  mounted  pressure 
t  r  ansducer s . 
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AN  EXAMINATION  OF  CURVE  SMOOTHING  USING  DIGITAL  FILTER  THEORY 
Paul  I.  King,*  and  Martin  L.G.  Oldfield** 

Abstract 

The  theory  of  digital  filters  is  used  to  explain  the  effect  of 
applying  a  smoothing  function  to  data  which  contains  unwanted  frequency 
components.  The  smoothing  of  a  curve  is  viewed  as  the  convolution  of 
the  raw  data  with  the  impulse  response  values  of  a  finite  impulse 
response,  non - r ecu r s i v e  filter.  By  the  convolution  theorem  the 
frequency  response  of  the  smoothed  curve  is  then  the  product  of  the 
response  of  the  raw  data  and  the  frequency  response  of  the  filter.  Two 
filter  impulse  responses  are  examined,  the  windowed  moving  average  and 
the  moving  least  square  polynomial.  The  windowed  moving  average  is  easy 
to  use  and  for  the  same  number  of  smoothing  coefficients  produces  a 
qualitatively  smoother  curve.  For  studying  transient  peak  data  the 
least  square  polynomial  is  the  best  choice. 


* •  1 n  t  roduc  t i on 

In  recent  years  the  development  of  high  speed  ana  1 og - t o-d i g i t a  I 
(A/D)  recording  devices  has  solved  the  problem  of  sampling  and  storing 
certain  measurements  of  physical  phenomena  at  appropriate  rates  for 
detailed  postmortem  analysis.  At  the  same  time  the  evolution  of  fast 
Fourier  transform  algorithms  along  with  the  development  of  digital 
filter  theory  has  given  engineers  and  scientists  powerful  mathematical 
and  computational  tools  for  the  analysis  of  complex  waveforms  inherent 
in  the  recorded  data. 

It  is  not  uncommon  nowadays  to  acquire  a  large  amount  of  data  in  a 
short  period  of  time  such  as  is  shown  in  Figure  1.  This  figure  shows 
the  results  of  pressure  measurements  made  in  a  turbulent  boundary 
layer.  The  turbulent  fluctuations  represent  a  numerically  noisy  signal 
and  it  is  difficult  to  extract  meaningful  information  from  this  raw 
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data.  A  common  technique  is  to  numerically  smooth  the  data  so  that 
qualitative  analysis  is  facilitated  and  numerical  differentiation  is 
made  possible.  What  follows  is  an  examination  of  smoothing  viewed  as  a 
process  of  filtering  raw  data  --  that  is,  a  process  in  which  one 
retains  (passes)  relevant  data  and  discards  (attenuates)  spurious  or 
irrelevant  data. 


RAW  BOUNDARV  LAYER  DATA 


Figure  1.  Raw  Data  from  High  Speed  Data  Acquisition 

The  filtering  (smoothing)  of  data  after  the  experiment  holds  many 
advantages  including  affording  one  the  time  to  examine  data  away  from 
the  hectic  environment  of  the  experiment  itself.  More  importantly  one 
can  structure  the  filter  network  (smoothing  function)  so  as  to  avoid 
the  problems  intrinsic  to  real-time  analog  filtering.  The  ability  to 
look  ahead  in  time  allows  the  digital  filter  to  avoid  the  time  delay 
and  phase  shift  common  to  analog  filters. 

This  paper  discusses  the  considerations  important  to  the  selection 
of  a  smoothing  function  for  a  particular  application.  In  doing  so  the 
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frequency  response  of  the  function  must  be  considered  and  thus  some 
familiarization  with  transform  anatysis  and  convolution  must  be  gained. 

In  particular  two  filter  functions  are  examined:  the  simple  running 
average  and  the  running  least  square  polynimial  fit. 

1  * •  The  Smoothing  Function  Viewed  as  a  Filter 

In  the  following  analyses,  the  data  are  assumed  to  be  a  function 
of  a  single  variable,  such  as  amplitude  vs  time  or  amplitude  vs 
distance.  For  discussion  purposes  it  is  assumed  that  amplitude  is  the 
ordinate  value,  and  time  or  distance  is  the  value  on  the  abscissa  axis. 

The  discussion  on  frequency  analysis  is  applicable  to  any  function  of 
a  single  variable. 

One  can  view  the  unsmoothed  data  as  a  synthesized  signal  composed 
of  oscillating  waves  of  varying  amplitude  and  phase,  superimposed  on  a 
time  varying  mean  signal  whose  form  and  magnitude  is  to  be  determined. 
The  synthesized  signal  is  sampled  in  uniform  time  or  spacial 
increments.  The  samples  form  a  sequence  x(nT)  or  simply  x(n)  where 

x(n)  =  x(t)  _oo<n<®  (I] 

t  =  nT 


y 


% 


and  T  is  the  sampling  interval  (Ref.  1.). 

The  problem  in  curve  smoothing  is  finding  an  appropriate  filter  to 
act  on  the  sampled  data  after  the  fact  and  produce  a  result  having  the 
desired  frequency  characteristics.  The  features  of  this  problem  are 
shown  in  block  diagram  form  in  Figure  2.  Here  one  assumes  linearity  as 
well  as  a  correspondence  between  the  input  and  output  signals.  This 
correspondence  is  defined  by  the  convolution  (smoothing  formula), 


m=<*> 

y(n)  =  y"  x(n-m)C 


(2) 
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x(n)  UNSMOOTHED  DATA  ^ 

FILTER 

^  y(n)  SMOOTHED  DATA 

W 

—  - 

INPUT--> — PROCESSING  -  -►--OUTPUT 


Figure  2.  The  Smoothing  Function  Viewed  as  a  Filter 

The  value  y(n)  is  the  smoothed  value  formed  by  the  convolution  of  the 

unsmoothed  values,  x(n-m),  with  the  convolution  variables  C  .  These 

m 

variables,  C  ,  are  explained  in  the  theory  of  digital  filtering  as  the 
impulse  response  of  a  linear  finite  impulse  response  (FIR)  filter  (Ref. 
1  )  . 

The  convolution  summation  shown  above  produces  a  moving  weighted 

average  of  the  input  data.  One  example  of  such  a  moving  weighted 

average  is  the  set  of  convolution  variables  C  =  constant  =  1/N,  where 

m 

N  represents  the  length  of  the  convolution  interval.  This  is  a  finite 
convolution  formed  by  a  truncated  smoothing  formula  composed  of  N 
terms,  N  usually  being  an  odd  integer.  The  weighted  value  y(n) 
replaces  the  original  data  value  x(n).  The  smoothed  values  are 


y(n)  =  1/N 


m  = 


(n-1 ) 
2 


x(n-m) 

(n-1 ) 

2 


(3) 


This  is  the  simple  running  average  of  which  more  will  be  said  later. 

As  shown  in  Equation  3,  y(n)  depends  only  on  values  of  the  input 
samples.  This  type  of  correspondence  defines  the  coefficients  of  what 
is  commonly  called  a  non-recursive  FIR  filter. 
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As  with  any  filtering  process,  analog  or  digital,  the 
characteristics  of  the  filter  must  be  analyzed.  All  filters  will 
distort  and/or  remove  information  from  the  input  signal.  It  is 
convenient  to  think  of  the  filter's  action  on  the  data  as  involving  a 
transfer  function  in  the  sense  that  information  is  processed  in  order 
to  produce  an  output.  The  filter  coefficients,  Cm,  are  the 
coefficients  of  the  transfer  function  and  are  thereby  related  to  the 
response  of  the  filter. 

III.  The  Fourier  Transform  and  Convolution 

To  examine  the  response  of  the  filter,  it  is  necessary  to  review 
certain  characteristics  of  discrete  data  as  observed  in  the  time 
(real)domain  and  the  frequency  domain  --  the  transformed  time  domain. 
The  two  regions  are  mapped  by  the  Fourier  transform.  A  continuous  time 
based  signal,  h(t),  is  transformed  to  the  frequency  domain  with 
amplitude  H ( f )  via  this  Fourier  transform, 


H(f) 


h(t)e"J2irft 


dt 


(4) 


H(f)  is  the  frequency  response  of  the  signal  h(t)  and  is  generally  a 
complex  value,  indicating  that  both  amplitude  and  phase  must  be 
determined  for  any  frequency.  If  the  continuous  signal  is  sampled  N 
times  with  sampling  interval  T,  the  tranform  has  discrete  values 
computed  as  (Ref.  3), 

H  (  n/NT  )  =  E1  h(kT)e"j2irnk/N  (5) 

K  =0 
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A  shorthand  notation  for  this  mapping  process  is 

h(kT)  <>  HCn/NT)  (6) 

The  Fourier  transform  is  used  to  study  the  effects  of  a  smoothing 
function  in  conjunction  with  a  powerful  analytical  tool,  the 
convolution  theorem. 

Convolution  of  two  continous  signals,  x(t)  and  h(t),  is  defined  as 

CD 

y  C  t )  =  x(t)h(t-T)dT  (7 

.OO 

where  y(t)  is  the  result  of  summing  the  product  of  signal  x(x)  with  the 
lagged  signal  h(t-x).  This  is  written  in  notational  form 

y(t)  =  x(t)  *  h(t)  (8) 

The  convolution  integral  is  best  understood  in  a  graphical  sense.  The 
signal  h(t)  is  folded,  or  rotated,  about  the  ordinate  axis  and 
displaced  an  amount  t.  The  overlap  region  of  x(t)  and  h(t-i)  is  then 
multiplied  point  for  point  and  the  products  summed  to  yield  y(t). 

Recall  that  this  is  similar  to  the  process  used  to  obtain  the  smoothed 
values  y(n)  in  Equation  2. 

If  one  is  processing  data  acquired  in  discrete  form,  the 
convolution  integral  (summation)  for  period  N  is  written 

N-l 

y(kt)  =  £  x(iT)h((k-i)T)  (9) 

i=0 

or  in  notational  form, 


y(kt)  =  x(kt)  *  h(kt ) 


(10) 
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A  mathematical  relationship,  the  convolution  theorem,  relates  the 
Fourier  transform  and  the  convolution  integral.  In  notational  form 
this  is  writ  ten 

X(t)  *  h(t)  <>  X(f)H(f) 


IV.  The  Smoothing  Function's  Relationship  to  Frequency  Response 


As  suggested  by  Equation  11,  the  Fourier  transform  of  the 
convolution  product  is  the  product  of  the  individual  transforms.  Thus 
convolution  in  the  time  domain  is  mathematically  equivalent  to 
multiplication  in  the  frequency  domain.  This  can  be  related  to  the 
f i  1  ter ing  process . 

Filtering  is  usually  represented  in  the  frequency  domain  as  the 
product  of  the  signal  response  with  the  filter  response  as  shown  in 
Figure  3.  From  the  convolution  thereom,  this  is  equivalent  in  the  time 
domain  to  convolution  of  the  signal  with  the  inverse  transform  of  the 
filter  frequency  response.  When  a  filter  is  subjected  to  an  impulse 
input,  the  Fourier  transform  of  its  time  response  to  the  impulse  is  the 
filter's  frequency  response.  Therefore,  the  inverse  transform  just 
referred  to  is  called  the  filter  impulse  response.  It  is  evident 
that  the  smoothing  process  can  be  analyzed  by  examining  either  the 
smoothing  function  time  characteristics,  or  its  frequency  response  in 
the  transformed  domain.  Furthermore,  if  one  can  specify  the  desired 
frequency  response  of  the  filter,  one  can,  in  theory,  invert  the 
response  to  obtain  the  time  domain  values,  C  .  These  ideas  form  the 
basis  of  digital  filtering  theory. 


Filtering  as  Viewed  in  the  Frequency  Domain 
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Stearns  (Ref.  4)  shows  that  the  frequency  response  of  a  numerical 

smoothing  function  is  the  discrete  Fourier  transform  of  some  function 

C(t)  which  has  been  sampled  and  stored  as  C(mT)  or  simply  C  .  He 

demonstrates  that  the  convolution  variables,  C_,  as  shown  earlier 

m 

(Equation  2)  do  indeed  represent  the  impulse  response  of  the  filter. 
Substituting  these  values  for  h(kT)  in  Equation  5,  the  discrete 
frequency  response  is 


N-l  -j2irnm/T 

C(n/NT)  =  £  Cme 

m=0 


(12 


Conversely,  if  the  frequency  characteristics,  C(n/NT),  are  known  or 
specified,  then  one  obtains  the  filter  synthesis  formula  (Ref.  3). 


N-l 

C  =  1/N  £  C(n/NT)ej2lrnm/T  (13 

n=0 

V .  Choosing  an  Appropriate  Filter  Function 

At  this  point,  it  is  instructive  to  examine  the  frequency 
response  of  some  candidate  smoothing  functions.  Figures  4a  and  5a 
represent  the  impulse  responses  of  two  common  filters,  the  moving 
average,  and  ideal  filter  respectively. 
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B.  IDEAL  FILTER 


Figure  5.  Time  Response  Required  for  Ideal  Filter 
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As  mentioned  earlier,  the  simplest  function  is  the  moving 
average,  which  in  the  time  domain  is  a  rectangular  truncation  function. 

The  function  and  its  transform  are  shown  in  Figure  4  (only  positive 
frequencies  are  shown).  The  sidelobes  in  the  frequency  domain  are  due 
to  truncation  in  the  time  domain  and  are  indicative  of  what  happens 
when  truncation  occurs  in  either  the  time  or  frequency  domain  (see 
Figure  5).  The  moving  average  thus  introduces  undesirable  frequency 
characteristics,  and  one  is  motivated  to  search  for  a  more  suitable 
funct ion. 

An  obvious  selection  would  seem  to  be  the  so  called  ideal  filter 
for  which  frequency  response  is  the  rectangular  truncation  function  as 
shown  in  Figure  5b.  The  impulse  response  of  the  ideal  filter  has 
sidelobes  which  extend  to  infinity  and  which  are  caused  by  the 
truncation  (sharp  cutoff)  in  the  frequency  function.  To  accurately 
represent  such  a  function  would  require  a  series  of  infinite  length. 

To  obtain  a  practical  convolution  series  one  usually  would  truncate  the 
infinite  series  after  n  terms  as  an  approximation  to  the  full  series. 
However,  truncation  after  the  nth  term  will  cause  the  well  known  Gibb's 
phenomenon,  an  overshoot  and  oscillatory  frequency  response  (Figure  6 
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In  summary,  there  are  two  problems  associated  with  choosing  a 
filter  function  which  will  yield  the  desired  results  and  which  is 
mathematically  realizable.  First,  it  is  always  necessary  for  any 
realizable  function  to  be  limited  in  time.  That  is,  it  must  be  of 
finite  length  (a  truncated  function);  the  values  Cm  must  equal  zero 
beyond  a  specified  upper  limit  on  m.  As  shown  above,  this  will  lead  to 
oscillatory  and  undesirable  frequency  characteristics  in  the  frequency 
domain.  Second,  it  is  necessary  for  the  sidelobes  of  the  impulse 
response  to  be  as  small  as  possible.  This  implies  that  the  frequency 
response  cannot  have  a  sharp  cutoff.  Figure  5  shows  that  a  sharp 
cutoff  necessitates  sidelobes  of  significant  magnitude  in  the  impulse 
response  coefficients.  These  two  problems  would  seem  to  preclude  the 
use  of  the  ideal  filter  function  or  simple  running  average  in  their 
usual  form. 

V I .  Windowing 

The  solution  to  these  problems  is  to  window  the  filter  function 
coefficients  (Ref.  1).  That  is,  weight  the  coefficients  in  some  manner 
to  reduce  the  end  effects  caused  by  truncation.  Two  approaches  are 
recommended.  The  first  (Ref.  4)  is  to  window  the  truncated  impulse 
response  series  shown  in  Figure  6.  This  series  represents  the 
approximate  impulse  response  of  the  ideal  bandpass  filter.  As  shown  in 
Figure  7,  by  using  an  appropriate  window  function  (discussed  below)  it 
is  possible  to  reduce  frequency  domain  ripple  (Gibb's  phenomeon). 
Windowing  also  diminishes  the  strength  of  the  sidelobes  which  results 
in  a  loss  of  sharpness  of  the  frequency  cutoff. 
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A  second  and  mathematically  simpler  method  is  to  window  the 
running  average  coefficients  C(T)=1/N.  The  same  benefits  accrue  as  in 
windowing  the  ideal  filter,  but  a  further  loss  in  cutoff  sharpness 
results.  However,  the  ease  of  the  method  outweighs  this  problem,  and 
for  the  purposes  of  discussion  this  method  will  be  the  one  referred  to 
for  the  remainder  of  this  paper. 

Many  window  functions  exist  (Ref.  1),  but  for  discussion  purposes 
only  the  two  most  often  used  for  data  smoothing  will  be  examined  here. 

The  first  is  the  cosine  bell  or  Hamming  window.  This  function  can  be 
expressed  in  discrete  form  as 


h(n)  =  a  +  ( 1-a)  cos2irn/N 

where  h(n)  represents  discrete  values  of  the  window  function.  The  best  a 
is  0.54  (Ref.  1). 

Using  this  a  and  multiplying  by  the  coefficients  of  the  running 
average  function,  one  obtains  the  coefficients  of  the  windowed  running 
average  function  as 


C  =  (1/.54NH.54  ♦  . 46cos2irm/N )  1 

m 

The  value  0.54  has  been  placed  in  the  denominator  for  the  purpose  of 
scaling  the  frequency  response  to  the  value  one  at  zero  frequency. 

Another  often  used  window  function  is  the  least  squares  polynomial 
developed  by  Savitzky  and  Golay  (Ref.  2).  This  method  involves  fitting 
an  mth  order  polynomial  through  the  data  using  the  least  squares 
criteria  for  each  interval  of  2n  +  l  points  in  the  smoothing  convolution. 
The  new  central  value  of  the  group  is  the  value  that  lies  on  the 
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fitted  curve  at  the  corresponding  data  location. 

Savitzky  and  Golay  (Ref.  2)  have  tabulated  the  values  of  the 
convoluting  variables,  C  ,  for  a  number  of  polynomials  and  their 
derivatives  beginning  with  a  quad r a t i c -cub i c  polynomial.  (The 
quadratic  and  cubic  polynomials  have  identical  convolution 
coefficients.)  This  polynomial  is  adequate  for  most  work  since  "any 
smooth  curve  will  look  more  or  less  like  a  quadratic  in  the  vicinity  of 
a  peak  or  like  a  cubic  in  the  vicinity  of  a  shoulder"  (Ref.  2).  The 
smoothing  coefficients  for  this  function  are 

(16) 

where  N  is  an  odd  number  representing  the  length  of  the  function  and  m 
varies  from  one  to  N.  One  advantage  of  the  least  squares  polynomial 
method  is  the  relative  simplicity  of  the  derivative  functions.  The 
smoothing  coefficients  for  the  first  derivative  quadratic  function  are 

Cm  =  m  -  1/2(N+1 )  (17) 

The  discrete  Hamming  and  quad r a t i c -cub i c  least  square  functions 
are  displayed  in  Figure  8.  The  two  values  of  N  shown,  65  and  129,  were 
chosen  for  mathematical  ease  in  computing  discrete  Fourier  transforms 
by  computer.  These  two  functions  possess  the  desired  impulse 
characteristics  of  a  mathematically  suitable  smoothing  function  -- 
finite  length  and  reduced  sidelobes.  The  discrete  Fourier  transform  of 
these  functions  are  shown  in  Figure  9.  As  indicated  in  the  figure,  the 
frequency  response  of  these  functions  are  also  finite  functions  with 
reduced  sidelobes.  Also  indicated  in  Figure  9,  for  two  functions  of 


170 


USAFA-TR-85-2 


the  same  length,  the  Hamming  function  has  a  lower  cutoff  but  attenuates 
the  higher  frequencies  more  than  the  least  squares  function. 


MAGNITUDE 
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The  results  of  using  65-point  functions  to  smooth  the  data  from 
Figure  1  are  shown  in  Figure  10,  in  which  the  raw  data  has  been 
omitted.  One  can  see  that  the  two  functions  produce  a  similar  mean 
curve;  however,  the  least  square  result  contains  more  ripples  while  the 
Hamming  window  result  is  attenuated  more  in  the  region  where  the  mean 
is  changing  rapidly.  In  Figure  9,  one  can  observe  that  the  65-point 
Hamming  result  is  similar  to  the  129-point  least  square  result.  This 
would  seem  to  indicate  that  an  N-point  Hamming  smoothing  is  similar  to 
a  2N-point  least  square  smoothing.  If  one  were  truncating  unsmoothed 
end  points,  the  Hamming  function  would  allow  for  a  larger  number  of 
retained  end  points.  This  is  demonstrated  in  Figure  11,  in  which 
smoothing  has  been  carried  out  on  the  data  of  Figure  1.  Using  the 
65-point  and  129-point  functions  for  each  smoothing  operation,  N-l 
points  have  been  truncated,  half  at  each  end.  (The  data  was  not 
smoothed  to  the  full  right  hand  end  in  order  to  conserve  computer 
time.)  The  results  of  the  65-point  Hamming  smoothing  are  qualitatively 
similar  to  the  129-point  least  square  smoothing. 
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129  POINT  HAMMING  FILTER 
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VI  I .Length  of  Smoothing  Function 

Once  one  has  decided  to  use  a  smoothing  function,  the  length 
(number  of  points)  must  be  determined.  This  will  be  governed  by  the 
type  of  data  one  is  examining  and  the  amount  of  data  that  one  can 
afford  to  lose  at  the  end  points.  Addressing  the  latter  point  first, 
it  is  recognized  that  for  a  2n+l  convolution  function  the  central  point 
(n+l  th  point)  is  replaced  by  the  convolution  summation,  and  at  either 
end  of  the  data  set,  n  points  cannot  be  smoothed.  For  example,  in 
using  a  nine  point  function,  eight  points  are  unsmoothed  for  each  pass 
of  the  smoothing  function,  four  at  each  end.  Many  theories  exist  on 
how  to  treat  these  end  points.  One  is  to  use  a  shorter  function  on  the 
end  points  than  that  used  on  the  central  portion  of  the  data.  A  second 
solution  is  simply  to  truncate  the  end  points  after  each  smoothing 
convolution.  There  is  no  general  answer  to  this  question;  it  must  be 
decided  based  upon  the  smoothing  effects  desired. 

The  length  of  the  smoothing  function  will  depend  on  whether  one 
wants  to  examine  the  mean  of  the  data,  the  peak  waveforms,  or  a 
combination  of  the  two.  A  filter  function  used  to  identify  the  mean 
would  reduce  the  standard  deviation  of  the  data  about  the  mean. 
According  to  a  study  by  Enke  and  Nieman  (Ref.  5)  this  filter  function 
should  be  as  long  as  possible  but  small  enough  so  that  not  more  than 
one  inflection  in  the  observed  data  is  included  in  any  convolution 
interval  of  2n+l  points. 

In  order  to  study  peaks  in  the  data,  the  smoothing  function  should 
cause  as  little  distortion  in  waveform  shape  as  possible  and  at  the 
same  time  increase  the  ratio  of  peak  values  to  signal  noise.  Enke  and 
Nieman  suggest  that  this  type  of  data  is  best  filtered  by  a  function  in 
which  width  at  half-peak  height  equals  that  of  the  peaks  in  the  data 
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and  in  which  the  function  length,  2n+l,  is  twice  that  of  the 
af foremen t i oned  width.  This  is  closely  approximated  by  a 
quad r a t i c -cub i c  function.  This  criteria  will  lead  to  shorter  functions 
than  those  used  to  simply  identify  mean  values  and  will  produce  a 
compromise  between  noise  reduction  and  signal  distortion.  To  choose  a 
function  with  which  to  smooth  data  in  which  both  mean  values  and  peaks 
are  of  interest  will  require  a  judgement  based  on  the  relative 
importance  of  the  signal  mix. 

VIII. Cone  1 u  s i on  s 

In  summary,  the  Hamming  function  would  seem  to  provide  for  the 
least  end  point  truncation  and  be  best  suited  for  smoothing  data  in 
which  the  mean  is  of  highest  interest.  The  least  squares  function  on 
the  other  hand  appears  to  have  a  sharper  frequency  rolloff  and  thus  is 
well  suited  to  data  in  which  peak  values  and  waveform  shape  are  of 
interest.  The  choice  made  to  smooth  the  boundary  layer  data  shown  in 
Figure  I  was  governed  by  a  consideration  of  the  loss  of  end  points  in 
the  data.  Thus  the  Hamming  function  was  chosen  to  provide  the  best 
smoothed  results  for  this  data. 

Regardless  of  the  choice  of  smoothing  functions,  it  is  important 
that  the  effects  of  the  smoothing  procedure  be  understood.  When  one 
views  the  procedure  as  one  of  filtering  the  data  with  a  low  pass 
filter,  it  is  evident  that  not  only  can  data  be  obscured  or  even  lost 
but  spurious  data  can  be  introduced  (the  Gibb's  phenomenon).  The  use 
of  a  windowed  function  offers  the  best  solution  to  the  smoothing 
problem.  This  paper  has  compared  two  popular  functions:  the  Hamming 
(or  windowed  retangular  function)  and  the  moving  least  square 
polynomial . 
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FLIGHT  TEST  MEASUREMENT  OF  THE  AERODYNAMIC  EFFECT  OF  THRUST  LEVEL  ON 

LIFT  AND  DRAG 
T.R.  Yechout* 

Abstract 

A  flight  test  technique  has  been  developed  under  NASA  Dryden 
sponsorship  to  define  the  aerodynamic  effect  of  thrust  level  on 
aircraft  lift  and  drag  characteristics.  Conventional  stabilized  "spee 
power"  tests  require  that  the  thrust  be  adjusted  for  each  test 
condition  and,  as  a  result,  the  effect  of  thrust  on  aerodynamic 
characteristics  cannot  be  easily  identified.  The  technique  utilizes 
quasi  steady-state  maneuvers  at  selected  power  settings  throughout  the 
Mach  range  of  the  aircraft  to  define  liic  and  drag  coefficient 
variation  as  a  function  of  angle  of  attack,  Mach  number,  and  power 
setting.  A  twenty  hour  verification  flight  test  program  was 
accomplished  using  a  Learjet  Model  35  aircraft.  Significant  power 
effects  were  identified  which  should  be  anticipated  on  any  aircraft 
with  jet  engines  mounted  on  the  aft  fuselage  above  the  inboard  wing 
section. 

I .  In t r oduc t i on 


In-flight  definition  of  the  aerodynamic  effect  of  thrust  level  on 
lift  and  drag  characteristics  has  not  been  accomplished  in  the  past. 
Normally,  lift  and  drag  measurements  are  made  using  a  series  of 
stabilized  points  throughout  the  aircraft  flight  envelope.  A  wide 
range  of  engine  power  settings  are  used  to  achieve  the  stabilized 
conditions  from  which  lift  and  drag  may  be  determined  given  an 
in-flight  thrust  and  airflow  model  along  with  normally  instrumented 
aircraft  parameters  such  as  weight  and  angle  of  attack.  Unfortunately 
the  flow  field  around  the  aircraft  may  be  significantly  altered  by  the 
airflow  through  engine(s)  which  will  result  in  the  lift  and  drag 
characteristics  being  directly  dependent  on  engine  power.  If  the 
stabilized  point  method  is  used  on  an  aircraft  with  significant  power 
effects,  use  of  the  resulting  data  to  predict  non- s t ab i 1 i zed 
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(i.e.  excess  thrust  not  equal  to  zero)  performance  characteristics  will 
be  susceptible  to  significant  error.  A  flight  test  technique  has  been 
developed  to  efficiently  evaluate  the  effect  of  engine  power  setting  on 
the  lift  and  drag  characteristics  of  an  aircraft.  The  technique 
utilized  quasi  steady-state  maneuvers  (level  accelerations  and 
decelerations)  at  selected  power  settings  throughout  the  Mach  range  of 
the  aircraft  to  define  lift  and  drag  coefficient  variation  as  a 
function  of  angle  of  attack,  Mach  number,  and  power  setting.  The 
technique  was  developed  and  evaluated  during  a  twenty  hour  flight  test 
program  using  a  Learjet  Model  35  aircraft.  This  effort  was  part  of  an 
overall  research  program  which  concentrated  on  modeling  aircraft 
performance  throughout  the  flight  envelope.  The  use  of  quasi 
steady-state  maneuvers  not  only  allowed  definition  of  power  effects  on 
lift  and  drag,  but  also  provided  a  very  time  efficient  approach  to 
overall  in-flight  aircraft  performance  definition  when  compared  to  the 
stabilized  point  method. 

I  I .  Concept 

Development  of  the  lift  and  drag  characteristics  from  quasi 
steady-state  maneuvers  began  with  consideration  of  the  forces  acting  on 
the  aircraft.  The  aircraft  force  balance  equations  resolved  parallel 
and  perpendicular  to  the  flight  path  (assuming  zero  sideslip,  wings 
level,  and  constant  mass)  are,  from  Figure  I, 


F 


x 


ma. 


U 


F  cos  (a+A)  -  F 


D  =  W  (r—  +  siny) 


(S 
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L  +  F  sin  (a+X)  =  +  cosy) 

g  9 


L 


Figure  1.  Aircraft  Force  Balance  Diagram 

As  discussed  in  Reference  l,  the  flight  path  load  factors  resolved 
along  the  x  and  z  wind  axes  are 

nx  =ir  +  siny 

nz  =-r +  C0SY 

The  force  balance  equations  may  then  be  expressed 
FgCOS  (a+X)  -  Fr  -  D  =  Wnx 


L  +  F_sin  (a+X)  =  Wn 
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and  the  relationships  for  lift  and  drag  coefficient  are 

Wn  -F  sin  (a+A) 

C  =  — - — 3 - 

L  ( 1/2 J  Y* p  M2S 

a 


(9) 


Fg  cos  (a+A)  -  Ff  -  Wnx 

( 1/2 )  Y*  p  H  2S 


(10] 


Equations  9  and  10  express  lift  and  drag  coefficients  in  terms  of  wind 
axes  accelerations,  the  engine  model  parameters  of  gross  thrust  and  ram 
drag,  and  normally  recorded  flight  test  parameters  such  as  angle  of 
attack,  ambient  pressure,  weight,  and  Mach  number.  The  equations  are 
compatible  with  quasi  steady-state  maneuvers  in  which  excess  thrust  is 
not  equal  to  zero  and  flight  path  accelerations  are  present.  Wind  axis 
accelerations  were  determined  from  accelerometers  mounted  along  the 
body  axis  of  the  aircraft  using  the  appropriate  angular 
transformations.  A  unique  inflight  thrust  and  airflow  prediction 
technique,  termed  "Thrust  Modeling,"  was  developed  as  part  of  the 
overall  program  to  define  the  gross  thrust  and  ram  drag  terms  in 
Equations  9  and  10.  The  technique  consisted  of  correcting  the  engine 
deck  predictions  of  thrust  and  airflow  to  match  the  performance  of  the 
actual  engines  installed  in  the  aircraft  using  a  three  step  approach: 

1)  Simplified  representation  of  engine  deck  predicted  thrust, 
fuel  flow,  and  airflow  in  corrected  form, 

2)  Correction  of  the  engine  deck  model,  developed  in  step  one,  to 
the  individual  characteristics  of  each  engine  based  on  a 
static  thrus  t  run , 

3)  In-flight  correction  of  thrust  and  airflow  predictions  based 
on  actual  test  fuel  flow,  an  accurate  specific  fuel 

consumption  prediction  and  a  balance  of  the  thrust  momentum 
equation. 
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The  final  equations  in  simplified  form  for  thrust  and  airflow  were: 


W  F 

Ftest  ^deck 


9test 


nW< 


(11) 


deck 


=  w, 


test 


test 


‘deck 


nWi 


- 1 


deck 


(12) 


where  n  was  the  ratio  of  static  thrust  run  specific  fuel  consumption  to 

engine  deck  specific  fuel  consumption  defined  as  a  function  of 

corrected  RPM.  A  complete  development  of  this  technique  is  presented 

in  Reference  2.  Thrust  and  air  flow  prediction  accuracies  were 

believed  to  be  three  to  five  percent  or  better  based  on  data  obtained 

from  the  Lear  55,  F-104G,  F - 1 1 1 ,  and  YF-12  programs  (Refs.  2,  3,  4,  and 

5).  This  prediction  technique  offered  several  advantages  over  the  most 

commonly  used  methods  as  discussed  in  Reference  2. 

Two  corrections  were  made  to  lift  coefficient  for  elevator  trim 

effects  resulting  from  1)  the  thrust  moment  about  the  center  of  gravity 

(c.g.)  and  2)  a  non-standard  c.g.  These  corrections  standardized  the 

lift  coefficient  data  to  a  common  baseline. 

In  the  first  case,  the  effect  of  the  associated  moments  created  by 

the  thrust  (F  )  and  ram  drag  (F  )  about  the  c.g.  were  removed.  From 
g  r  6 

Figure  2,  this  moment  is  given  by 


(13) 
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The  distance  Ztfirust  is  a  function  of  c.g.  and  airframe  geometry,  while 
hf  and  £  ^  are  also  functions  of  angle  of  attack. 


Figure  3.  Moment  Arms  for  C.G.  Standardization 

C.  was  also  standardized  to  a  particular  c.g.  location.  From 
lT 

Figure  3,  this  correction  begins  with  a  moment  balance 

ALwing  =  Ltai 1  'tail  (test  c.g.)  (17; 

Lwing  (A  “  Ac-9-  )  =  z(Ltail  +  ALtail')  (18) 

(standard  c.g.) 

where  Ltai|  is  the  change  in  tail  lift  required  for  a  standard  c.g. 
Since 

1  =  'tail  +  Ac.g. , 

Ltai  1  'tail  '  Lw ing  Ac-9-  = 

(Ptail  +  Ltail  +  ^ *tail  +  Ac.g. )ALtall.  (19) 


With  the  total  aircraft  lift  (L)  given  by 
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L  Lwing  +  Ltail 

for  the  test  condition,  Equation  19  becomes 


(20 


-L(Ac.g. )  =  (^tail  +  Ac‘g*  •)ALtail, 


(2 


and 

-L ( Ac .g .  ) 

^Ltail  “  *tail  +  Ac.g.  (2 


I n  coe  f  f i c i en  t  form 


AC, 


'c.g. 


-C.  (Ac.g.) 

..hi _ 

(etau+  ^-g*3 


(2 


and  the  standardized  lift  coefficient  corrected  for  thrust  moment 
effects  and  to  a  standard  C.g.  is 


C 


L 


S 


+  AC, 


"c.g, 


(2 


A  correction  to  the  drag  coefficient  was  made  for  skin  friction 
variation  as  a  function  of  Reynolds  number  to  standardize  the  drag 
coefficient  data  to  a  particular  altitude.  Sch 1 icht ing' s  formula  for 
the  skin  friction  coefficient  assuming  turbulent  flow  (Ref.  6)  was 
used . 


_ .455 _ 

( J.°g  LnRe  ) 2  *  58  ( 1  +  .144M2)-65 


R  =  —  f  5  characteristic  length 

e  p  p  s  viscosity  coefficient 
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The  drag  coefficient  due  to  skin  friction  is  then 

f Wetted  Area  \ 

-  cf  V  i  J 


SF 


(27; 


where  the  drag  on  the  aircraft  is  obtained  by  adding  the  drag  on  each 
of  the  components  shown  in  Table  1.  The  Reynolds  number  calculation  of 
Equation  26  requires  the  charac ter i s t ic  length,  and  the  drag 
coefficient  calculation  of  Equation  27  requires  the  applicable  wetted 
area.  The  skin  friction  drag  contributions  were  then  standardized  to 


an  altitude  of  25000  feet  by  computing  Cf 


and  Cr 


and 


25000' 


JSF. 


25000' 

defining  the  incremental  change  in  drag  coefficient  due  to  skin 


friction  variation  for  off-standardized  conditions  as 


AC 


0SF  =  CDc,e.  -  C0 


SF  SF 

25000 


(28 


This  methodology  was  used  for  each  of  the  aircraft  components;  and  the 

total  skin  friction  drag  correction,  Cn  ,  was  obtained  by  summing 

°SF 

the  contribution  of  each  component.  total 


AC 


SF 

total 


ACn  +  AC 
USF 

fuselage 


SF 

wing 


+AC 


D 


SF 

h.tall 


♦  AC 


+  AC 
tall 
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+  AC 
SF 

pylon 


0 


SF 

nacelles 


+  ACn  +  AC 
°SF 

ventral 

fin 


SF 

tank 


+  AC 


D 


SF 

tank  fin 
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The  standardized  drag  coefficient,  , 


is  then 


C 


+  AC, 


SF 

total 


(30) 


An  altitude  of  25000  feet  was  chosen  for  standardization  since  it  was 
approximately  in  the  middle  of  the  altitude  envelope  of  the  aircraft. 


Fuselage 

Wing 

Hor izontal  ta i 1 
Vertical  tail 
Engine  pylons 
Engine  nacel les 
Ventral  fin 
Tip  tanks 
Tank  fins 


Table  1:  Aircraft  Components  Used  to  Calculate  Skin  Friction  Drag 


and  Cp  versus  angle  of  attack  characteristics  were  defined  from  a 
series  of  test  points  obtained  during  acceleration  and  deceleration 
maneuvers.  These  characteristics  were  defined  as  a  function  of  power 
setting  and  Mach  number.  The  needed  lift  coefficient  range  was 
obtained  through  variation  of  the  we igh t- pr es su r e  ratio  (W/tf)  as 
discussed  in  the  Test  Procedure  section.  By  determining  the  lift  and 
drag  characteristics  as  a  function  of  power  setting,  the 
power -dependen t  effects  could  be  defined  when  comparing  data  for  the 
same  Mach  number  and  angle  of  attack. 

III.  Test  Procedure 

Quasi  steady-state  acceleration/deceleration  maneuvers  provided 
the  necessary  data  to  define  lift  and  drag  characteristics.  These 
maneuvers  were  conducted  at  nearly  constant  altitude  using  the  altitude 
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hold  mode  of  the  autopilot.  Normally  less  than  a  60  foot  excursion 
from  the  start  altitude  was  experienced  during  a  maneuver.  Eight 
"cardinal"  power  settings  were  evaluated  consisting  of  95,  90,  85,  80, 
75,  70,  60,  and  50  percent  low  pressure  fan  RPM  (Nj).  The  Nj  was 
chosen  as  the  variable  to  represent  power  because  of  the  relatively 
high  bypass  ratio  of  the  engines  and  the  resulting  high  correlation  to 
engine  airflow.  An  acceleration/deceleration  was  conducted  at  a 
cardinal  power  setting  by  holding  Nj  to  within  1/2  percent  during  a 
maneuver.  A  range  of  the  weight-pressure  ratio  parameter  (W/<5)  within 
the  aircraft  envelope  was  designated  to  provide  a  lift  coefficient 
variation  for  a  given  Mach  number  so  that  Mach  effects  could  be 
defined.  Eight  values  of  W/g  were  evaluated  as  shown  in  Table  2. 

These  eight  values  of  W/s  provided  eight  evenly  spaced  points  on  a 
constant  Mach  drag  polar  in  the  mid-Mach  range.  At  each  value  of  W/ 6 , 
an  acceleration/deceleration  sequence  was  performed,  which  included 
maneuvers  at  all  cardinal  power  settings  above  idle.  As  W/ s  increased, 
the  number  of  available  power  settings  decreased  because  the  idle  RPM 
increases  with  altitude.  For  example,  at  40,000  feet  only  the  95,  90, 
and  85  percent  power  settings  could  be  evaluated.  As  a  result,  the 
largest  amount  of  data  was  obtained  for  the  higher  power  settings. 


Nominal 

w/s 

Altitude 

22,000 

10.000 

ft 

40,000 

23,000 

ft 

47,000 

26,000 

ft 

53,000 

29,000 

ft 

60.000 

32,000 

ft 

67,000 

35,000 

ft 

73,000 

38,000 

ft 

80,000 

40,000 

ft 

Table  2:  Performance  Modeling  Maneuvering  Sequences 


Figure  4. 

Typical  Maneuvering  Sequence 

Sequence 

Power 

Settina 

Maneuver 

Data 

Recorded 

1 

95 

Accel 

Yes 

2 

90 

Oecel 

Yes 

3 

70 

Oecel 

Yes 

4 

90 

Accel 

Yes 

5 

95 

Accel 

No 

6 

85 

Oecel 

Yes 

7 

70 

Oecel 

Yes 

8 

85 

Accel 

Yes 

9 

95 

Accel 

No 

10 

80 

Oecel 

Yes 

11 

70 

Oecel 

Yes 

12 

80 

Accel 

Yes 

13 

95 

Accel 

No 

14 

75 

Oecel 

Yes 

15 

70 

Oecel 

Yes 

16 

75 

Accel 

Yes 

1S0 
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A  typical  maneuvering  sequence,  illustrated  in  Figure  4,  assumes 
that  the  drag  curve  and  engine  idle  level  are  as  shown  for  a  particular 
W/ 6  configuration.  A  sequence  began  by  slowing  the  aircraft  to  an 
acceptable  minimum  speed  (for  the  Lear  35  this  was  an  airspeed  slightly 
above  stick  shaker  speed)  at  an  altitude  based  on  the  target  value  of 
W/6.  A  95  percent  acceleration  was  then  performed.  When  the 
acceleration  had  slowed  to  approximately  a  quarter  knot  per  second,  the 
throttles  were  retarded  to  90%  and  a  deceleration  was  performed  until  a 
stabilized  point  was  approached.  The  sequence  then  continued  as  shown 
in  Figure  4.  Altitude  adjustments  were  made  at  convenient  times  in  the 
sequence  to  maintain  W/ 6  within  approximately  1  percent  as  weight 
decreased.  Although  not  specifically  shown  in  the  Figure  4  diagram,  a 
high  power  setting  for  sequence  9  was  used  to  accelerate  past  the  last 
stabilized  condition  so  that  the  deceleration  as  shown  in  sequence  10 
could  be  obtained.  Although  a  stabilized  "speed  power”  point  was 
generally  not  obtained,  the  Mach  number  for  which  the  drag  and  net 
thrust  curves  intersected  could  easily  be  estimated  based  on  the  Mach 
region  in  which  the  acceleration  and  deceleration  for  a  particular 
power  setting  were  terminated.  The  general  guideline  used  was  to 
accelerate  far  enough  past  the  last  stabilized  condition  so  that  the 
engine  RPM  would  achieve  stabilization  on  the  subsequent  deceleration 
before  reaching  the  Mach  number  of  the  last  stabilized  point.  Data 
were  taken  periodically  throughout  an  acce 1 er a t i on /dece 1 er at i on  rather 
than  continually  to  keep  the  volume  of  data  to  a  manageable  level. 
Ideally,  approximately  a  twenty-second  burst  of  data  was  recorded  as 
the  aircraft  passed  through  each  .05  Mach  increment.  The  actual  test 
sequence  performed  at  each  W/ 5  condition  depended  directly  on  the 
location  of  the  drag  curve  with  respect  to  the  net  thrust  levels.  For 
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example,  if  two  cardinal  power  settings  were  located  between  engine 
idle  and  the  bottom  of  the  drag  curve,  then  at  least  one  deceleration 
would  be  performed  at  each  of  these  power  settings.  The  maneuver 
sequence  was  designed  to  acquire  the  needed  data  in  a  time-efficient 
manner  and  also  to  be  easily  accomplished  by  the  flight  crew.  It 
clearly  met  these  objectives.  For  planning  purposes,  approximately  45 
minutes  were  required  to  accomplish  a  maneuvering  sequence  at  one  value 
of  W/<S  for  this  aircraft. 

IV .  Results 

The  versus  angle  of  attack  characteristics  for  the  Lear  35 

s 

fell  into  two  distinct  categories.  Above  .65  Mach,  power  effects  were 

negligible  but  distinct  Mach  effects  were  identified.  A  summary  of  the 

standardized  lift  coefficient  characteristics  in  this  high  Mach  region 

is  presented  in  Figure  5,  in  which  an  increase  in  Mach  number  resulted 

in  an  increase  in  as  well  as  the  slope  .  The  extrapolated 

portions  of  each  curve  are  identified  by  the  uniform  dashed  lines  as 

indicated.  Below  .65,  Mach  effects  were  negligible,  but  power  effects 

were  found  as  presented  in  Figure  6.  At  power  settings  above  60%,  a 

small  but  significant  increase  in  was  observed.  At  70%  power,  an 

s 

approximate  .01  increase  in  C,  resulted  throughout  the  angle  of  attack 

s 

range  when  compared  to  the  data  from  the  60%  power  curve.  As  power  was 

increased  to  75%  and  above,  an  additional  increase  of  approximately  .01 

over  the  70%  curve  was  found.  The  data  scatter  experienced  when 

defining  the  C.  versus  a  curves  presented  in  Figures  5  and  6  was  +  .02 
s 

maximum  (based  on  C.  )  with  approximately  95  percent  of  the  data 

Ls 

falling  within  +.01  of  the  defined  curves. 
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Figure  6.  Lift  Coefficient  Characteristics 

M  <  .65 
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The  power  effects  on  are  thought  to  be  directly  related  to  the 

close  proximity  of  the  engine  inlets  above  the  inboard  upper  wing 
surface  as  profiled  in  Figure  7.  Either  of  twp  effects  could  be 
present.  First,  the  flow  field  around  the  wing/nacelle  may  be  fairly 
normal  at  high  engine  speed;  however,  at  low  engine  speed,  inlet 
spillage  reduces  the  lift  over  the  inboard  section  of  the  wing  by 
retarding  the  flow.  Second,  above  60%  power,  the  increased  airflow 
through  the  engine  may  alter  the  flow  field  in  the  engine  nacelle/wing 
root  area  so  that  the  overall  circulation  around  the  inboard  wing 
section  is  increased,  resulting  in  a  corresponding  increase  in  lift. 
This  increase  in  lift  does  not  continue  with  increasing  power  settings 
above  75%  but  rather  remains  constant  at  approximately  the  75%  value. 
The  increased  airflow  through  the  engine  with  increasing  power  may 
produce  an  increase  in  lift  on  the  forward  portion  of  the  inner  wing 
but  may  also  result  in  flow  starvation  and  separation  near  the  trailing 
edge,  producing  an  offsetting  effect.  Obviously,  a  flow  field  survey 
in  the  engine  nacelle/wing  root  area  is  needed  to  help  explain  these 
power  effects.  Flow  tufting  in  this  area  would  be  an  excellent  first 
step  in  understanding  the  power  effects  observed  and  could  easily  be 
accomplished  concurrently  with  the  quasi  steady-state  maneuvers  used 


Figure  7.  Engine/Airframe  Configuration 
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As  with  lift  coefficient,  Cp  versus  angle-of  attack 

characteristics  fell  into  two  distinct  and  consistent  categories.  For 

.6  Mach  and  above,  power  effects  were  not  observed  but  Mach  effects 

were  identified.  A  summary  of  the  standardized  drag  coefficient 

characteristics  in  the  high  Mach  region  is  presented  in  Figure  8  where 

an  increase  in  Mach  number  generally  resulted  in  an  increase  in  Cp  for 

s 

a  given  angle  of  attack.  As  shown  in  Figure  8,  the  largest  increase  in 

Cp  with  Mach  number  was  projected  above  four  degrees  angle-of  attack, 
s 

For  .55  Mach  and  below,  Mach  effects  were  not  significant  but  power 

effects  were  found.  As  presented  in  Figures  9  and  10,  Cp  generally 

decreased  as  power  decreased  with  approximately  a  45  drag  count  band 

between  90%  and  50%  power  in  the  mid  angle-of -at  tack  region.  The  95% 

power  curve  intersected  and  crossed  over  the  90%  curve  at  two  locations 

and  dropped  below  the  90%  curve  in  the  mid  angle-of-at tack  region  as 

shown.  The  data  scatter  experienced  when  defining  the  Cp  versus  a 

curves  presented  in  Figures  8,  9,  and  10  was ±  .003  maximum  (based  on 

Cp  )  with  approximately  95  percent  of  the  data  falling  within  .001  cf 
s 

the  defined  curves. 
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As  with  the  lift  coefficient  curves,  the  complex  flow  interaction 

in  the  nacelle/wing  root  area  must  be  analyzed  to  understand  these 

characteristics.  Normally  it  would  be  expected  that  lower  drag  would 

occur  at  higher  power  settings  due  to  reduced  inlet  spillage.  This 

trend  is  seen  in  the  mid  ang 1 e-of -at  tack  region  for  90%  and  95%  power 

(Figure  9).  However,  this  is  obviously  not  the  only  factor  affecting 

the  drag.  Another  possible  interaction  may  be  an  increased  pressure  on 

the  aft  facing  wing  and  fuselage  surfaces  (a  drag  reduction)  resulting 

from  increased  inlet  spillage  at  the  lower  power  settings.  The  close 

proximity  of  the  engine  nacelle  to  these  surfaces  makes  this 

possibility  quite  feasible.  In  Figure  9,  the  crossover  experienced  by 

the  95%  curve  in  the  higher  ang 1 e-of -at t ack  region  (lower  speed) 

indicates  that  the  increased  pressure  phenomenon  may  become  predominant 

as  the  projected  area  of  aft  facing  wing  and  fuselage  surfaces 

increases  with  ang 1 e-of -at  tack  and  as  propagation  of  the  inlet  spillage 

air  also  increases  with  lower  speed.  Another  contributing  factor  may 

be  increasing  flow  starvation  and  separation  near  the  trailing  edge  of 

the  wing  with  increasing  power,  as  discussed  in  the  lift  section.  This 

would  account  for  the  increased  drag  observed  with  increased  power. 

The  absence  of  power  effects  on  drag  at  .6  Mach  and  above  is  probably 

due  to  the  low  propagation  of  inlet  spillage  air  at  higher  speeds. 

Again,  a  flow  field  survey  in  the  wing  root/nacelle  area  would  help 

clarify  the  causes  of  the  identified  drag  characteristics. 

A  selected  number  of  steady-state  (speed  power)  points  were  flown  for 

comparison.  Excellent  agreement  was  obtained  between  the  steady-state 

points  and  the  curves  presented  in  Figures  5  through  10.  The  maximum 

deviation  for  C,  was  .01  and  for  Cn  was  .001  which  was  clearly  within 
Ls  us 

the  observed  data  scatter.  The  maximum  magnitudes  of  the  thrust 
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moment,  c.g.  and  skin  friction  corrections  were  also  determined  for  the 

flight  test  conditions  evaluated.  The  maximum  values  for 

AC.  , AC.  AndACn  were  .003,  .005,  and  .0012 

Shrust  c.g.  USF 

mome n  t 

respectively,  which  indicate  that  these  corrections  had  a  relatively 
small  but  still  significant  effect  on  the  data. 

V .  Cone  1  us i on  s 

An  in-flight  technique  was  developed  to  define  the  aerodynamic 
effect  of  thrust  level  (power  effects)  on  aircraft  lift  and  drag 
characteristics.  A  time  efficient  test  approach  was  used  which  relied 
primarily  on  level  acceleration  and  deceleration  maneuvers.  Flight 
test  results  for  the  Learjet  Model  35  aircraft  showed  that  power 
effects  were  very  significant  and  must  by  considered  if  the  data  is  to 
be  used  for  predictions  of  aircraft  performance  under  non-steady  flight 
conditions.  The  power  effects  found  on  the  Lear  35  were  believed  to  be 
directly  related  to  the  complex  flow  field  in  the  wing  root  area 
resulting  from  the  overwing  mounting  of  the  engines.  The  approach 
developed  in  this  program  should  generally  be  applicable  to  a  large 
range  of  both  jet  and  propeller  powered  aircraft. 

Nomenclature 

ax  acceleration  along  x  wind  axis 

az  acceleration  along  z  wind  axis 

A  moment  arm  about  test  c.g.  for  L„. 

°  wing 

Cq  drag  coefficient 

Cp  drag  coefficient  corrected  for  skin  friction  variation  (Eq.  7) 

s 

lift  coe  f  f i c i en  t 

C,  lift  coefficient  corrected  for  thrust  moment  effect  and 

S 

non-standard  c.g.  (Eq.  4) 
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lift  coefficient  corrected  for  thrust  moment  effect 
skin  friction  coefficient 
gross  thrust 
ram  drag 

forces  along  x  wind  axis 
forces  along  z  wind  axis 
acceleration  of  gravity 
moment  arm  for  F 


A  tail 


'thrust 


moment  arm  about  test  c.g.  for 

Mach  number 

aircraft  mass 

low  pressure  fan  RPM 

x  wind  axis  load  factor 

z  wind  axis  load  factor 

ambient  pressure 

dynamic  pressure 

Reynolds  number 

wing  reference  area 

aircraft  weight 

air  flow 

fuel  flow 

longitudinal  axis 

normal  axis 

moment  arm  about  standard  c.g.  for  L t a j  j 

moment  arm  for  F^ 
angle  of  attack 
f 1 ight  path  angle 


■r 's' \ 


<“LV 


V"  ,W.  W.N-V-, 
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y'  ratio  of  specific  heat  for  air 

X  thrust  inclination  angle 

6  pressure  rat i o 

p  air  dens i t  y 
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FLIGHT  CONTROL  SYSTEM  DESIGN  FOR  A  COMPUTER  CONTROLLED  AIRCRAFT  WITH 

LIMITED  SENSORS 
Thomas  P.  Webb* 

Abstract 

A  complete  flight  control  system  for  a  small  computer  controlled 
aircraft  was  designed  using  only  yaw  rate,  heading,  lateral  load 
factor,  airspeed,  altitude,  and  rate  of  climb  feedback.  This 
multi-input  multi-output  control  problem  was  done  using  the  classical 
root  locus  technique  on  a  linearized  system  model.  The  performance  of 
the  flight  control  system  was  then  checked  using  a  12  degree-of - f reedom 
nonlinear  simulation.  The  simulation  results  revealed  surprisingly  good 
performance,  considering  the  limitation  on  sensors. 


1 .  1 n  t  roduc  t i on 

The  Department  of  Electrical  Engineering  at  the  United  States  Air 
Force  Academy  is  attempting,  through  one  of  its  senior  design  courses, 
to  design,  build,  and  fly  a  computer  controlled  aircraft.  The 
Department  of  Aeronautics  was  asked  to  help  design  the  flight  control 
system  to  be  implemented  by  the  on-board  digital  computer.  The  project 
involved  building  and  testing  a  wind  tunnel  model  of  the  aircraft  to 
determine  its  aerodynamic  characteristics,  performing  mass  tests  on  the 
actual  aircraft  to  determine  inertia  characteristics,  developing  a  12 
degree-of - f reedom  nonlinear  aircraft  simulation  program,  and  designing 
the  actual  flight  control  system.  This  report  describes  only  the  last 
task . 

I  I .  Aircraft  Description 

The  aircraft  acquired  by  the  Electrical  Engineering  Department  is 
an  off  the  shelf  hobby  radio  controlled  airplane  called  the  "Big  Stick" 
sold  by  Hobby  Shack  in  kit  form.  This  particular  aircraft  was  chosen 
for  its  large  size  and  docile  handling  qualities.  The  aircraft  is 


♦Major,  USAF,  Assistant  Prof.,  Dept,  of  Aeronautics,  USAFA. 
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configured  for  normal  radio  controlled  operation  to  allow  for  initial 
testing,  manual  backup/over ide  for  safety,  and  manual  takeoffs  and 
landings.  The  aircraft  is  propeller  powered  by  a  2.5  brake  horsepower 
(3HP)  two- s t r oke -eye  1 e  gasoline  Quadra  35  engine.  The  aircraft  (see 
Fig.  1)  has  a  wingspan  of  8.73  feet.  The  estimated  weight  with  full 
fuel  and  computer  on  board  is  30  pounds.  The  tricycle  landing  gear 
configuration,  as  shown  in  the  picture,  was  later  modified  to 
conventional  (tail  wheel)  for  structural  reasons  and  to  facilitate 
operation  on  grass. 


Figure  1.  The  Big  Stick  Airplane 

The  aircraft  is  controlled  by  conventional  ailerons,  rudder,  elevator, 
and  throttle.  Drawings  for  the  .122  scale  wind  tunnel  model  are 
contained  in  Appendix  A. 
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III.  The  Control  Problem 

The  purpose  of  the  flight  control  system  is  to  make  the  aircraft 
fly  an  arbitrarily  specified  (and  not  necessarily  straight)  path  given 
continuous  information  on  the  current  state  of  the  aircraft  through  a 
limited  number  of  sensors.  The  design  of  the  flight  control  system  was 
formulated  as  a  multi-input  multi-output  feedback  control  problem.  The 
actual  parameters  to  be  controlled  were  specified  as  altitude  (h), 
airspeed  (V),  and  heading  (  <|j)  . 


Figure  2.  The  Control  Problem 

Referring  to  Figure  2,  the  control  problem  can  be  visualized  as 
one  of  driving  the  values  of  h,  V,  and  <l>  to  those  of  hc,  V  ,  and  1J1  , 

respectively,  where  c  denotes  the  commanded  value. 

The  flight  control  system  makes  inputs  to  the  aircraft  by 
adjusting  the  settings  or  deflections  of  the  elevator  (t$e),  throttle 
(6.J.),  ailerons  (6a),  and  rudder  (<5r>.  These  control  settings  depend  on 
the  sensor  measurements,  which  contain  information  about  the  actual 
state  of  the  aircraft,  and  the  commanded  values  of  altitude,  airspeed, 
and  heading.  This  process  is  depicted  in  Figure  3. 


P{ 


FLIGHT  CONTROL  SYSTEM 


MEASUREMENTS 


Figure  3.  Control  Process 


The  control  system  design  problem,  stated  simply,  was  to  find  an 
algorithm  to  convert  the  commands  and  sensor  measurements  into  the 
proper  control  settings  to  make  the  aircraft  fly  according  to  the 
cor  <nand s  . 

Due  to  cost  constraints,  the  sensors  available  for  the  project 
were  limited  to  a  yaw  rate  gyro,  a  lateral  accelerometer,  a  heading 
indicator,  an  altimeter,  and  an  airspeed  meter.  This  is  a  very  limited 
set  of  measurements  considering  the  job  required;  therefore,  it  was 
feared  that  design  of  a  satisfactory  control  system  might  prove 
impossible.  For  instance,  note  that  there  are  no  position  gyros.  This 
means  that  the  control  system  is  required  to  roll  the  aircraft  in  and 
out  of  turns  with  no  feedback  whatsoever  as  to  the  bank  angle  of  the 
aircraft.  Almost  all  three-axis  autopilots  in  use  today  have  position 
gryos  for  yaw,  pitch,  and  roll  angle  measurements. 

For  design  purposes,  the  measurements  available  to  the  flight 
control  system  were: 

n  lateral  I oad  factor 
v 

r  yaw  rate 


head i n g 
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h  altitude 
V  airspeed 

h  rate  of  climb  (to  be  derived  by  numerically 
differentiating  h) 

I V .  Des i gn  Procedure 

"Classical"  control  theory  methods  (LaPlace  tranforms  and  root 
loci)  were  used  in  the  design  of  this  control  system  as  opposed  to 
"modern"  optimal  control  techniques  (although  the  state-space  matrix 
representation  was  borrowed  from  modern  theory). 

Two  important  simplifications  were  made  and  carried  throughout  the 
entire  design  procedure.  The  first  was  that  the  control  system  was 
continuous  instead  of  discrete.  (Recall  that  the  control  system  is  to 
be  implemented  by  a  digital  computer.)  This  assumption  is  not  too 
unreasonable  provided  the  cycle  time  of  the  computer  is  quite  a  bit 
faster  than  the  aircraft  response.  The  second  simplification  was  that 
of  perfect  sensors.  This  means  that  the  measurements  of  the  aircraft's 
state  provided  to  the  control  system  are  true  and  uncorrupted  by  noise. 
This  simplification  may  or  may  not  be  valid,  depending  on  the  quality 
of  the  sensors.  As  is  the  case  in  many  feedback  control  design 
problems,  the  sensors  were  of  minor  concern  during  the  design  process, 
but  their  performance  will  make  or  break  the  flight  control  system  when 
it  is  implemented  in  the  actual  aircraft. 

The  design  procedure  consisted  of  four  major  steps: 

1.  Determining  the  control  system  structure 

2.  Formulating  the  complete  system  linear  model 

3.  Selecting  the  control  system  gains 

4.  Checking  the  control  system  performance  in  a 
nonlinear  simulation. 
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A.  Control  System  Structure 

The  control  system  structure  or  framework  was  arrived  at  by 
paralleling  the  way  a  pilot  f 1 ys  an  aircraft  in  instrument  conditions 
without  attitude  information.  (This  type  of  flying  is  referred  to  as 
"needle,  ball,  and  airspeed"  flying  and  is  normally  only  done  as  an 
emergency  procedure  following  an  attitude  indicator  failure.)  The 
structure  is  shown  in  the  block  diagram  in  Figure  4. 


Figure  4.  Flight  Control  System  Structure  (Basic) 

In  the  block  diagram  (Figure  4)  the  squares  with  K' s  in  them 
function  to  multiply  the  input  signal  by  the  constant,  K,  which  is 
called  a  gain.  The  circles  represent  summers  that  subtract  and  add 
signals.  The  subscript,  c,  denotes  a  commanded  value,  as  before.  The 
e  subscript  denotes  an  error  value  which  is  the  difference  between  the 
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commanded  value  and  the  actual  value.  The  trim  control  values  were 
added  to  allow  the  control  system  to  operate  relative  to  some  reference 
steady-state  condition  and  also  to  allow  for  compensation  of  items 
impossible  to  account  for  in  the  math  model,  such  as  wing  warp, 
propeller  torque,  etc. 

To  see  how  the  system  works,  let's  look  at  the  aileron  (<5  ) 

control  loop  which  is  driven  by  heading  ( ii) .  Suppose  is  015  degrees 

(.262  rad.)  but  the  actual  heading,  v,  is  0  degrees  or  north.  then 

is  .262  rad.  This  will  result  in  a  commanded  yaw  rate  (r  )  of 

.  x  K  ■.  rad/s.  If  the  aircraft  currently  has  zero  yaw  rate,  rg  will 

also  be  v  x  K  rad/s.  From  the  diagram,  6  will  now  be  r  x  K  more 
e  0  a  e  r 

than  .  This  will  cause  the  aircraft  to  bank  to  the  right  and 

atrim 

develop  a  positive  yaw  rate.  As  r  approaches  r  ,  rp  gradually 
decreases  and  the  aileron  deflection  is  reduced.  The  aircraft  is  now 
turning.  As  approaches  i  ,  rc  decreases  and  rp  goes  negative.  This 
deflects  the  ailerons  in  the  opposite  direction  to  gradually  roll  the 
aircraft  out  as  $  is  approached.  It  is  instructive  to  note  that  r  is 
not  exactly  equal  to  the  rate  of  change  of  .  However,  it  is  close  for 
small  bank  angles.  (The  same  situation  exists  for  the  pilot  flying 
needle,  ball,  and  airspeed.)  Heading  is  usually  controlled  much  more 
effectively  through  bank  angle  but  that  measurement  is  not  available. 
The  other  loops  work  in  much  the  same  fashion.  The  rudder  is  used 
solely  to  keep  the  aircraft  coordinated,  that  is,  to  null  out  any 
lateral  load  factor.  The  pilot  does  this  the  same  way  by  watching  the 
ball  (skid/slip  indicator). 

The  control  structure  was  later  modified  to  include  a  yaw  damper 
when  the  aircraft  equations  showed  a  highly  oscillatory  dutch  roll  mode 
(lightly  damped  oscillations  in  yaw).  As  can  be  seen  in  Figure  5,  this 
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modification  consists  of  feeding  the  yaw  rate  back  through  the  rudder 
to  damp  out  the  yaw  oscillations.  The  washout  component  is  a  high  pass 
filter  that  prevents  the  rudder  from  fighting  the  commanded  yaw  rate. 
Although  the  yaw  damper  gain,  Kyo  ’  was  ‘nc^udec*  ‘n  the  following 
analyses,  the  washout  filter  was  not. 


With  the  control  structure  defined,  the  problem  remaining  was  to 
determine  the  values  of  the  gains. 

B.  Linear;zed  System  Model 

The  method  used  to  determine  the  control  gains  required  a 
linear  system  model,  that  is,  a  system  described  by  linear  differential 
equations.  The  linear  model  for  the  aircraft  was  derived  using  the 
methods  of  Reference  1.  The  necessary  aerodynamic  data  was  obtained  by 
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wind  tunnel  testing  conducted  in  the  Air  Force  Academy's  low  speed  wind 
tunnel  by  Cadets  First  Class  Charles  T.  Myers  and  Daniel  A.  Draeger  and 
by  estimation  methods  contained  in  Reference  5.  The  required  mass  data 
was  obtained  by  Cadet  First  Class  Thomas  A.  O'Berg  through  inertia 
testing  of  the  actual  aircraft.  All  of  the  aircraft  data  used  are 
contained  in  Appendix  B.  The  actual  aircraft  model  is  assembled  in 
Appendix  C. 

The  aircraft  equations  of  motion  were  linearized  about  a 
steady-state  flight  condition  of  straight  and  level  at  7500  feet, 
standard  day,  and  a  true  airspeed  of  73.33  ft/s  (50mph).  The  usual  two 
independent  sets  of  coupled  equations  were  obtained  and  are  listed 
below  in  matrix  form. 

Long i tud i na I : 


1 

-3.733  0 

.01 189 

-1.341 

c 

- 

-.2033 

0 

e 

1  3 

3 

0 

0 

.3 

0 

3 

■j 

= 

0  -32.2 

-.2566 

M.Oi 

; 

J 

* 

-.0257 

5.63 

a 

.9193  0 

-.31169 

-4.642 

1 

- . 305  7 

3 

n 

1 - 

o 

vx 

0 

-73.33 

c 

0 

3 

q 

=  pi tch  rate  in 

r  ad  /  s 

5  r 


9  =  pitch  angle  change  from  steady-state  in  rad. 
u  =  change  in  velocity  component  along  longitudina 
axis  in  ft/s 

a  =  angle  of  attack  change  in  rad. 
h  r  altitude  change  from  steady-state  in  ft. 
g  =  elevator  deflection  change  from  trim  in  deg. 
gy  =  throttle  change  from  trim  in  brake  horsepower 
( BHP) 

Lateral-directional : 


-.5366  -.000278  -.9875  .439 

-33.18  -8.55  2.41  0 


Pi 


0  .00274 

10.14  .49 
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where:  g  =  sideslip  angle  in  rad. 
p  =  roll  rate  in  rad/s 
r  =  yaw  rate  in  rad/s 
<P  =  roll  angle  in  rad. 

ip  =  heading  change  from  steady-state  in  rad. 
g a  =  aileron  deflection  change  from  trim  in  deg. 

$  =  rudder  deflection  change  from  trim  in  deg. 

Note  that  velocity,  altitude,  and  heading  (u,  h,  ip  )  are  now  pertubation 
values  from  the  reference  flight  condition.  For  example,  an  altitude 
of  7000  feet  would  be  represented  by  h  =  -500  ft.  Figure  5  is  still 
perfectly  valid  if  the  command  values  are  also  pertubation  values. 

Each  of  these  equations  can  be  written  in  the  standard  state 
variable  form  as: 

x  =  [A]  x  +  [B]  u  (3) 

where  x  is  defined  as  the  state  vector  and  u  is  defined  as  the  control 
vector.  The  components  of  the  state  vector  (such  as  h  and  r)  are 
called  states.  [A]  and  [B]  are  referred  to  as  the  system  and  control 
matrices,  respectively. 

These  matrix  equations  describe  the  aircraft  only.  It  is 
necessary  to  combine  the  control  system  and  the  aircraft  or  "close  the 
loop".  Referring  to  Figure  5,  it  can  be  seen  that  the  control  inputs 
will  be  linear,  algebraic  combinations  of  the  measurements  and  the 
commands.  (As  mentioned  earlier,  we  are  ignoring  the  washout  filter.) 
The  measurements,  in  turn,  are  all  either  states  themselves  or  linear, 
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algebraic  combinations  of  the  states.  Equation  1  shows  that  this  is 
true  for  the  measurement  h.  This  is  also  true  for  n^  .  V  is  closely 
approximated  by  u.  It  is  possible,  then,  in  both  the  longitudinal  and 
1  a t er a  1 -d i r ec t i ona 1  cases  to  express  the  controls  as: 


u  =  -  [F]  x  +  [ B ' ]  uc 


(4) 


where  [F]  is  called  the  feedback  matrix  and  [B1  ]  is  called  the  input 

T 

matrix.  The  vector  uc  is  the  command  which  is  defined  as  [hc,Vc]  for 

j 

the  longitudinal  case  and  [ ^ ,  O]  for  the  1  a t er a  1 -d i r ec t i ona 1  case. 

From  Figure  5  and  Equation  1  and  2  the  feedback  and  input  matrices 
were  determined  to  be  as  follows: 

For  the  longitudinal  model: 


[F] 


0  7  3 . 3  3 

0  0 


0  -73.33K^  KhKh- 

Kv  0  0 


(5: 


[B1  ] 


KhKh 


0 

K 


v 


(6) 


n  is  lateral  load  factor  which  is  equal  and  opposite  to 
no^g r av i t a t i ona 1  lateral  acceleration,  normalized  to  the  acceleration 
of  gravity,  g.  It  can  be  shown  from  Reference  1  and  Appendix  C  that  for 
the  linearized  approximation, 

n  =  1.222H  +  .00Q634p  -  .0285  r  -  .00624 6  ,  which  is  a  linear, 
algebraic  combination  of  the  states.  (6r  is  a  combination  of  the 
states  only,  since  there  are  no  command s r t ha t  feed  into  the  rudder.) 
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For  the  I  a t e r a  1 -d i r ec t i ona I  model: 


rH  1.22K 


1  - .  0  0  6  2  4  K 


B’]  =  K  0 
0  0 


.  000634K 


1  -  .  00624K 


-  .  0  2 85K 


1  -  .  00624K 


By  substituting  (4)  into  (3),  we  obtain 


t  =  [A]  Y  +  [B]  [-[F]  x  +  [ B 1  ]  TTcJ. 


Combining  terms  gives 


x  =  [A-  BF]  x  +  [BB']  u 


The  complete  linearized  system  model  (or  closed-loop  system), 
then,  consists  of  two  independent  equations  of  the  form  of  Equation  10 
--  one  for  longitudinal  motion  and  one  for  1  a t er a  1  - d i r ec t i ona 1  motion. 
The  vectors  and  constant  matrices  have  been  defined  for  each  case. 

C.  Determination  of  Control  System  Gains 

The  core  of  the  feedback  control  problem  is  the  selection  of 
the  control  gains.  In  this  project,  that  means  finding  values  of  , 
K^,  Ky,  K ^ ,  Kr ,  Kny,  and  Ky^  that  give  satisfactory  response  to 
aircraft  heading,  altitude,  and  airspeed  commands  in  the  presence  of 
disturbances  sur-h  as  wind  gusts. 

The  response  of  any  linear  dynamic  system  is  characterized  by  the 
roots  of  its  characteristic  equation  (also  called  system  poles  or 
eigenvalues).  Root  loci  were  used  in  this  project  to  select  the  gains 
A  root  locus  is  a  complex  plane  plot  showing  how  a  pole  varies  as  a 
gain  is  changed.  By  way  of  review,  the  sign  of  the  real  part  of  the 
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pole  indicates  the  stability  of  its  associated  mode  (-  for  stable,  + 
for  unstable,  0  for  neutral)  and  the  imaginary  part  indicates  the 
oscillation  frequency  (aperiodic  if  real).  Complex  poles  occur  in 
conjugate  pairs.  The  magnitude  of  the  pole  (distance  from  origin) 
indicates  the  speed  of  the  associated  mode.  For  example,  in  Figure  6, 
poles  1,  2,  and  5  are  aperiodic.  Poles  1  and  5  are  stable  while  2  is 
unstable.  The  mode  associated  with  pole  5  will  die  out  faster  than  the 
one  associated  with  pole  1.  The  complex  conjugate  pairs  3  and  4 
represent  oscillatory  modes.  Mode  3  is  stable  and  4  is  unstable. 


IMAGINARY 


REAL 


Figure  6.  Poles  on  the  Complex  Plane 


It  can  be  shown  (Ref.  4)  that  the  characteristic  equation  of  a 
system  described  by  Equation  3  is  det(A[I]  -  [A])  =  0  where  \  is  an 


arbitrary  scalar  number  and  [I]  is  an 
the  equation  for  >  are  the  poles.  We 
the  closed-loop  system,  Equation  10. 
solving  det('[I]  -  [A-BF]  )  =  0.  The 


identity  matrix.  The  solutions  of 
are  interested  in  the  poles  of 
These  can  be  determined  by 
F  matrix,  of  course,  depends  on  the 
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gains. 

The  root  loci  were  constructed  by  varying  the  gains  and  solving 
for  the  poles  after  each  change.  This  process  had  to  be  computerized. 
Separate  programs  were  written  on  an  Apple  microcomputer  for  the 
longitudinal  and  1 ateral -d i rec t i ona 1  cases.  The  results  are  discussed 
in  the  next  section. 

D.  Nonlinear  Simulation 

A  12  degree  of  freedom  nonlinear  Big  Stick  simulation  program 
was  written  for  the  Burroughs  6900  computer  at  the  Air  Force  Academy  by 
Cadet  First  Class  Daniel  A.  Draeger.  A  hard  copy  is  included  in 
Appendix  D.  This  simulation  provides  a  much  more  accurate  mathematical 
model  of  the  aircraft  than  the  linearized  equations  which  were  used  to 
determine  the  control  gains.  Basically,  the  program  numerically 
integrates  the  nonlinear  aircraft  equations  of  motion  from  Reference  1, 
modified  to  include  the  control  system,  and  plots  out  any  of  the  state 
variables  versus  time.  The  nonlinear  equations  do  not  decouple  into 
longitudinal  and  1  a t er a  1  - d i r ec t i ona 1  sets. 

The  simulation  was  run  to  see  how  the  control  gains,  selected 
under  the  linear  assumption,  would  actually  perform  in  the  real, 
nonlinear  world.  The  effects  of  such  elements  as  control  deflection 
limits  and  changes  in  air  density  could  be  observed.  The  simulation 
was  also  used  to  check  the  limiting  values  for  r£  and  hc.  These  are 
cutoff  values  which  had  to  be  incorporated  into  and  to  prevent 
the  aircraft  from  stalling  itself  out  or  entering  a  steep  dive  in  the 
case  of  a  big  change  in  hc  or  rolling  inverted  when  t|i  changed.  A  side 
benefit  of  the  simulation  is  that  it  provides  a  check  on  the  previous 
calculations.  Performance  should  be  close  to  the  linear  prediction 
around  the  steady-state  condition. 
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The  simulation  results  for  the  selected  gains  are  presented  in  the 
nex t  section. 

V .  Results 

A.  Root  Loci 

In  multi-loop  feedback  systems,  such  as  the  one  being  dealt 
with  in  this  report,  the  gains  affect  the  system  poles  in  an 
interrelated  and  complex  fashion.  A  "shotgun"  (trial  and  error) 
approach  was  used  to  initially  find  a  neighborhood  of  gains  that 
appeared  to  give  reasonable  poles.  The  gains  were  then  varied  in  a 
more  systematic  fashion  to  refine  the  gain  selection.  Gains  were 
selected  on  the  basis  of  the  speed  and  stability  of  the  resulting 
poles.  This  was  done  separately  for  the  longitudinal  and 
1 ateral -di rect ional  cases. 

1 .  Longi tudinal  case 

The  three  longitudinal  gains  selected  were: 

- . 2  deg/ft/s,  Kh=.2  ft/s/ft,  and  Ky=.l  BHP/ft/s.  Figure  7  shows  a 
segment  of  the  longitudinal  root  loci  plot  that  indicates  how  the  poles 
are  affected  by  the  gains  around  the  selected  values.  Only  the  upper 
left  quadrant  of  the  complex  plane  is  shown  since  any  values  in  the 
right  half  plane  are  unstable  and  unacceptable  and  since  the  bottom 
half  plane  is  a  mirror  image  of  the  top.  The  actual  values  of  the 
poles  plotted  are  contained  in  the  computer  printout  in  Appendix  E. 
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B.  Simulation 

Simulations  with  the  selected  gains  were  run  using  the 
nonlinear  simulation  program  of  Appendix  D.  Limitations  on  control 
deflections  and  on  some  of  the  gains  were  added  for  reasons  discussed 
in  the  last  section.  Four  flight  maneuvers  were  simulated.  The 
initial  conditions  for  each  maneuver  were  the  steady-state  reference 
condition  (straight  and  level  with  41  =0,  V=7 3.33  ft/s,  and  h=7500  ft). 
The  four  maneuvers  were: 

(1)  level  turn 

(2)  straight  cl imb 

(3)  level,  straight  acceleration 

(4)  combination  turn,  descent,  and  deceleration 
1 .  Non  1  inear  1 imi ts 

The  following  control  limits  were  used  in  the  simulation 
based  on  estimated  aircraft  limits: 


ntrol 

min . 

max. 

6e 

-15  deg 

15  deg 

<ST 

0  BHP 

3  BHP 

6  a 

-15  deg 

15  deg 

6  r 

-15  deg 

15  deg 

Limits  on  commanded  rate  of  climb  (hc)  were  selected  as  ±11.5  ft/s. 

This  corresponds  to  climb  and  descent  angles  of  nine  degrees.  The 
climb  would  require  about  2  BHP  at  steady-state.  The  commanded  turn 
rate  limit  was  selected  as  ±.12  rad/s.  This  equates  to  a  bank  angle  of 
about  15  degrees  and  a  turn  radius  of  about  625  feet. 

2.  Level  turn  simulation 

A  level  turn  to  a  heading  of  090  deg  (1.57  rad)  was 
commanded.  The  aircraft  response  is  plotted  in  Figures  10-13.  The 
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slow  oscillatory  mode  predicted  in  root  loci  Figure  9  can  be  seen  in 
Figures  10  and  11  as  the  aircraft  overshoots  the  commanded  heading  and 
corrects  back.  Note  how  this  1  a t er a  1 -d i r ec t i ona 1  mode  couples  over 
into  the  longitudinal  variables  h  and  u  in  this  nonlinear  simulation 
(Figures  12  and  13).  The  aircraft  does  stabilize  on  its  new  heading  in 
less  than  a  minute  after  overshooting  about  20  degrees.  The  altitude 
is  held  within  25  feet  and  the  airspeed  within  2  ft/s. 

3.  Straight  climb  simulation 

A  straight  climb  to  8000  feet  was  commanded.  The 
aircraft  response  is  plotted  in  Figures  14  and  15.  Figure  14  shows 
that  the  aircraft  enters  a  steady  climb  of  approximately  10  ft/s. 
Approaching  8000  ft,  it  levels  off  nicely  with  little  or  no  overshoot. 
Notice  from  Figure  15  that  the  control  system  is  unable  to  keep  the 
airspeed  at  the  commanded  value  during  the  climb.  Instead,  the 
airspeed  stabilizes  at  a  value  of  about  6  ft/s  below  the  commanded 
value.  Due  to  the  fact  that  there  is  no  acceleration  feedback,  some 
steady-state  error  in  velocity  will  always  exist  when  not  at  the 
reference  condition.  (In  the  control  literature,  this  is  known  as 
"type  zero"  behavior  and  is  predictable  from  the  linear  analysis.) 

When  the  aircraft  levels  off,  the  control  system  is  able  to  increase 
the  airspeed  back  up  to  approximately  the  commanded  value. 

4.  Level  acceleration  simulation 

An  acceleration  to  90  ft/s  was  commanded.  The  aircraft 
response  is  plotted  in  Figure  16.  As  in  the  straight  climb,  a 
steady-state  error  of  about  6  ft/s  exists. 

5.  Mu  1 t i -command  simulation 

To  see  how  the  aircraft  would  react  to  simultaneous 
commands,  inputs  of  ^  =-1.05  rad  (300  deg),  h  =7000  feet,  and  V  =60 
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ft/s  were  made.  The  results  are  shown  in  Figures  17-20.  As  Figure  17 
shows,  the  heading  again  overshoots  but  returns  to  the  commanded  value 
From  Figure  18  we  can  see  that  the  aircraft  levels  off  at  a  value 
slightly  below  7000  feet.  Figure  19  shows  that  the  velocity  goes 
through  some  gyrations.  It  stabilizes  at  about  70  ft/s  in  the  descent 
and  then  reduces  to  about  62  ft/s  when  the  aircraft  pulls  out  of  its 
dive.  Figure  20  is  a  plot  of  the  throttle  activity  during  this 
maneuver.  It  shows  the  throttle  at  idle  from  the  5  to  13  second  point 
The  power  comes  back  in  to  about  .05  BHP  during  the  descent  and  then 
comes  up  to  about  .8  BHP  after  the  aircraft  levels  out. 


222 


r  -V  ** 


USAFA-TR 


o 

o 

<o 

o 

to 


O 
J-  O 


STATE  VARIABLE  (u,  ft/s) 

76.00  78.00  80,00  82.00  84.00  86.00 


USAFA-TR-85-2 


Figure  16.  Level  Acceleration  Simulation  -  u 
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VI  .  Reconnmendat  i  on  s 

Based  on  the  previous  analysis,  the  following  recommendations  are 
made  for  the  Big  Stick  flight  control  system: 

1.  That  the  control  structure  of  Figure  5  be  used. 

2.  That  the  following  control  gains  be  used: 

Kh  =  .2  (ft/s)/ ft 

=  -.2  deg/(ft/s) 

Ky  =  .1  BHP/(ft/s) 

K  ^  =  .1  (rad/s)/rad 

=  1  deg/(rad/s) 

Kny  =  30  de8/8 

KyD  =  -20  deg/(rad/s) 

3.  That  hc  be  limited  to  ±11.5  ft/s  and  rc  be  limited  to  +  12 
rad/s . 

<i.  That  a  washout  filter  be  included  in  the  yaw  damper  feedback 
loop  to  pass  the  dutch  roll  frequency  (4  rad/s)  and  attenuate 
the  steady-state  commanded  yaw  rate. 

5.  That  the  trim  control  settings  be  determined  from  flight  test 
at  the  steady-state  flight  condition  of  straight  and  level  at 
73.33  ft/s  and  7500  feet. 

In  implementing  the  flight  control  system,  the  same  sign 

conventions  must,  of  course,  be  used  as  are  used  in  the  aircraft 

equations  of  motion.  These  are  (from  Appendices  C  and  D  and  Reference 
1): 

s  positive  for  trailing  edge  down 

(>t.  positive  for  trailing  edge  left 

.  positive  for  right  up  and  left  down,  *  =  (*  +6a  )/2 

3  aL  aR 

r  positive  to  right 
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pos  i  t  i  ve  to  r  ight 
n  pos itive  to  right 
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APPENDIX  B-  Big  Stick  Aerodynamic  Coefficients,  Stability  Derivatives 
and  Mass  Data 

I.  Miscellaneous  Aircraft  Data: 


Weight  30  1  bp 

Reference  area  (wing,  S)  15.6  ft^ 

Chord  (c)  1.784  ft 

Span  ( b )  8.73ft 

Propeller  efficiency  ( n )  .7 


2.  Wind  Tunnel  Data  (Reference  2): 


CQ  =  .061 
o 

K  =  .086 


CL  =  4.0/rad 
a 

C  =  +. 005/deg 

6e 

c.o  =  -0093 

C„  =  -.29/rad 
m 

a 

C  =  -.01/deg 

m . 

6e 

Cn  =  .0003/deg 


o 

.00235/deg 
- .46/rad 

.025/deg 

.00105/deg 

.17/rad 

-.002/deg 

-.00144/deg 
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3.  Estimated  Aerodynamic  Data  (Reference  1  and  5): 


c*  = 
u 

CL  =  CD 
u  u 

c  -  o 

fflT 

a 

o 

3 

II 

-4.0 

% 

=  0 

T6 

Cm  ■ 
q 

-11.0 

-.046 

'V 

1.66 

V 

-.03 

V 

4.16 

%  = 

-.078/rad 

CT 

x 

u 

-.22 

c*  = 

p 

-.36 

V 

-.004 

CJ>  = 

r 

.096 

Cv  * 

.18 

CT  = 

.066/BHP 

T 

♦wind  tunnel  data  not  used  because  model  had  no  wing  dihedral 


4.  Mass  Data  (Reference  4); 

rolling  moment  of  inertia  Uxx)  =  1.7  slg-ft2 

pitching  moment  of  inertia  (Iyy)  =  6.8  slg-ft2 

yawing  moment  of  inertia  d2z)  =  9.3  slg-ft2 

product  of  inertia  (Ixz>  =  0 

notes  these  inertia  terms  are  relative  to  body  axes 


APPENDIX  C~  Aircraft  Linearized  Equations  of  Motion 

The  methods  of  Reference  l  were  used  to  linearize  the  aircraft 
equations  of  motion  about  a  steady-state  condition  of  coordinated, 
straight  and  level  flight  at  7500  feet  and  73.33  ft/s  (50  mph).  The 
angle  of  attack  required  for  this  condition  was  7.2  degrees  relative  to 
body  axes.  Since  the  linearized  model  uses  stability  axes 
(longitudinal  axis  parallel  to  the  steady-state  relative  wind),  the 
inertia  terms  lxx,  \ zz,  and  IX2,  of  Appendix  B  had  to  be  transformed 

from  body  axes  to  stability  axes. 
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The  remaining  necessary  coefficients  and  derivatives,  listed  below, 
were  obtained  by  linearizing  the  drag  polar  around  the  steady-state 
condition  and  by  assuming  that  the  propeller  thrust  acts  through  the 
aircraft  center  of  gravity. 


C 


=  .073 


CD  =  .259/rad  c  =  0 

a  T1 

c4  =  -57«  C»T  *  0 


CT  =  .073 
*1 


From  Reference  1,  the  linearized  longitudinal  equations  of  motion  are 
(using  the  definitions  from  Table  C-l): 


u  =  -g  e 


V 


+  XT  u  + 


V 


V* 


+  X 


Vt 


u 


In  addi t ion, 


-  U,q  =  Z,,u  +  Z  a  +  Z*a  +  Z_q  +  Z.  fie 
1  u  a  a  q  6 

e 

q  =  Muu  +  Mj  u  +  MQa  +  My  a  +  M^a  +  M^q  +  fie 
the  following  approximations  were  added: 


e=  q 

h  =  Uy  e  -  U1a 

The  values  were  substituted  in  and  the  equations  were  manipulated 
algebraically  to  obtain  first  order  matrix  form. 


-3.733 

0 

.01189 

-1.341 

0 

q 

-.2033 

0 

1 

0 

0 

0 

0 

0 

0 

0 

= 

0 

-32.2 

-  .2566 

10.01 

0 

u 

+ 

-.0257 

5.63 

.9193 

0 

-.01169 

-4.642 

0 

a 

-.0057 

0 

. 

0 

73.33 

0 

-73.33 

0 

h 

0 

0 
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Table  C-l 


qlSCL. 


-q iS(Cq^  +  2Cq^ ) 


Xu  - 


mil. 


-(sec-1 ) 


^SCCy  +  2Ct  ) 


x  . - !H - li - (sec-1) 

u  mU, 


-Qls(CD  -  CL  ^ 
a  1 


X  =- 

a 


-(ft  see-2) 


"qlSCD. 


— - (ft  sec-2deg-1) 


q1S(.CL  +  2C,  ) 


u 


Zu  =  - 


•  (sec-1) 


mil. 


qiS(CL  +  Cq  ) 
a  i 


.(ft  sec-2deg-1) 


^e  = 


qlS^Cmu  +  “m^ 


M  - - — (ft-1sec-1) 


Iyyui 


M 


qlS^CmT  +  2CmT  5 

_ _ ^ - 1 — (ft-1sec-1 ) 

Tu 


yy  i 


q ,  ScC_  0 

=  1  a  _ /sec-2) 


yy 

qlS^Cm. 


V  » 


Jt _ (sec-2) 


a  I 


yy 


qlS^Cm 


- 


2Iyyul 


.(sec-1 ) 


-(ft  sec-2) 


q^Sc2C 


m 


Mq  - 


2Iyyui 


-(sec-1) 


Q  i  SC .  C 
a 


Zo  = 


2mU. 


-(ft  sec-1) 


q  ^ScCp 


O  1 

(sec-  deg-  ) 


e  I 


yy 


qlscT 


QlSC,  c 

q  _ i , 

Z  - - - - (ft  sec  1 ) 


XT 


-(ft  sec-2BHP-1 ) 
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From  Reference  1,  the  linearized  1  a  ter al -d i r ec t i ona 1  equations  of 
motion  are  (using  the  definitions  from  Table  C-2): 


Ui®  +  Uxr  =  g0  +  Yg0  +  Y  p  +  Yrr  +  Yfi  «a  ♦  Y.  5 

a  r 

P  -  All  =  Lb6  *  LPP  *  Lrr  *  L,  «,  .  L  Sr 

a  r 

r  "  BiP  *  V  +  nt  0  +  n  p  +  Nrr  +  n6  6a  +  n.  6 

p  a  r 


In  addition,  the  following  approximations  were  added: 

*  P 

i  =  r 

The  values  were  substituted  in  and  the  equations  were  manipulated 
algebraically  to  obtain  first  order  matrix  form. 
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q  i sc  > 


(ft  sec"2) 


qiSbCyp 

Y  = -  (ft  sec"1) 

p  2mU, 


qlsbc 


Y_  = 


2mU, 


(ft  sec”1) 


q  i  SC 


1  Yj 


(ft  sec“2deg_1) 


a  m 


qlSC5 


(ft  sec"2deg"1) 


qjSbC  it 


B 


L8  = 


(sec-2) 


XX 


q^b2^ 


L_  = 


p  .1, 

—  (sec  A) 


2IWWU 


XX  1 


q1Sb2c‘r 


L_  = 


21  WU 


(sec”1 ) 


xxul 


^ ^SbC  jt 


O  1 

(sec"  deg"  ) 


'a  I 


xx 


q^SbC  jt 


ry  _  1 

(sec”  deg”  ) 


xx 


qlSbCn£ 


NB  = 


(sec-2) 


zz 


qlsbCn 


Nt  =  -  (sec-2) 

•  t 


B  I 


zz 


N_  = 


P  21Z2U1 


P  -1. 

—  (sec  A) 


q1Sb2C„_ 


N_  = 


(sec-1 ) 


21. .U 


zz  1 


qlsbc 


(sec"2deg-1 ) 


’a  I 


zz 


qlSbC 


*y  1 

(sec"  deg"  ) 


‘r  I 


zz 
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APPENDIX  D-  Simulation  Program 

This  nonlinear  simulation  program  was  written  in  FORTRAN  77  for  the 
Burroughs  6900  computer  at  the  Air  Force  Academy  by  Cadet  First  Class 
Daniel  A.  Draeger.  The  program  numerically  integrates  the  six  aircraft 
equations  of  motion  for  u,  v,  and  w  (velocity  components  in  body  axes) 
and  for  p,  q,  and  r  (rotational  rates  in  body  axes).  Six  kinematic 
equations  are  also  integrated  to  get  the  Euler  angles  0  ,  9  ,  and0  ;  and 
displacements  in  the  earth  axes:  altitude  (h),  north  distance  (N),  and 
east  distance  (E).  The  results  of  simulation  runs  are  delivered  in 
graphical  form  using  plotting  routines  not  included  in  tbe  listing. 

The  aircraft  data  are  read  in  from  separate  data  files,  a  feature  which 
allows  the  program  to  be  used  to  simulate  other  aircraft.  The  twelve 
equations  solved  by  the  program  are  listed  below  in  general  form 
(Ref erence  1 ) . 

m  (u-vr+wq)  =  forces  in  x  direction 

m  (v+ur-wp)  =  forces  in  y  direction 

m  (w-uq+vp)  =  forces  in  z  direction 

‘xxP  -  'xzf  -  'xzrq  +  (,zz-Iyy)rq  =  moment  about  x  axis 

I y yq  +  (Ixx-,zz)Pr  +  1 X2  (p  -r  )  =  moment  about  y  axis 

*zzr  ‘  *xzp  +  ^yy_Ixx^  pq  +  ^z^  s  moment  about  z  axis 


♦  =  (q  sin  0 +  r  cos  0  )  sece 

0  =  q  cos  0  -  r  sin  0 

0  =  p  +  q  sin  0  tan  0  +  r  cos  0  tan  9 


11  =  (cos  i/-  cos9 )U  +  (cos  'I'  sin  9  sin  0  -  sin  *  cos0)V  +  (cos  *  sin  9  cos  0  +  sin  *  sin0)W 
E  =  (sin  0  cos  0 JU  +  (sin  <l>  sin 9  sin0  +  cos  <l>  cos0)v  +  (sin  1  sin  0  cos  a  -  cos  sin0)w 
h  =  (sin9)u  +  (-cos  0  sin^Jv  +  (-cosecos0)W 
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100  4RESET  FREE 

200  FILE  4(TITLE»"CNSTNT.''.KIND»"DISK*,FILETYPE*7) 

300  FILE  5(TITLE»"SYS DAT.", KIND- " DISK", FILETYPE-7) 

400  FILE  6 (K  IND*"PR INTER" ) 

500  FILE  7<KINO-"PEMOTE"> 

550  FILE  8  (KIND«"REMQTE“ ) 

600  4INCLUDE  "PLOT/ INCLUDE/FORTRAN. “ 

700  4INCLUDE  "  (AER04S7)  DFAN/PLOTTER" 

800  C  *»**»»»#***»*»*»»******»*#***#*#*******###**************#*********#** 

900  C  #**FOLLOWING  IS  A  LIST  OF  VARIABLE  NAMES  AND  USES  IN  MAIN  PROGRAM**#* 

1000  C  DT  ■  STEP  SI2E  IN  NUMERICAL  INTEGRATION 

1100  C  T  »  TIME  THE  INTEGRATION  IS  CURRENTLY  AT 

1110  C  X  •  STATE  VAR.  ARRAY-  U. V.W.P.Q.R.PSI .THE, PHI ,N,E,H 

1120  C  DX  ■  DERIVATIVES  OF  X 

1130  C  DE.BHP.DA.DR  »  AIRCRAFT  CONTROLS 

1140  C  DETRM.BHPTRM.ETC.  -  TRIM  SETTINGS  FOR  CONTROLS 

1150  C  AKHDOT.AKV.ETC.  ■  AUTO  PILOT  CONTROL  GAINS 

1160  C  IX  »  VARIABLE  OR  CONTROL  TO  BE  PLOTTED 

1170  C  ALTC, VELC.PSIC  ■  COMMANOS  TO  FLY  (ALT, VEL, HEADING) 

1200  C  TMAX  ■  THE  TIME  LIMIT  FOR  THE  INTEGRATION 

1300  C  PLTSPC (5)  ■  SELECTION  FOR  LP  (1),  HP  (2),  BOTH  (3) 

1400  C  NSPP  «  NUMBER  OF  POINTS  PLOTTED 

1500  C  NPTS  •  CURRENT  POINT  NUMBER  COUNTER 

1600  C  NDIM  »  12.  THE  DIMENSION  OF  X  AND  DX 

1700  C  ICNT  ■  CURRENT  LOOP  NUMBER  COUNTER 

1800  C  XP  >  ARRAY  ST0RIN6  THE  DESIRED  X  POINTS 

1900  C  TP  «  ARRAY  STORING  THE  DESIRED  T  POINTS 

2000  C 


DT 

T 

X 

DX 

DE.BHP.DA.DR 

DETRM.BHPTRM.ETC. 

AKHDOT.AKV.ETC. 

IX 

ALTC. VELC.PSIC 
TMAX 
PLTSPC (5) 
NSPP 
NPTS 
NDIM 
ICNT 
XP 
TP 


»»*t**»»»»* »*•*••»»••••• »••*•»•*•*•»**»•••*•••••»••»••»•»»••••••«*«*• 


2200 
2300 
2400 
2500 
2600 
2650 
2700 
2800 
2900 
3000  C 
3100  C 
3200  C 
3300  C 
3400 
3410 
3420 
3450 
3500 
3600  1 

3650  1 

3700  2 

3800  2 


DIMENSION  X (12) ,0X112) . ER < 12) , XP ( 1 . 250 > , TP ( 1 . 250) .PLTSPC (5) 
CHARACTER*6X TITLE  <4) , YTITLE !4> .LABEL (4) ,FLINE(4) .SLINE (4) ,TLINE<4> 
COMMON  / BLK 12/  CDO  ,  AKD  .  CLO . CLA . CLAD , CLQ , CLDE . CMO . CMA . CMAD , CMQ , 

1  CMDE,CY8,CYR.CYP,CYDR,CLB,CLP.CLR,CLDA,CLDR,CNB, 

2  CNP , CNR ,CNDA .CNDR.S.CB.B, ANP .AIX.AIZ.AIXZ.AM.AIY, 

3  DA.DE.DR.BHP 

DATA  LABEL/"  " , "FLT  HI", "STORY  " , "  "/ 

DATA  XTITLE/"  TIME"."  <SECO’."NDS)  "/ 

DATA  YTITLE/"  STAT'VE  VARI" ,  "ABLE  "/ 


READ  IN  THE  DATA 


READ  <5,»)  (XU)  ,1-1.12) 

READ  (5.*)  DETRM,  BHPTRM,  DATRM,  DRTRM 

READ  (5.*)  AKHDOT, AKV,  AKH,  AKR.  AKPSI.  AKNY ,  AKYD 

READ  (5,100)  DT , TMAX 

READ  (5.200)  IX,  PLTSPC  <5  > 

FORMAT  (2 (FI 0. 4) ) 

FORMAT  (6 (F9. 3) , / , 6 (F9, 3) ) 

FORMAT  (2(15)) 

FORMAT  ("  YOUR  INPUT  DATA  NASi  ") 
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3900  WRITE  (7,225) 

4000  WRITE  (7,110)  < X < I ) . I* 1 , 12) 

4010  WRITE  (7,110)  DETRM,  BHPTRM,  DATRM,  DRTRM 

4020  WRITE  (7,110)  AKHDOT.AKV,  AKH,  AKR,  AKPSI,  AKNY ,  AKYD 

4050  WRITE  (7,100)  DT.TMAX 

4100  WRITE  (7,200)  I  X  , PLTSPC (5) 

4200  READ  (4,#)  CDO . AKD , CLO , CLA , CLAD , CLQ , CLDE , CMO . CMA , CHAD , CMQ , CMDE 

4300  READ  (4,#)  CYB . CYR . CYP , CYDR . CLB , CLP , CLR , CLDA , CLDR . CNB , CNP . CNR 

4400  READ  (4,»)  CNDA . CNDR , S , CB , B . ANP , AI X . AI Z , A I XZ , AM, AI Y 

4500  WRITE  (7,410)  CDO, AKD, CLO, CLA, CLAD, CLQ, CLDE. AM 

4600  WRITE  (7.420)  CMO, CMA, CHAD, CMQ, CMDE. CYB, CYR, CYP, CYDR 

4700  WRITE  (7,430)  CLB. CLP. CLR. CLDA. CLDR, CNB. CNP. CNR 

4600  WRITE  (7,440)  CNDA , CNDR , S , CB , B , ANP , AI X . AI Z , AI X 2 , AI Y 

4900  410  FORMAT ( 1 X . 'CDO* ' ,F6. 3 ,  '  AKD«',F6.3.'  CL0«'.F6.3,'  CLA*  '  ,F6. 3. / , 

5000  1  '  CLAD*  '  ,F6, 3, '  CLQ»',F6.3,'  CLDE*  '  .  F6. 3 ,  '  AM*',F6.3) 

5100  420  FORMAT (IX,  CMO* ' ,F6. 3,  '  CMA«’,F6.3.'  CHAD*  '  ,F7. 3.  '  CMQ*  '  ,F7 . 3 , / . 

5200  1  '  CMDE*' .F6.3,  '  CYB*',F6.3,'  CYR*',F6.3,'  CYP*',F6.3, 

5300  2  '  CYDR* ' , F6 . 3 ) 

5400  430  FORMAT (IX, 'CLB*’.F6.3,  '  CLP*',F6.3,'  CLR*‘,F6.3,'  CLDA* ' , F6. 3 , / , 

5500  1  '  CLDR*' ,F6.3. '  CNB«',F6.3,'  CNP*',F6.3,'  CNR*'.F6.3) 

5600  440  FORMAT (IX, ‘CNDA*’ ,  F6.3,  '  CNDR*  '  ,  F6. 4 ,  '  S«',F6.3,'  CB*',F6.3,/, 

5700  1  '  B* '  ,F6. 3 , '  ANP«',F6.3.'  AIX*',F6.3,'  AIZ*'.F6.3.'  AIXZ*'.F6.3, 

5750  2  / , '  AIY*  '  ,F6. 3) 

5775  WRITE  (7,101) 

5780  101  FORMATUX,  'PLEASE  ENTER  ALTITUDE,  VELOCITY,  AND  HEADING  COMMANDS' 

5785  1  /,'  (UNITS  ARE  FT..  FT/SEC.  AND  RADIANS).') 

5790  READ  (8,#)  ALTC ,  VELC,  PSIC 

5800  C  ##****#*#*«###***♦#*#*******#*#»****♦*#♦»*#****#*»**♦#*******#«**#* 

5900  C  INITIALIZATION  BLOCK 

6000  C  HIIHIHHmmtHHtmmiHHIItHHHmHiHIHHMHimtH 

6100  NSPP  *  INT ( (TMAX/DT) /250. >  +  1 

6150  PLTSPC ( 1 )  *  8 

6155  PLTSPC (2)  »  6 

6160  CALL  PLOTS (1,0,0) 

6170  IF  (IX.GT.12)  GO  TO  230 

6200  XP(l,l)  *  X(IX) 

6210  GO  TO  240 

6220  230  IF  (IX.EQ.13)  XP(l.l)  *  DE 

6230  IF  (IX.EQ.14)  XP  < 1 , 1 )  ■  BMP 

6240  IF  (IX.EQ.15)  XP <  1 , 1 )  *  DA 

6250  IF  (IX. EQ. 16)  XP(1,1)  *  DR 

6300  240  TP(l.l)  *  0. 

6400  NPTS  -  0 

6500  ICNT  «  0 

6600  NDIM  *  12 

6700  DO  250  I  ■  1,12 

6800  250  ER  ( I )  -  0. 

6900  T  *  0. 

7000  C 

7100  C  LOOP  THROUGH  INTEGRATION  USING  RUNGA-KUTTA  METHOD 

7200  C  ####»»#«»#*###**##**»*###»**##»#*»##**#####*###** 

7210  C  •*  CALCULATE  CONTROLS  FOR  EACH  INTEGRATION  STEP  *♦#**•**»»**#**§# 
7220  300  VEL  *  SQRT  (  X(1)*X(1)  +  X (2) *X  (2)  ♦  X ( 3 ) *X ( 3 )  ) 

7230  ANY  «  COS (  X<8)  )*SIN(  X  <  9 )  )  -  (DX (2)  +  X(l)*X(6)  -  X(3)*X(4)  ) 

7232  1  /  32.17 

7235  HDOTC  «  AKH  ♦  ( ALTC-X (12)  ) 

7236  IF  (HDOTC. GT.  11.5)  HDOTC  ■  11.5 

7237  IF  (HOOTC.LT. -11. 5)  HDOTC  —11.5 

7240  OE  *  DETRM  +  AKHDOT#  (HDOTC  -  DX(12>  ) 

7250  BHP  -  BHPTRM  ♦  AKV  ♦  (VELC  -  VEL) 
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7253 
7236 
7257 
7260 
7270 
7273  C 

7280 

7281 

7282 

7283 

7284 

7285 

7286 

7287 
7300 
7400 
7500 
7600 
7700 
7710 
7800 
7810 
7820 
7825 
7830 
7835 
7900 
8000 
8100 
8200 
8300  C 
8400  C 
8300  C 
8600 
8700 
8800 
8900 
9000  C 
9100  C 
9200  C 
9300 
9400 
9300 
9600 
9700 
9800 
9900 
10000 
10100 
10200 
10300 
10400 
10300 
10600 
10700 
10800 
10900 
11000 
11100 
11200 


RC  ■  AKPSI  #  (PSIC  -  X  (7)  ) 

IF  (RC.GT.O. 12)  RC  ■  . 12 
IF  (RC.LT.-O. 12)  RC  «-.12 
DA  -  AKR  •  (RC  -  X (6)  )  ♦  DATRM 
DR  «  -AKNY  #  ANY  -  AKYD  «  X (6)  +  DRTRN 
**  CHECK  CONTROL  LIMITS  f ************* 
IF  (DE.8T.15.)  DE  ■  13, 

IF  (DE.  LT.  -13. )  DE  —15. 

IF  (BHP.BT.3.)  BMP  -  3. 

IF  (8HP.LT.0.)  BHP  ■  0. 

IF  (DA.GT.13.)  DA  ■  15. 

IF  (DA.LT.-13.)  DA  —15. 

IF  (DR. GT. 15. )  DR  ■  15. 

IF  ( DR.  LT.  - 15.  )  DR  —15. 

CALL  RNGKTA  (NDIM.  T.  DT.  X,  OX .  ER) 

ICNT  »  ICNT  ♦  1 

IF  (ICNT.NE.NSPP)  BO  TO  300 

NPTS  »  NPTS  +  1 

ICNT  -  0 

IF  (IX. GT. 12)  SO  TO  320 
XP ( 1 , NPTS)  ■  X(IX) 

GO  TO  340 

320  IF  (IX. EQ. 13)  XP(l.NPTS)  ■  DE 
IF  ( IX. EQ. 14)  XP(1,NPTS)  ■  BHP 
IF  (IX. EO. 15)  XP ( 1 . NPTS)  ■  DA 
IF  (IX.EQ.16)  XP ( 1 , NPTS)  ■  DR 
340  TP(l.NPTS)  ■  T 

IF  (NPTS. SE. 250)  80  TO  400 
IF  (T.GE.TMAX)  60  TO  400 
GO  TO  300 


CPLOT  PLOTTING  ROUTINE  CALL 


400  CALL  CPLOT  (NPTS . I . TP , XP.PLTSPC , XTITLE, YTITLE , LABEL, FLINE, 
1  SLINE.TLINE) 

STOP 

END 


SUBROUTINE  RNGKTA  WRITTEN  BY  DR  D  K  SCHMIDT,  AtAE  DEPT,  PURDUE  UNIV 


SUBROUTINE  RNGKTA (ND IH ,T ,DT, X ,DX ,ERR> 

REAL  X (NDIM) ,DX (NDIM) , ERR (NDIM) . A (4) , B ( 4 ) ,C<4> 

DATA  A, B. C/0. 3, 0.29289321 9. 1.707106781,0.  16666667, 2. , 1 . , 1 
B  2. .0.5,0.292893219.1.707106781,0.5/ 

DT2  •  DT/2. 

CALL  DXDT  (X.DX) 

J  -  1 

101  DO  102  I  ■  1 .NDIM 

R1  -  DT  t  DX ( I ) 

R2  ■  A ( J )  «  <R1  -  B  <  J )  •  ERR ( I ) > 

X ( I )  ■  X ( I )  ♦  R2 
R2  »  3.  *  R2 

ERR ( I )  »  ERR ( I )  ♦  R2  -  C(J)  •  Rt 

102  CONTINUE 

IF  (J.EQ.4)  RETURN 
J  ■  J+l 

IF  (J.NE.3)  T  ■  T  ♦  DT2 
CALL  DXDT(X.DX) 

GO  TO  101 
END 
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11300 

C 

Ml 

11400 

C 

SUBROUTINE  DXDT,  DETERMINES  DX  FROM  X 

11500 

C 

•*»#**#*****»*********t#*t#*»#*t**t*#*t**»»t*t«t»*****«tt**ttt*tt 

11600 

SUBROUTINE 

DXDT  <  X , DX ) 

11700 

DIMENSION 

X  ( 1 2) , DX ( 1 2 ) 

11800 

COMMON  / BL K 1 2/  CD0,AKD,CL0,CLA, CLAD, CLQ,CLDE,CMO,CMA,CMAD, CMS, 

11900 

l 

CMDE,CYB,CYR,CYP,CYDR,CLB,CLP,CLR,CLDA,CLDR.CNB 

12000 

2 

CNP,CNR,CNDA,CNDR,S,C,B,ANP,A1X,AIZ,AIXZ,AM,AIY 

12030 

3 

DA.OE.DR.BHP 

12100 

C»*#*####FOLLOWlNG 

IS  A  LIST  OF  VARIABLE  NAMES  AND  USES******««t#t# 

12200 

##«**#«##**# 

12300 

c 

V 

8 

SCALAR  VALUE  OF  VELOCITY  VECTOR 

12400 

c 

ALPHA 

■ 

ANGLE  OF  ATTACK.  RADIANS 

12500 

c 

BETA 

8 

SIDESLIP  ANGLE.  RADIANS 

12600 

c 

V2 

3 

TWO  TIMES  THE  VELOCITY 

12700 

c 

VND1 

8 

ALPHADOT  TIMES  MEAN  CHORD  DIVIDED  BY  V2 

12800 

c 

VND2 

a 

Q  TIMES  MEAN  CHORD  DIVIDED  BY  V2 

12900 

c 

VND3 

a 

R  TIMES  SPAN  DIVIDED  BY  V2 

13000 

c 

VND4 

a 

P  TIMES  SPAN  DIVIDED  BY  V2 

13100 

c 

CL 

a 

COEFFICIENT  OF  LIFT 

13200 

c 

CD 

a 

COEFFICIENT  OF  DRAG 

13300 

c 

CY 

■ 

SIDE  FORCE  COEFFICIENT 

13400 

c 

CLL 

a 

ROLLING  MOMENT  COEFFICIENT 

13500 

c 

CM 

a 

PITCHING  MOMENT  COEFFICIENT 

13600 

c 

CN 

8 

YAWING  MOMENT  COEFFICIENT 

13650 

c 

RHO 

8 

AIR  DENSITY.  SLG/FT**3 

13700 

c 

Q 

a 

DYNAMIC  PRESSURE 

13800 

c 

XLIFT 

a 

LIFT 

13900 

c 

DRAG 

a 

DRAG 

14000 

c 

SIDEF 

a 

SIDEFORCE 

14100 

c 

ROLMOM 

a 

ROLLING  MOMENT 

14200 

c 

PITMOM 

a 

PITCHING  MOMENT 

14300 

c 

YAWMOM 

a 

YAWING  MOMENT 

14400 

c 

THRUST 

a 

THRUST 

14300 

c 

SA 

a 

SIN  ALPHA 

14600 

c 

CA 

a 

COS  ALPHA 

14700 

c 

SB 

a 

SIN  BETA 

14800 

c 

CB 

a 

COS  BETA 

14900 

c 

STHE 

a 

SIN  THETA 

15000 

c 

CTHE 

a 

COS  THETA 

13100 

c 

SPHI 

a 

SIN  PHI 

15200 

c 

CPHI 

a 

COS  PHI 

15300 

c 

TTHE 

a 

TAN  THETA 

13400 

c 

SP9I 

a 

SIN  PSI 

13500 

c 

CPSI 

a 

COS  PSI 

13600 

c* 

»•« 

13700 

V  -  SORT (X(1)*X(1)+X(2)*X(2)+X(3)*X(3)) 

13800 

ALPHA  ■  ATAN < X (3) /X ( 1 ) ) 

13900 

BETA  ■  A S I N ( X (2) /V) 

16000 

DALPHA  »  DX (3) / X ( 1 ) 

16050 

N  ■  AM*32, 

2 

16100 

V2  ■  2#V 

16200 

VND1  ■  DALPHA*C/V2 

16300 

VND2  ■  X  (3) *C/V2 

16400 

VND3  ■  X  (6) #8/V2 

16300 

VND4  «  X(4>*B/V2 

16600 

ct 

►  #  « 

16700  CL  ■  CL0+CLA*ALPHA+CLAD*VND1+CLQ*VND2*CLDE»DE 

16800  CD  »  CDO+AKD* (CL*CL) 

16900  CY  ■  CY8*BETA*CYR»VND3tCYP*VND4+CYDR*DR 
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17000 

17100 

17200 

17300  C»»» 
17350 
17400 
17450 
17500 
17600 
17700 
17800 
17900 
18000 
18100 

18200  C*o 
18300 
18400 
18500 
18600 
18700 
18800 
18900 
19000 
19100 
19200 
19300 

19400  C### 
19300 
19600 
19700 
19800 
19900 
20000 
20100 
20200 
20300 
20400 
20500 
20600 
20700 
20800 
20900 
21000 
21100 
21200 
21300 
21400 


CLL  ■  CLB«BETA*CLP*VND4*CLR*VND3+CLDA*DA+CLDR*DR 
CM  *  CM0*CMA#ALPHA+CMAD*VNDl+CHQ«VND2+CHDE*DE 
CN  ■  CNB*BETA+CNP»VND4+CNR*VND3+CNDA*DA+CNDR»DR 


RHO  ■  0.002377  #  (1.  -  6.87E-06  *  X<12)  >  «*  4.2559 
Q  ■  . 5#RH0* ( V«V) 

IF  <CL. 8T. 1 . 53)  CL  ■  .1»<1.0/CL> 

XLIFT  ■  CL*Q»S 
DRA6  ■  CD*Q*S 
SIDEF  «  CY*Q*S 
ROLMOM  ■  CLL*Q*S*B 
P1TM0H  ■  CM*fl*S*C 
YAWMOM  -  CN*Q*S#B 
THRUST  ■  ANP*8HP#530/V 


SA  >  SIN (ALPHA) 
CA  «  COS (ALPHA) 
SB  ■  SIN (BETA) 

CB  >  COS (BETA) 
STHE  «  SIN ( X (8) ) 
CTHE  >  COS (X (8) ) 
SPHI  ■  SIN(X  (9) ) 
CPHI  ■  COS < X (9) ) 
TTHE  »  STHE/CTHE 
SPSI  ■  SIN(X (7) ) 
CPSI  •  COS(X (7) ) 


DX(1)  -  ( XLIFT*SA-DRAB*CB»CA-SIDEF*SB*CA-N«STHE+THRUST) /AH 
1  -X (3) *X ( 3) +X (6) *X (2) 

DX (2)  «  (-DRA8#SB+SIDEF#CB+M*CTHE»SPHI) /AM  -X(6)#X(1) 

1  +X (4) *X (3) 

DX (3)  ■  <-XLIFT*CA-DRA8*CB*8A-SIDEF#SB#SA*N#CTHE*CPHI>/ 

1  AM  -  X (4 ) #X (2) +X (5 ) »X ( 1 ) 

DX (4)  ■  (ROLMOM* (AIY-AIZ) *X(5)*X(6)+AIXZ*(DX(6)+X(4)tX(3)))/AIX 
DX (5)  >  (PITHOM- < AIX-AIZ) *X (4)*X(6)“AIXZ*(X(4)*X(4)-X(6)*X(4))I 
1  /  A I Y 

DX (6)  ■  ( YANHQH* < AI X-AI Y) *X (4)*X(5)*AIXZ*(DX(4)-X(S)«X(6)))/AIZ 

DX (7)  ■  ( X (5) *SPHI+X (6) *CPHI ) /CTHE 

DX (8)  ■  X (5) *CPHI-X (6) *8PHI 

DX (9)  ■  X (4 ) +X (3) *SPHI*TTHE  +  X (6 > *CPHI*TTHE 

DX ( 10)  ■  CPSI*CTHE*X ( 1 ) * ( -SPSI *CPHI*CPBI*8THE*SPHI ) *X (2) ♦ 

1  (SPSI*SPHI*CPSI *STHE*CPHI ) *X (3) 

DX(ll)  ■  SPSI*CTHE*X(1)*(CP9I*CPNI*SPSI*STHE»SPHI)«X(2>* 

1  (-CPSI*SPHI+SPSI*STHE*CPHI)*X<3) 

DX (12)  ■  STHE* X  < 1 ) -CTHE*SPHI*X (2) -CTHE*CPHI*X (3) 

RETURN 
END 


Data  File  "CNST1JT" 


)  66 
Id. 


»*.12»i.U*1.66*i‘.16..0L5*.OC93*-.29,-i.O*“11.0»* 

-.CG4*.00235»-0.0?£**.3€*.096*.025*.00105*.17»“. 

00!<i4»l5.65*l.'73<>*f.72,.70,l.*,9.3*0.0*.9IlT*6.S 


Data  File  "SYSDAT" 


0.  *09. 5*C»0»0.0*O.CCC*. 130*0. 00 *0.0*7 500. 

*  1.11*  0.*  0. 
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APPENDIX  E.  System  Poles  for  Root  Loci  Plots 


Longitudinal  Poles 


?  HD  =  0  0E5  m  SEC 
i'H  =  I.  ct  5EC  -'FT 
KV  -  BHf’  :  "7  EEC 

i  5**5  ♦  9.630*2  S**4  *  20.824344  5**3  *  5.4'2'2025  S**2  *  2.28162392  5  *0 
POLES  ARE: 

-4,19655335  l  -1.04669128 

-4.19655835  i  1.0466*128 

-.118901651  i  -.3283/616 

-.11890165!  .)  .  32317616 

0  J  0 


*  HD  =  -.!  DCS  '  FT.  SEC 
KH  =  .2  FT.  SEC  /FT 
kV  =  ,1  8HP  '  FT. 'SEC 
1  3* *5  *  9.15366437  S**4  *  25. 

POLES  ARE: 

-.508715358  1  -.330638958 

-.508715358  l  .330638958 

-3.351*8519  i  0 

-4.412666*1  I  0 

-.17208205  JO 


►Hi  =  -, 15  DEG  ■  FT 'SEC 
■-H  =  .2  ft  SEC  'FT 
=  .1  &HP  r  T  SEC 
i  S**5  ♦  ». !32'1715s  5**4  *  25. 

FILES  APE: 

:.  ;"6l ■  bS-  I  0 

-4,  ’92*5812  > 

-.6l’»644!4  I  -.45',u62806 

- . 643864414  :  .45’9i2s06 

-.i'5»7'*;2  i  v 


1HD  =  -.2  DEG  ■  FT/SEC 
KH  =  .2  Ff  SEC  FT 
TV  =  , 1  SHF  FT/SEC 
1  S**5  *  9. 1  1  187873  5**4  ♦  25 


4759382  S**3  *  22.7784477  5**2  *  8.90190019  5  *  .979279342 


4378578  5**3  *  26.1930904  5**2  *  12.2120383  S  ♦  1.46793395 


.3997775  5**3  ♦  29.6077331  S»*2  ♦  15.5221765  S  +  1.95593185 


POLES  ARE: 
-2.00669996  1  0 
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-5. 23439389  t  0 

-.946614841  j  -.57618756 

-.846614841  l  .57818756 

-.I7'555202  j  0 


KHD  =  -.25  DEG  !  FT/SEC 
KH  =  .2  FT/SEC  /FT 
KV  =  .1  BHP  /  FT/SEC 

1  S**5  ♦  9.09098592  S»*4  ♦  25.3616971  S«3  ♦  33.0223758  S«*2  +  18.8323146  S  +  2.44327303 

POLES  ARE: 

-5.4207603  i  0 
-1.13625785  j  -.883970426 

-1.13625785  )  .883970426 

-1.21935259  j  0 

-.178357333  j  0 


KHD  =  -.3  DEG  /  FT/SEC 
1H=  .2  FT/SEC  /FT 
(TV  =  .1  8HP  /  FT/SEC 

1  S*«5  ♦  9.0700931  S»#4  ♦  25.3236168  S«3  ♦  36.4370185  S»»2  ♦  22.1424527  S  ♦  2.92995751 
POLES  ARE: 

-1.15545131  )  -1.26088437 

-1.15545131  j  1.260B8437 

-5.57560036  j  0 

-1.004782  j  0 
-.178808123  jO 


KHD  =  -.2  DEG  /  FT/SEC 
KH  =  .1  FT/SEC  /FT 
KV  -  .1  BHP  /  FT/SEC 

1  S*»5  ♦  9.  11187873  S«*4  *  25.4081346  S»3  ♦  29.6212938  S»«2  ♦  14.1536072  S  ♦  .977965923 

POLES  ARE: 

-.0823008552  )  0 

-.994005835  j  -.525117323 

-.994005835  j  .525117323 

-1.79062917  i  0 

-5.25093704  j  0 


KHD  =  -.2  DEG  /  FT/SEC 
KH  *  .15  F  T  -  SEC  /FT 
KV  =  .1  BHP  /  FT/SEC 

1  S**5  ♦  9.11107873  S**4  ♦  25.403956  S»»3  +  29.6145134  S*»2  ♦  14.8370918  S  ♦  1.46694880 
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-.128241961  )  0 

-.915390605  j  -.551649739 

-.915390605  )  .551649739 

-1.91013555  j  0 

-5.24272001  j  0 


KHD  =  -.2  DEG  /  FT/ SEC 
KH  -  .2  FI/ SEC  /FT 
KV  =  .1  BHP  /  FT/SEC 

1  5**5  ♦  9.11187873  S**4  ♦  25.3997775  S**3  +  29.6077331  S**2  +  15.5221765  S  ♦  1.95593185 

POLES  ARE: 

-2.006*9996  j  0 

-5.23439389  j  0 

-.846614841  j  -.57618756 

-.846614841  j  .57618756 

-.177555202  j  0 


KHD  =  -.2  DEG  /  FT/SEC 
KH  =  .25  FT/SEC  /FT 
KV  =  .1  BHP  '  FT/SEC 

1  S**5  ♦  9.  11187873  S**4  ♦  25.3955989  S**3  ♦  29.6009527  S**2  ♦  16.2064611  S  ♦  2.44491481 

POLES  ARE: 

-2.08932918  j  0 

-5.22595506  j  0 

-.783362962  )  -.600383671 

-.783362962  j  .600383671 

-.229868565  j  0 


KHD  =  -.2  DEG  /  FT/SEC 
KH  -  .3  FT/SEC  /FT 
KV  =  . I  BHP  /  FT/SEC 

1  3**5  ♦  9.11187873  S**4  ♦  25.3914204  S»»3  ♦  29.5941724  S**2  ♦  16.B907457  S  ♦  2.93389777 

POLES  ARE: 

-2. 16250743  j  0 

-5.21739973  j  0 

-.724015404  i  -.625789061 

-.724015404  j  .625789061 

-.293’40771  j  0 


KHD  =  -.2  DEG  /  FT/SEC 
KH  =  .2  Ff/SEC  /FT 
KV  *  .04  BHP  '  FT/SEC 

1  S»*5  ♦  8.77316074  S*»4  ♦  22.5915516  S«3  ♦  23.3513994  S*»2  ♦  10.8772586  S  ♦  1.02598105 
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POLES  ARE: 

-.1232)3703  )  0 

-.720433944  j  -.535504229 

-.720433944  )  .535504228 

-1.97301432  i  0 

-5.23544482  j  0 


‘HE  =  -.2  DE6  /  FT /SEC 
1H  =  .2  FT/SEC  /FT 
RV  =  .04  BHP  /  FT/SEC 

1  3**5  ‘  9.88404473  S«4  ♦  23.5274249  S**3  ♦  25.434844  S«2  4  12.4255444  S  4  1.33594445 

POLES  ARE: 

-.145187511  ]  0 

-.741329041  j  -.553845598 

-.'41329041  j  .553845598 

-1.98295914  j  0 

-5.23524197  j  0 


KHD  =  -.2  DEE  ■  FT/SEC 
KM  =  .2  FT/SEC  /'FT 
KV  =  ,08  BHP  /  FT/SEC 

1  5*45  ♦  8.99897273  S*»4  4  24.4437022  S»«3  4  27.5222885  S**2  4  13.9738705  S  4  1.44594825 

POLES  APE: 

-1.9941143  )  0 

-5.23483885  j  0 

-.803528034  i  -.54734352 

-.807528034  ,  .54734352 

-.142941521  j  0 


KHD  =  -.2  DEG  /  FT/SEC 
KH  =  .2  FT/SEC  /FT 
KV  =  . 1  BHP  /  FT/SEC 

1  S**5  4  9.  11107873  S*»4  4  25.3997775  S**3  ♦  29.4077331  S«2  4  15.5221745  S  4  (.95593185 

POLES  ARE: 

-2.00449997  ]  0 

-5.23439388  j  0 

-.846414839  j  -.576187561 

-.844414839  j  .576187541 

-.  177555202  )0 


*HD  =  -.2  DEG  /  FT/SEC 
KH  =  .2  FT/SEC  /FT 
KV  =  .12  BHP  /  FT/SEC 

1  5*45  4  9.22478473  S444  4  26.3358527  S443  4  31.6931776  S442  4  17.0704824  S  4  2.26591545 
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-0l£5  -PE: 

-.390089733 
-5. 27392534 
-2.02097073 
-.l89',09889 


>  -.580455401 
l  .580455401 
)  0 
j  0 
)  0 


K.HD  =  -.2  DE6  /  FT/SEC 
PH  =  .2  FT/SEC  /FT 
TV  =  .14  BNP  /  FT/SEC 

1  5**5  +  9.33769073  S**4  *  27.271928  S**3  *  33.7784222  S*»2  ♦  18.4187883  S  + 


POLES  ARE: 

-.933523349  i  -.580293241 

-.953523349  j  .580293241 

-5.2334313  j  0 
-2.03724597  <  0 


-.199964744  |0 


KhB  ®  -.2  DE6  !  FT  SEC 
PH  =  .2  FT/ SEC  /FT 
KV  =  .16  BHP  /  FT/SEC 

1  S**5  ♦  9.45059474  S**4  +  28.2080033  S**3  *  35.8440448  S**2  *  20.1470943  S 

FOLES  ARE: 

-5.23290958  j  0 

-.976525323  )  -.575809024 

-.976525323  j  .575809024 

-2.05591145  ]  0 


2.57589904 


♦  2.B8588265 


USAFA-TR-86-2 


POLE?  ARE: 

-0.64578947  )  0 

-1.64450809  i  -J. 95216212 

-1.64450809  j  3.95216212 

-.101569701  !  -.179452201 

-.101569701  )  .179452201 


m  =  30  DE6/0 
kR  =  1  DE6/  RAD/S 
KPSI  =  .1  RAD/S  /RAO 
KVD  =  -20  DEG/  RAO/S 

1  S**5  *  12.5275614  S**4  ♦  54.4966084  S*«3  ♦  189.418558  S**2  ♦  38.7851718  S  ♦  7.51334065 

POLES  ARE: 

-8.71023749  j  0 

-1.80600033  j  -4.14745607 

-1.80600033  ]  4.14745607 

-.10266164  j  -.177803379 

-.10266164  j  .177803379 


KH¥  =  40  DEG.'o 
KR  =  1  DEG/  RAD'S 
kPSI  *  .1  RAD/S  /RAD 
m  -  -20  DEG/  RAD/S 

1  3**5  ♦  12.9819845  S**4  ♦  60.6427901  S**3  *  213.897442  S**2  *  44.1676258  S  *  8.419B2831 

POLES  ARE: 

-8.78719367  j  0 

-1.99373372  j  -4.35188056 

-1.99373372  J  4.35188056 

-.10366168  )  -.176271039 

-.10366168  j  .176271039 


W  -  50  DEG; a 
KR  =  1  DEG/  RAD/S 
KPSI  *  .1  RAD/S  /RAD 
KVD  =  -20  DEG/  RAD/S 

1  S»»5  *  13.5188503  S«»4  ♦  67.9042413  S**3  ♦  242.818953  S**2  ♦  50.5271121  S  +  9.49083428 

POLES  ARE: 

-8.88083769  j  0 

-2.21442524  j  -4.56544296 

-2.21442524  j  4.56544296 

-.104581081  j  -.174843122 

-.104581081  j  .174843122 


KNY  =  60  DES /g 
KR  *  1  DEG'  RAO/S 


USAFA-TR-85-2 


1PSI  =  :  I  »f!0 •  g  /RAD 
I  'D  =  -:v  DEE;  RAD/S 

1  5**5  *  14.1628327  S**4  *  76.6148123  S**3  +  277.513218  S**2  ♦  58.1562149  S  +  10.775618 

POLES  ARE: 

-8.99750778  i  0 

-2.47723309  j  -4.78718744 

-2.47723309  )  4.78718744 

-.105429385  i  -.573509154 

-.105429385  i  .173509154 


‘O-  =  30  DEC/ a 
i.P  =  ,c  DEC"  RAD/S 
*  PCI  =  .1  RAD/S  .'RAD 
KfD  =  -20  DEC/  RAD/S 

i  5**5  +  12.9555614  S**4  *  55.9238229  5**3  +  190.287381  S**2  ♦  8.80488741  5  ♦  4.50800439 

POLES  ARE: 

-9,14473222  i  0 

-1.90550947  t  -4.11124422 

-1.70550947  i  4.11124422 

-.0199051288  i  -.153661036 

-.01 W51288  i  .  153651036 


IN'i  =  30  DES/a 
IP  =  .6  DE8/  PAD/5 
KPSI  *  .1  RAD. 3  /RAD 
HD  =  -20  DEG/  PAD/S 

i  5**5  *  12.7615614  5**4  *  55.2102156  S**3  +  189.85297  S**2  *  23.7950296  S  *  6.01067252 

POLES  ARE: 

-8.92673673  i  0 

-1.85720673  j  -4.12897346 

-1.85720673  i  4.12897346 

-.0602056227  i  -.170952076 

-.0602056227  t  .170952076 


Mir  =  30  DEQ/o 
1R  =  1  DEG/  RAD/S 
IPS!  =  .!  RAO/S  PAD 
=  -20  DES/  RAD/S 

1  5**5  ♦  12.5275614  5**4  ♦  54.4966084  5**3  ♦  189.418558  S**2  ♦  38.7851718  S  *  7.51334065 

POLES  ARE: 

-8.71023749  j  0 

-1.80600033  i  -4.14745607 

-1.80600033  i  4.14745607 

-.10246164  i  -. 177803379 


USAFA-TR-85-2 


10266 164  i  .  177803379 


KNi  =  30  DEG/a 
KR  =  1.2  DEG/  RAD/S 
KPSI  =  .1  RAD/S  /RAD 
K/0  =  -20  DEG/  RAD/S 

i  5**5  *  12.2935614  S**4  +  53.7330012  S**3  +  188.984146  S»*2  ♦  53.775314  S  ♦  9.01600878 

POLES  ARE: 

-8.49525597  j  0 

-1.75174853  j  -4.16677264 

-1.75174853  )  4.16677264 

-.1474042  )  -.173834725 

-.1474042  j  .173834725 


KNV  =  30  DEG/ a 
KR  -  1.4  DEG/  RAD/S 
KPSI  =  . 1  RAD/S  /RAD 
KYD  =  -20  DEG/  RAD/S 

1  S**5  ♦  12.0595614  S**4  +  53.069394  S»»3  ♦  188.549734  S**2  ♦  68.7654562  S  ♦  10.5186769 

POLES  ARE: 

-8.28179457  j  0 

-1.69432034  j  -4.18701812 

-1.69432034  j  4.18701812 

-.194563085  j  -.15620246 

-.194563085  i  .15620246 


kN(  =  30  DE6/g 
KR  =  1  DEG/  RAD/S 
KPSI  =  .06  RAD/S  /RAD 
k*D  =  -20  DE6/  RAD/S 

1  5**5  ♦  12.5275614  S*»4  ♦  54.5434084  S**3  ♦  189.556599  S»»2  ♦  38.85B25  S  ♦  4.50800439 

POLES  ARE: 

-8.70522368  )  0 

-1.90588843  j  -4.14995057 

-1.8058B843  j  4.14995057 

-.105280445  j  -.119153661 

-.105280445  j  .119153661 


kKt  =  30  DEG/g 
KR  =  1  DEG/  RAD/S 
KPSI  *  .08  RAD/S  /RAD 
KYD  *  -20  DEG /  RAD/S 

1  5**5  ♦  12.5275614  S**4  ♦  54.5200084  S»*3  ♦  189.487579  S**2  *  38.8217109  S  ♦  6.01067252 


USAFA-TR-85-2 


POLES  APE: 

-8.7077312  i  0 
-1.80594354  j  -4.14870324 

-1.80594354  j  4.14870324 

-.103971572  j  -.151346198 

-.103971572  )  .151346198 


KNV  =  30  DES/g 
KR  *  1  DES/  RAO/S 
KPSI  =  .1  RAD/S  /RAD 
KVD  =  -20  DES/  PAD/S 

1  S«5  ♦  12.5275614  S*»4  +  54.4966084  S»*3  *  189.418558  S*«2  ♦  38.785171B  S  ♦  7.51334065 

POLES  ARE: 

-8.71023749  i  0 

-1.80600033  ]  -4.14745607 

-1.80600033  t  4.14745607 

-.10266164  j  -.177803379 

-.10266164  ]  .177803379 


KNV  =  30  DES '0 
KR  =  1  DES/  RAD/S 
KPSI  *  .12  RAD/3  /RAD 
K»P  =  -20  DES'  RAD/S 

1  S*»S  ♦  12.5275614  S»»4  ♦  54.4732084  S«»3  ♦  189.349537  S«#2  6  38.7486326  S  ♦  9.01600878 

POLES  ARE: 

-8.71274251  i  0 

-1.80605881  j  -4.14620905 

-1.80605881  ]  4.14620905 

-.101350644  i  -.200805122 

-.101350644  j  .200805122 


KNV  =  30  DES/a 
KR  =  1  DES/  RAD/S 
KPSI  =  .14  RAD/S  /RAD 
KVD  -  -20  DEB/  RAD/S 

1  S«5  ♦  12.5275614  SM  ♦  54.4498084  S«3  ♦  189.280516  S»#2  ♦  38.7120935  S  ♦  10.5186769 


,V 


POLES  ARE: 

-8.71524628  )  0 

-1.80611899  j  -4.14496219 

-1.80611899  )  4.14496219 

-.100038582  )  -.22143087 

-.100038582  )  .22143087 


KNV  =  30  DEG/ q 


mmsmm 


'y'./vv/v.'  * 

■ -■  '•> y&y-  v 


'•  «"»  V  *  *  N 


mm 


USAFA-TR-85-2 


=  i  DE'j'-  ^«D/H 
kf'St  =  .i  PAt>  S  -PAD 
1  ‘  L*  =  .'  D6B/  *  •*  0 

i  3**5  t  8. 6t437:45  5**4  ♦  29.0232275  5**3  ♦  184.017261  9*»2  ♦  60.6243416  £  ♦  7.51334065 

POLES  APE : 

-7.S33*2fc<*  -  v 

-.244°83415  i  -4. 71o75074 

244*86415  1  4. ’1675079 

-,l’0236/65  >  -.118374156 

-.170256765  1  .116374156 


>&»  *  30  £ES/c 
t?  -  1  EE&;  ’■AO/ 5 
i fSI  *  .1  RAD/S  RAD 
t  5  =  -U«  DES/  RAC'S 

;  4»*5  »  U. 595966*  9**4  +  41.759917?  S«3  *  186.71791  S**2  ♦  49. ’047567  5  *  7.51334065 

ROLES  ARE: 

-8. 1915929  i  C 
-1.06527569  1  -4.5061798’ 

-1.06527563  ]  4,50617987 

-.13676133°  i  -.155157179 

-.13676133°  j  .155157179 


»«»  r  70  CES/3 
=  1  DEC;  RAC'S 
k°3!  =  .1  RAO/S  /RAD 
>'&  =  -20  DEC-'  PAO/S 

1  2**5  ♦  12.5275614  5**4  ♦  54.4966084  S**3  *  189.418558  S**2  ♦  38.7851718  S  ♦  7.51334065 

ROLES  APE: 

-8.'lo23749  j  0 

-1.3' ’60003!  1  -4.14’45607 

-1.80600037  j  4.14745607 

-.10266164  1  -. 177303379 

-.■0:64164  i  . 17’8033’9 


t  N  ■  =  30  DEC  1 
’R  -  1  Dca'  RAD  5 
/PS!  -  .!  RPC'S  'RAD 
!  :  -  -3"  LEG ■  RAD <5 

1  it*!  t  14,4591559  9**4  ♦  67,2772988  S**3  ♦  192.119206  S**2  *  27.8655869  S  ♦  7.51334065 
ROLES  A9£; 

-6>leCRgC3e  , 

-2.42565318  1  -3.65940816 

-1. 42565618  ;  3.65°4''8l6 

-.  7689’7i.j4«  ;  -,190’97648 

•.0o897’t043  ;  .1°>'7°7648 


>•2 


USAFA-TR-85-2 


•:N;  =  :0  jEG’j 
-  I  0 lG-’  PAD  c 
IPS!  *  .1  Rr.5'5  /PSD 
MD  -  -4'.'  DEG’  RAD/S 

i  S»*5  +  16.3907504  3**4  *  ?9.969°8°3  S**I  +  194.81*855  S*»2  *  16.946002  S  *  7.51334065 


POLES  APE: 

-10.5453974 

-2.90603603 

-2.88603603 

-.0366404783 

-.036640478! 


!  0 

i  -3.08809945 
]  3.03809945 
;  -.196309277 
:  ,196I0»277 


it  U.t.  QOVIANMEMT  PflIMTINO  OFEICC:1«K-  STKXO  /  21074 


243 


'  •  .  4  .  «  .  «V  *  ■  .  ■  .  ,•  ,*  .  * 


1.  Please  discontinue  the  Digest. 

2.  Change  my  address  to:  (print  new  address) 


END 


FILMED 


9-85 


DTIC 


